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Introduction 


The work presented in this thesis concerns the optics research field that is commonly 
called micro or nano-photonics. The term “photonics” emerged at the end of last cen- 
tury to describe the fundamental study of light and its applications to technology, 
whose purpose is to accomplish functions traditionally undertaken by electronics, 
such as telecommunications, information processing, and so on. In fact, in the effort 
to further progress the high density transistor integration and thus the modern elec- 
tronic system performances, scientists are now turning to photons instead of electrons 
as the information carriers in communication. Photons show several advantages. Light 
can carry a larger amount of information per second. In addiction, at the single photon 
level, quantum information networks and quantum cryptography scheme can be pur- 
sued. Furthermore, photons are not as strongly interacting as electrons, thus reducing 
energy losses. However, in spite of the numerous advantages of photons, all optical 
integrated circuits have yet to be commercially available on large scale due to the 
difficulties in designing efficient optical components building blocks. 

Nowadays, the suited objects to advance this research field are made up of pho- 
tonic crystal devices, because of their versatility and ability to control the propagation 
of light at the nanoscale. Since the introduction of the photonic crystals, a great deal 
of groundbreaking phenomena concerning light-matter interaction at optical frequen- 
cies have been theoretically predicted and experimentally verified. They led to large 
technological developments, especially giving rise to innovative optoelectronic and 
photovoltaic applications. For instance, microscopic light manipulation is achievable 
through photonic crystals that exhibit finite energy band-gap for light. In particular, 
two-dimensional photonic crystal cavities fabricated on semiconductor slab are the 
state of the art devices to strongly localize electromagnetic (e.m.) fields in volumes 
below a cubic optical wavelength of light, thus acting as nano-resonators. In fact, they 
provide an efficient mechanism for controlling and inhibiting spontaneous light emis- 
sion, or for reducing almost to zero the light group velocity. A noteworthy advantage 
of two-dimensional photonic crystal cavities is that they can be included in all-optical 
integrated circuits fabrication. The high quality factors enable these systems to reach 
the strong coupling regime between photonic localized modes and two level quantum 
emitters. However, the fundamental request to maximize these effects is the accurate 
positioning of the light source in the spot where the localized mode is maximum. 
Moreover, when systems of interacting resonators are conceived, a basic condition to 
be satisfied is the precise control on the overlap between adjacent elements to induce 
an effective coupling. The sub-wavelength imaging of the (e.m.) fields localized in 
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optical nano-resonators is therefore of the utmost relevance in the road map to inte- 
grate photonic nano-resonators in chipset architectures. In order to achieve this goal 
we need to overcome the resolution limitations that are imposed by light diffraction 
to any standard microscope. 

In Chapter (1) we review the more commonly methods employed to achieve a sub- 
wavelength spatial resolution. In particular, we describe the principle of operation of 
the scanning near-field optical microscope. 

In Chapter (2) we describe the basic properties of light confinement in photonic 
crystal nanocavities and also in optical nano-resonators based on disordered arrange- 
ments of light scatterers, where multiple scattering regime gives rise to randomly 
placed nanocavities. In particular, we illustrate how near-field techniques can experi- 
mentally probe the e.m. optical fields localized at the nanoscale and how finite ele- 
ments numerical calculations can simulate the behavior of light in such nano-resona- 
tors. 

In Chapter (3), by exploiting the scanning near-field probe perturbation imaging 
technique, we report the sub-wavelength mapping of the electric or magnetic localized 
field component of light using dielectric tapered probes or aperture metal coated 
probes, respectively. The advent of artificial metamaterials operating at optical fre- 
quencies, in which the magnetic interaction with light can be as relevant as the electric 
response, makes it straightforward to simultaneously detect both the e.m. field com- 
ponents. We develop a plasmonic-based near-field probe to strongly enhance the col- 
lection efficiency and, more interestingly, to achieve an ultra-bright and sub-wave- 
length simultaneous detection of both the resonant electric and magnetic fields. 

Photoluminescence based imaging methods, however, require optically active 
samples with the constrain of spectral and spatial matching between photonic mode 
and the light emitters, which may suffer from bleaching or blinking. Therefore, a pure 
optical method that can be applied on any kind of high quality factor nano-resonators 
to retrieve the confined modes distributions is actually missing. In order to achieve 
this goal we investigate the localized nanophotonic modes by developing a different 
approach. The experimental method, presented in Chapter (4), combines scanning 
near-field optical microscopy with resonant scattering spectroscopy, and it is called 
Fano-imaging. This technique largely extends the investigation of nanoscale localized 
light states, since it is applicable to nano-resonators based on any kind of material, 
even where light sources cannot be embedded. Moreover, resonant scattering experi- 
ments exhibit spectral Fano resonances, which correspond to the interference between 
light directly scattered from the sample and light scattered after being resonantly cou- 
pled to the localized mode. From the detailed analysis of Fano lineshapes it is possible 
to retrieve a deep sub-wavelength imaging of both the electric field intensity and the 
electric field spatial phase distribution, polarization resolved. Thus, we obtain unprec- 
edented local information about the resonant light states. 

In Chapter (5) we deeply investigate, both theoretically and experimentally, sys- 
tems composed by coupled nano-resonators, called “photonic molecules’. In fact, 
light behavior in system based on multiple aligned photonic crystal nanocavities re- 
sembles the molecular interaction where the resulting normal modes exhibit energy 
splitting and spatial delocalization. This condition is achieved by an evanescent pho- 
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ton tunnelling. These structures represent a large research topic also for quantum op- 
tics. However, a fundamental requirement to create proper quantum-optics devices is 
the design and control of adjacent nanocavity modes at the target wavelengths, within 
an accuracy that is not experimentally obtainable due to the fabrication tolerances. 
The compensation of the fluctuations related to the structural disorder can lead to the 
control of the resonance wavelength of each resonator and of the tunnelling coeffi- 
cient. In our analysis, we compare the interaction strength and the mode symmetry 
character of photonic molecules aligned along different lattice symmetry directions 
and composed by two or three photonic crystal nanocavities. Moreover, we theoreti- 
cally evaluate the proper set of parameters to efficiently act on the coupling strength 
at the fabrication level or even with post-growing techniques. In particular, we de- 
velop a laser-assisted local oxidation of the dielectric environment in which the pho- 
tonic cavities are fabricated. This oxidation induces a smooth and irreversible spectral 
shift of the resonant modes confined close to the laser spot. Therefore, the spatial 
selectivity of the post-fabrication technique is exploited not only to adjust the resonant 
wavelength of a given nano-resonator to a target value, but, more strikingly, to modify 
the coupling strength in photonic molecules. Finally, by comparing the case of two 
and three nano-resonators and next-nearest-neighbour coupling in array of photonic 
molecules. 

The last part of the thesis deals with the engineering of light states localized in 
strongly scattering disordered media, as reported in Chapter (6). Light behavior in 
complex disordered systems attracts a lot of attention by fundamental physics as well 
as by technological applications involved in imaging through turbid media such as 
fog, clouds or living tissues. The occurrence of localized states in disordered media is 
a well-established phenomenon traced back to Anderson localization for electrons. 
However, the interaction between adjacent light sates driven by disorder has still to 
be completely understood. In Chapter (6) we demonstrate the possibility to engineer 
the confinement and the mutual interaction of modes in a two-dimensional disordered 
photonic structure. On one hand, the strong light confinement is achieved at the fab- 
rication stage by an optimization of the design parameters. On the other hand, exploit- 
ing the local laser oxidation, we probe the interaction between overlapping localized 
modes, thereby paving the way for the creation of open transmission channels in 
strongly scattering media. 
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Chapter 1 


Optics at the nanoscale 


Until the second half of last century, the spatial resolution achievable by an optical 
microscope, even without any imperfections, was hampered to the diffraction limit. 
Only the development of new kinds of microscopes allowed for investigating object 
and light emitters down to a sub-wavelength scale. Here we report the basic principles 
of nanoscale optics. Firstly, by describing the limit imposed by light diffraction. Af- 
terwards, by theoretically investigating and reviewing the methods largely employed 
to overcome this limitation and to achieve a sub-wavelength spatial resolution. We 
focus on detecting the optical near-field employing scanning microscopy, in particular 
we summarize the possible implementations and we introduce the most commonly 
used near-field probes. 


1.1 The diffraction limit 


The time when humans began to exploit the optical properties of materials to en- 
large the vision of small or distant objects is traced back to centuries ago. As the sci- 
entific method was developed, both optics fundamental laws and technical achieve- 
ments quickly enlarged the understanding of natural phenomena together with the in- 
strument assisted eye perception. Likewise, the first prototype of an optical micro- 
scope (1610) is attributed to Galileo Galilei, however, only since the pioneering de- 
velopment of van Leeuwenhoek (1671) it was possible to observe red blood cells and 
bacteria gaining a few um size resolution [1]. In the 1873 thanks to the advances of 
the theory of light waves Abbe, working at the Zeiss microscope factory, clarified that 
also an ideal microscope without any imperfections in the lenses or misalignments 
would not be able to overcome the physical resolution limit imposed by diffraction 
[2]. In particular, the Abbe’s limit accounts for the smallest radial resolution (Ary,) 
obtainable in the image plane by an optical system that collects monochromatic light 
with vacuum-wavelength equal to A: 


0.614 
n sin 


Min[A4r,,] ~ (1.1) 

where n is the refractive index of the medium between the source and the micro- 
scope, @ is the half-angle of the maximum light cone accepted by the lens, as reported 
in Fig. (1.1 a), and it is related to the numerical aperture of the optical system that is 
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defined as NA = n sin@. The resolution of an optical microscope evaluated in 
Eq. (1.1) is defined as the shortest distance between two distinguishable points on a 
specimen that can still be observed as separate entities, following the Rayleigh crite- 
rion [3, 4]. In order to understand the basic imaging process let us consider the image 
formation of a point-like source (whose dimension are D « J), as schematically re- 
ported in Fig. (1.1 a). The incoming light is diffracted by the point-like specimen (S), 
as always happens when light encounters an obstacle, than it is once again diffracted 
by lens L, which mimics a generic microscopy objective. In the case of a finite size 
object, rays leaving with the same angle cross again exactly at the back focal plane 
Èp, and the finite object image appears in the conjugate plane X; [5]. 

Abbe found by experimentation that a larger NA resulted in higher resolution, 
even though the apparent cone of incident light filled only a small portion of the lens. 
In fact, the entire size of L interacts with the diffracted light. Since NA can not be 
large enough to collect all the diffracted light (0 < 77/2), the image recorded in X; 
does not exactly correspond to the original object and it is blurred with the typical 
Airy disk intensity pattern shown in Fig. (1.1 b) [5]. The Airy disk is defined as the 
diffraction pattern created by a circular pupil of diameter D, uniformly illuminated. It 
consists on a central disk of angular radius 1.227/D, surrounded by many alterna- 
tively dark and bright rings [6]. Therefore, in any finite optical system that detects 
e.m. radiation the image of a point-like light source is not a delta-like function but it 
is broadened as an Airy disk, which thus represents the point spread function. Abbe 
used the Airy disk properties to retrieve the net imaging resolution stated in Eq. (1.1), 
that is valid for every detector even with the highest collection efficiency. Moreover, 
any finite specimen can be schematized by a collection of point-like sources whose 
size is smaller than 4, consequently every time we use standard microscopy systems 
to perform images of small objects we inevitably loose information. 

Although modern oil or solid immersion microscopes can achieve large equiva- 
lent NA up to 2.0 in the visible spectral range [7], Eq. (1.1) does not allow resolving 
an optical image with a sub-wavelength resolution (Ax « 4), thus closely resembling 
and being related to the Heisenberg’s uncertainty principle [8]. For instance, in the 
telecom window near-infrared spectral range the largest transparent NA is about 1.4 
and gives a spatial resolution of Ax ~ 700 nm. 


1.1.1 Evanescent waves 


As Maxwell equations were derived (1861), they entirely described the e.m. field 
behavior. But it took almost a century to discover a peculiar optical phenomenon that 
has put fundamental limits to any image reconstruction: this is the existence of e.m. 
evanescent waves in any diffraction process [9-11]. 

Evanescent e.m. waves are exact solutions of Maxwell wave-equations that do not 
propagate in space but exhibit an exponential decay of the amplitude. When total in- 
ternal reflection occurs, they commonly describes light behavior in the low refractive 
index medium, as reported in Fig. (1.2 a). Evanescent waves do not deliver a net en- 
ergy transfer. Only if frustrated total internal reflection takes pace it is possible to 
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generate propagating waves in an high refractive index material close to the evanes- 
cent waves source (closely resembling the electrons tunnel effect), as shown in 
Fig. (1.2 b) [5]. Classical wave optics based on Huygens-Fresnel principle has a fun- 
damental deficiency since it overlooks the role of evanescent waves in diffraction. In 
fact, the classical solution of diffraction problems are accurate in the radiation zone, 
but fails to reproduce what happens in close proximity to the light source. 

In order to analyze the reason why, by collecting only propagating radiation with 
any kind of microscope (even an ideal one with 6 = 7/2), the image is always 
blurred, we exploit the argument based on the angular spectrum decomposition that is 
a noteworthy technique for describing optical fields in homogeneous media. Using 
this Fourier based characterization we can estimate the series expansion of an arbitrary 
field in terms of propagating and evanescent waves with variable amplitudes and di- 
rections. The description of the e.m. field radiated by a light source or scattered by an 
object, as shown in Fig. (1.3), evaluated by the angular spectrum representation in a 
arbitrary z = const. plane is [12]: 


A(r,t) = ff, A (kx ky, 0) exp(ik,z + ikyx + ikyy —iwt)dkydky (1.2) 


where A(r, t) represents either the electric (E) or the magnetic (H) field and A (k) 
its Fourier transform; while w = 2rrncA7* is the angular frequency. Since the wave- 
number is k = /k2 + k2 + k2, the e.m. fields described by Eq. (1.2) have different 
behaviors in two complementary regimes in which the wavenumber component k, is 
either real or imaginary. 


2 

ki+k3<(2) > keR (1.3) 
2 

k? + k? > (=) > kec (1.4) 


Hence, for a purely dielectric medium with no losses, turns out that the factor 
exp(ik,z) is either an oscillatory or an exponentially decaying function if Eq. (1.3) 
or Eq. (1.4) respectively holds. Therefore, the angular spectrum is a superposition of 
both plane and evanescent waves. Remarkably, only if all the spatial frequencies 
(ky, ky) are collected it is possible to reconstruct without any distortions the field 
coming from the specimen by exploiting Eq. (1.2). 

On the other hand, if the imaging system is sufficiently separated from the object, 
the contribution of the evanescent waves is negligible and only the low spatial fre- 
quencies, which satisfy Eq. (1.3), have to be considered. In other words, the image 
collected is a low pass filtered representation of the original field. By filtering out the 
high in plane spatial frequencies described by Eq. (1.4) we loose information about 
the object smallest details. In fact, dealing only with propagating waves the largest 
detectable in plane wavenumber k,, = 2rrnA I sin0, where we consider also the prac- 
tical limitation in collecting a given range of angles, as highlighted in Fig. (1.1 a). 
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Accordingly to the properties of conjugate variables of Fourier optics, which resem- 


bles the Heisenberg principle, the radial in plane uncertain Ar,, = y (4x)? + (Ay)? 
and the field finite angular spectrum Ak, are related by: 


Hence, to achieve the minimum spatial uncertain the largest angular spectrum is 
needed, which corresponds to Ak ,, = 21NAA~, giving for the ideal resolution: 


Min[4r,,]> aa (1.6) 


Therefore, even if an objective may collect 0 = 1/2 it would not be able to inter- 
act with the evanescent waves and the specimen image would consequently still be 
incomplete. This limit is almost four times better with respect to Abbe’s prediction of 
Eq. (1.1). In fact, Abbe required the paraxial approximation (small @) and applied the 
arbitrary Rayleigh’s criterion to the case of two parallel dipoles oriented perpendicular 
to the optical axis. Nevertheless, the two formulations must not be perceived as op- 
posed but conveying the same straightforward result. Meaning that any standard mi- 
croscope has a limited spatial resolution, which is about Ar,, = 4/2. 

By performing and experiment in which also the evanescent waves are detected, 
an infinite bandwidth of spatial frequencies can, in principle, be achieved. But this 
kind of measurement requires to probe the e.m. fields in close proximity to the sample 
surface, where evanescent waves are not yet totally suppressed. From this analysis it 
is clear that two spatial regimes, which show very different light behaviors, are rele- 
vant for any imaging purpose: 


Near-field r<K<a 
Far-field A<r 


where r is the distance between the specimen and the probe that collects the e.m. 
fields. The intermediate-field (r ~ A), also called “Fresnel-field”, is the vaguely de- 
fined as the region in between the near-field and the far-field. In this region, near-field 
behavior ceases to be important and far-field effects become dominating. 

The two extreme spatial regimes can be highlighted by studying the emission of 
a dipole, that is the smallest radiating e.m. unit source, which represents the first ap- 
proximation of any antenna emission [13]. In the far-field the magnetic and the elec- 
tric fields are, respectively: 


B=k?(n\p) PE? (1.7) 


E= Ban (1.8) 
where n = (r/r) is the unit vector of real space and the electric dipole moment 


is represented by p = f rp(r)d?r with p(r) being the charges density. In the far- 
field region we find propagating spherical waves where both field radiation intensities 
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decrease as the square of distance. As in every propagating light wave, both fields 
described by Eq. (1.7) and Eq. (1.8) are strongly connected since each of them is gen- 
erated by a variation of the other. Moreover, any absorption of the radiation does not 
feed back to the emitter. While in the near-field (kr « 1) they become: 


B = ik(naAp)r? (1.9) 
E= Bn(n-p)-p]r* (2.0) 


Differently from the propagating regime, in the near-field region the electric and 
magnetic fields exist independently of each other, and one type of field can dominate 
the other. In particular, the magnetic induction is a factor kr smaller than the electric 
field. Moreover, the near-field components depend on the detailed properties of 
source. For instance, Eq. (1.10) coincides with the electrostatic dipole, and both decay 
very rapidly as a function of the distance. Consequently, in the since a finite emitter 
can be regarded as an ensemble of Hertzian dipoles, near-field components carry the 
information needed to determine the spatial distribution of any finite size emitter. In 
addiction, in the near-field region absorption of radiation does strongly affect the load 
on the emitter. For instance, the magnetic induction in a transformer can be seen as a 
simple model of near-field e.m. interaction. 

In conclusion, it is straightforward that every time light is emitted or finds an 
obstacle and scatters there is a loss of information on propagating from the near-field 
region to the far-field. Optical techniques that allow near-field probing could improve 
the imaging hallmarks, reaching a sub-wavelength spatial resolution. However, the 
higher the spatial frequencies of an evanescent wave are, the faster the field decay 
along z will be. Therefore, unavoidable practical limitations always make the band- 
width finite. 


1.2 Modern microscopy 


Many fascinating and disruptive optical techniques have been proposed and de- 
veloped through the last decades revolutionizing the field of optical microscopy by 
gaining a sub-wavelength spatial resolution. At first, we briefly review the most suc- 
cessful methods that exploit specific optical schemes by collecting only far-field ra- 
diation and afterwords we focus the attention on near-field detection that is the basic 
research topic discussed in this thesis. 

Confocal microscopy was the first and simplest invention for increasing optical 
resolution and especially the signal-to-noise ratio [14]. It takes advantage of a spatial 
filtering performed by a pinhole placed in the conjugate plane of the specimen, so that 
light impinging the sample and light following the inverse path by passing through 
the pinhole have the same focus on the object plane as reported in Fig. (1.4 a). This 
geometry allows eliminating most of the out-of-focus light and prevents nearby scat- 
terers from contributing to the detected signal. The confocal microscope must detect 
only light coming from one point at a time, therefore a scanning system is required 
for reconstructing two or three-dimensional images. The drawback is that the lateral 
resolution does not overcome the diffraction limit although in the vertical direction a 
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narrow depth of field is achieved with respect to standard microscopy [15]. Neverthe- 
less, this non-invasive technique has gained popularity in biology, biophysics, semi- 
conductor inspection and materials science, particularly for investigating thick speci- 
mens or buried samples [16]. 

On the other hand, the invention of pulsed laser radiation propelled the field of 
nonlinear optics enabling, for instance, the invention of multiphotons microscopy [17, 
18]. With respect to confocal microscopy it offers the advantages of deeper penetra- 
tion and less photodamage but has the disadvantage of requiring a non-linear medium 
and may damage the sample due to the high power laser excitation. 


Recent advances in ultra-sensitive instrumentation have allowed for the detection, 
identification, and dynamic studies of single molecules in the condensed phase. This 
measurement capability provides a new set of tools to explore the frontiers in many 
scientific disciplines, such as chemistry, molecular biology, molecular medicine, and 
nanostructured materials. Optical detection of single molecules has been achieved by 
both frequency modulated absorption and laser-induced fluorescence. Since the latter 
exhibits lower background and higher signal-to-noise ratio, it is the more employed 
one [19, 20]. Highly fluorescent molecules have been synthesized that can be specif- 
ically attached to biological entities. This chemically specific labelling and the asso- 
ciated discrimination of different dyes based on their fluorescence emission allows 
scientists to visualize the interior of cells and study biochemical reactions in live en- 
vironments. 

Another noteworthy technique is the stimulated emission depletion (STED) mi- 
croscopy, which provides super resolution by selectively deactivating fluorophores 
thus enhancing the imaging in a localized area [21, 22]. STED works by depleting 
specific regions of the sample while leaving a center focal spot active to emit fluores- 
cence, leading to an improved resolution with respect to confocal microscopy as re- 
ported in Fig. (1.4 b). The focal area can be engineered by altering the properties of 
the pupil plane of the objective lens. In this content a sub diffraction method was also 
developed by using spatially structured illumination in a wide-field fluorescence mi- 
croscope employing a series of different excitation light patterns [23]. 

Photoactivated localization microscopy (PALM) and stochastic optical recon- 
struction microscopy (STORM) have also proved to be super resolution microscopy 
techniques allowing to image specific proteins down to few nanometers and trajecto- 
ries of single molecules [24-26]. They are both based on collecting under a fluores- 
cent microscope a large number of images as a function of time, each containing just 
a few isolated stochastically activated fluorophores, resulting in a promising approach 
to distinguish neighboring molecules by separating their fluorescent emission in time. 
All these results recently gave to Betzig, Hell and Moerner the honour to be jointly 
awarded with the 2014 Nobel Prize in Chemistry [27]. 


A completely different approach was recently brought to the realm of device en- 
gineering by Pendry in his seminal work [28]. He proposed the realization of a “su- 
perlens” based on a metamaterial slab. Metamaterials are artificial material with both 
negative dielectric permittivity and negative magnetic permeability. As if the refrac- 
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tive index was negative, these engineered materials allows optical properties unavail- 
able in naturally occurring substances, such as negative refraction, inverted Snell’s 
law and Cherenkov radiation and optical cloaking [29-34]. The superlens is, in prin- 
ciple, able to gain a perfect resolution by retrieving all the information carried by 
evanescent waves, since in the transmission process they emerge from the metamate- 
rial enhanced in amplitude after being resonantly coupled to surface plasmons, as 
shown in Fig. (1.4 c). Note that there is not any violation of energy conservation, since 
evanescent waves transport no energy, but nevertheless it is a surprising result con- 
firmed by experiments that in the visible spectral range show a resolution down to 
A/6 [35, 36]. 


Instead of photons, by using electrons techniques to magnify objects down to the 
nanometer range have been developed by exploiting the shorter De Broglie wave- 
length of electrons, such as scanning electron microscope (SEM) or transmission elec- 
tron microscope (TEM) [37, 38]. They involve an electron beam that generates a char- 
acteristic response from the sample surface (cathodoluminescence, back-scattered or 
secondary electrons) or is transmitted through an ultra-thin specimen, respectively. 

Significant advances in instrumentation are pursued by the wide class of scanning 
probe microscopes, where a small probe tip is brought into intimate contact with the 
sample surface without damaging it. The presence of the specimen manifests itself as 
a tunnel current flowing between the tip and the sample in scanning tunneling micro- 
scopes (STM) or as a deflection in the tip-height during the scanning in atomic force 
microscopes (AFM) [39, 40]. Both techniques push the spatial resolution beyond the 
Nanoscale and retrieve the image of single atoms. The former case provides a topo- 
graphic picture of the surface by keeping the tunnel current constant and monitoring 
displacements of the metal tip given by the voltages applied to piezoelectric drives. 
Unfortunately, electron microscopes have to be operated in vacuum, which limits their 
application in life sciences, and requires special sample preparation. The AFM is a 
higher sensitivity tool designed to investigate both conductors and insulators down to 
the atomic scale. The typical AFM setup reported in Fig. (1.5 a) shows that the canti- 
lever deflection is measured using a laser spot reflected from the top surface of the 
cantilever into an array of photodiodes. With this microscope images are obtained by 
measuring the force on a sharp tip (insulating or not) created by the proximity to the 
surface such as van der Waals forces, capillary forces, chemical bonding, electrostatic 
forces, magnetic forces and so on. This force is kept small and at a constant level with 
a piezoelectric based feedback mechanism to adjust the tip-to-sample distance, so that 
when the tip is moved sideways it follows the surface contours as is highlighted by 
the trace B in Fig. (1.5 b). Although these scanning techniques achieve a deep sub- 
wavelength resolution in mapping any sample topography, they cannot investigate or 
detect the light behavior in the same media and consequently retrieve the noteworthy 
spectral and dynamic properties of the sample. 

Therefore, all the described methods are complementary to those presented in this 
work, since, on one hand, PALM and STED are based on standard optical microscopy 
instrumentation and rely on fluorescent labeling of the sample, on the other hand, 
TEM and AFM provide a high-resolution topography but cannot probe the sample 
optical response. 
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1.3 Scanning near-field optical microscopy 


If AFM shear-force probe microscopy is combined with optical detection huge 
advantages can be achieved in the characterization of light emitting nanostructures 
and fluorescent molecules. The most striking example is near-field scanning optical 
microscopy (SNOM), in which the apex of a tapered or pulled optical fiber acts like a 
AFM tip, thus providing both optical and topographical information of the sampling 
surface. This integrated microscope works in close proximity with the sample surface, 
thus being strongly interacting with optical near-field and evanescent waves. Sub- 
wavelength imaging based on SNOM was experimentally tested for the first time in 
the visible spectral range more than 30 years ago by Pohl et al. employing aperture 
probes [41]. The principle of SNOM imaging operation relies on the point-by-point 
investigation by raster-scanning the local probe over the sample surface and recording 
at every position the corresponding optical signature. The local probe acts as an opti- 
cal antenna, converting localized e.m. fields into radiation, and vice versa. Therefore, 
SNOM detection gives ultra-sensitive optical measurement beyond the diffraction 
limit and provides light behavior at the nanoscale, together with the topographical 
capability of the AFM imaging. Although the morphological performance is not as 
good as in cantilever based AFM, SNOM proved to be the state of the art technique 
for investigating nanostructures such as quantum wells, quantum dots, nanoresonators 
or for single molecules since it combines the excellent spectroscopic and temporal 
selectivity of classical optical microscopy with the sub-wavelength spatial resolution 
[42-44]. 

The concept of near-field optical microscopy was originally proposed in by Synge 
1928 in a prophetic article, where he proposed an apparatus similar to the actual 
SNOM implementations [45]. A sub-wavelength aperture in an opaque plate placed 
in close proximity to the sample surface creates an illumination spot not limited by 
diffraction. Scanning the source along the surface and recording the transmitted, or 
reflected, light produces an image with subwavelength resolution. In fact, the resolu- 
tion of such an image should be limited only by the size of the aperture and not by 
diffraction, as was also pointed out years later by Bethe studying the diffraction by 
small holes [46]. Only after the invention of the STM an efficient positioning tech- 
nology was available at the nanoscale. Therefore an optical microscope similar to 
Synge’s proposal was realized, firstly in the microwave by Nicholls and Ash, then in 
the visible spectral range simultaneously by Pohl et al. and Lewis et al. due to the 
obviously different fabrications difficulties demonstrating a sub-wavelength imaging 
with a resolution of 1/60 and 4/20, respectively [41, 47, 48]. The key innovation was 
the fabrication of a sub-wavelength optical aperture at the apex of a sharply pointed 
transparent (quartz) probe tip whose sidewalls were coated by metal. Such sharp ap- 
erture replaces the Synge’s opaque screen for experimental purposes. In addition, a 
feedback loop was implemented based on proper piezoelectric quartz tuning forks 
(later replaced by piezoelectric bimorph bender element) maintaining the tip-sample 
distance constant (at few nanometers) while raster scanning the sample surface [49]. 
SNOM technique was systematically advanced and extended over the years to various 
configurations that provide a sharp confined photon flux in the near-field region, as 
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summarized in Fig. (1.6 a-d). Aperture probes are used in near-field illumination or 
detection, as reported in Fig. (1.6 a-b), respectively, or in the case of a sufficiently 
strong signal they can be even employed in illumination/collection geometry, Fig. 
(1.6 c). The last one is the best suited for improving the spatial resolution, since it 
provides an almost background free illumination. These configurations are suited of 
employing aperture or dielectric probe which, by exploiting the frustrated total inter- 
nal reflection that occurs at the fiber end, allow the conversion of evanescent waves 
with high spatial frequencies to e.m. fields propagating through the optical fiber and 
finally measured by a conventional far-field detector. 

In Fig. (1.6 d) SNOM setup is based on sharply pointed apertureless tip (a-SNOM) 
that locally enhances light scattering. Both illumination and collection are in the far- 
field, therefore a-SNOM requires a sophisticated method for filtering out the huge 
background. Moreover, in order to enhance light scattering apertureless probes are 
commonly fabricated on high refractive index material, such as silicon, that would 
drastically perturb strongly localized optical modes, which are the basic topic of the 
research presented in this thesis. 

In combination with time-resolved excitation and detection schemes, near-field 
microscopy has the potential to give direct access to processes on ultra-short time and 
length scales. Near-field microscopy has therefore all the advantages of conventional 
spectroscopy with a high spatial resolution solely defined by the probe size. 


1.3.1 Probes design and fabrication 


The fundamental parameter of a SNOM setup is the near-field probe, as it defines 
the spatial resolution, the amount of detected signal and the induced perturbation on 
a given sample. In this section, we review the basic near-field probe designs, fabrica- 
tion method and optical properties focusing on those that will be mentioned in the 
thesis concerning the experiments discussion. 

In e.m. theory an antenna is a device used to create e.m. waves with a well-defined 
radiation pattern, or to receive e.m. waves from a remote source in order to extract 
some encoded information, by measuring changes in their intensity, or by exploiting 
the transmitted power [50]. So antennas main ability relies in establishing efficient 
coupling between the near-field and the far-field by enhancing and concentrating the 
incoming fields and simultaneously allowing to control the directional emission of 
light sources placed in the near-field zone [51, 52]. Moreover, optical antenna can 
enhance light absorption, leading to advantageously photodetection and to photovol- 
taics applications [53]. Although antenna theory has been developed for the radio- 
frequency and the microwave range it holds great promise for inspiring new concepts 
in the optical frequency range. Therefore, the local SNOM probe can be viewed acting 
as an optical antenna, converting the investigated near-field into propagating radiation 
and vice-versa (because of reciprocity), with a net energy transfer efficiency that de- 
fines the antenna throughput. Standard aperture probes need a trade-off between spa- 
tial resolution versus throughput, since the transmission through a metal screen with 
a sub-wavelength hole of radius a scales as (a/A)*, with the same wavelength de- 
pendence of the Rayleigh scattering [46, 54, 55]. 
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Many different near-field probe designs have been conceived and have been 
tested in the last decades. They can be grouped in four categories: i) dielectric, i.e. 
transparent, probes that are also called uncoated tips, ii) metal coated aperture probes, 
iii) pointed probes (used in a-SNOM setup), iv) probes functionalized by placing a 
nano emitter or a plasmonic nanoantenna at the tip apex. 


Dielectric probes 


Dielectric tips are the fundamental class of near-field optical probes since they 
represent the key components for the fabrication of more complex probes. Tapering 
of dielectric optical fibers can be done by chemical etching, or by local heating and 
subsequent pulling. 

A well established method for fabricating smooth dielectric tapered probes is the 
tube etching method [56, 57]. We exploit this recipe to produce by ourselves the die- 
lectric probes employed the SNOM experiments that will be described in the thesis. 
Glass tip formation occurs inside a cylindrical cavity formed by the polymercoating 
of an optical glass-core fiber which is not stripped away prior to etching in aqueous 
48% hydrofluoric acid (HF). The two basic steps of the etching mechanism are re- 
ported in Fig. (1.7 a-b). In the first one, the acid largely corrodes the lateral side of the 
fiber core with respect to the central part since at the rim of the glass cylinder HF 
supply occurs out of a larger volume as compared to the central region. As soon as 
the preliminary taper is formed, step convection driven by concentrations gradients 
that are influenced by the etching process itself and the gravitational removal of the 
reaction products delivers HF to the upper cone region. After about 90 minutes of 
tube-etching the fiber is rinsed with ethanol and successively the jacket is removed 
without damaging the cone. The solution is covered by an organic overlayer (isooc- 
tane) to protect the fiber against acid vapor. Moreover, since the etching process takes 
places inside a hollow cylinder formed by the protective jacket, it is less sensitive to 
environmental influences such as vibrations or temperature drifts. The procedure 
gives tapered probes with large cone angles (about 30°) and small apex (less than 
100 nm size), resulting in a strong decrease of the loss region with respect to pulled 
probes [58], and low sidewall surface roughness, since the meniscus formed by the 
acid along the fiber does not move during the etching. Dielectric tips with reproduci- 
ble shapes and optical performance are formed in a high yield, as the one reported in 
the SEM image of Fig. (1.7 c). Since glass is almost transparent in the visible and 
near-infrared spectral range, these tapered dielectric probes show a high throughput 
(more than 50%), thus being suited for the high-contrast illumination-collection ge- 
ometry. 

However, the drawback is the low spatial resolution in the direct imaging. This 
feature can be simply explained if we consider how the fundamental mode of the op- 
tical fiber interacts with the tapered apex and it is focused in the near-field. The ta- 
pered, conical part of the probe can be represented as a series of glass disks with de- 
creasing diameters and infinitesimal thicknesses. With each step part of the radiation 
is reflected, and the transmitted mode becomes less confined, since the field extends 
more and more into the surrounding medium (air) and exhibits increasing losses. By 
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reciprocity, it is clear that dielectric probes collect not only evanescent waves in the 
near-field by frustrated total internal reflection at the tip apex, but also a considerable 
amount of the propagating e.m. far-field injected through the upper part of the cone. 
Therefore, by collecting also large K,, high throughput but poor spatial resolution are 
expected for uncoated probes and the use of such a probe is not always convenient. 


Coated aperture probes 


Tapered rather than pulled glass based tips, after metal coating in high vacuum 
evaporator by sputtering, can be modeled to result in small well-defined apertures at 
the apex, showing good homogeneity and negligible leakage of light along the cone 
[12]. Typically the metal is an aluminum or sometimes a gold layer of about 100 nm 
thickness, but the former is preferred for the smallest skin depth in the visible spec- 
trum. The aperture is formed by metal shadowing evaporation (that is rotating the 
probe along the cone axis during the coating process) and it is typically between 
100 nm and 200 nm diameter size for balancing high resolution and good throughput, 
see Fig. (1.8 a) [59, 60]. 

A more sophisticated fabrication technique allows for obtaining a well-defined 
flat circular aperture with controllable diameter size down to 20 nm, by milling the 
coating with a focused ion beam (FIB) [61]. FIB is the state of the art method to build 
high quality nanoscale structures; however, it is a rather expensive and inherently slow 
technique. Significantly less expensive and reproducible procedures have been put 
forward for the reliable fabrication of aperture probes as controlled electrolysis, aper- 
ture punching, selective resin coating or by integrating a waveguide into the AFM 
cantilever making the most of established silicon engineering techniques [62—65]. 

The lateral metal coating avoids the leakages in the outer space thus confining the 
e.m. fields at the aperture. In fact, by approaching the tapered tip apex light is pro- 
gressively confined, but exceeding a transverse diameter size smaller than the wave- 
length of the modes propagating in the fiber, all the sustained modes are back reflected 
or absorbed and dissipated by the metal and cannot propagate any further. For this 
reason, the sub-wavelength aperture probes introduce a diameter-size dependent cut- 
off causing only evanescent waves to escape, as highlighted in Fig. (1.8 b) [12, 60]. 
Of course, this characteristic gives an experimental extremely small throughput with 
respect to dielectric probes (with a ratio of about 107°) but allows to obtain high 
spatial resolution since only the near-field is excited/detected [66]. 


Apertureless probes 


Also by using narrow pointed apertureless probes, made of metal or high refrac- 
tive index dielectric such as silicon, it is possible to retrieve optical information during 
the scan by exploiting the strong field enhancement in the proximity of the antenna 
that for gold tips is predicted to be very large, up to 50 [67]. Of course, the only avail- 
able configuration is collecting in the far-field the scattered light tip-sample system 
that carries information about the near-field interaction, as shown in Fig. (1.6 d). The 
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illumination can be provided by standard microscopes as well as by evanescent waves 
from total internal reflection [68]. The optical signal is measured by inducing a mod- 
ulation in the tip oscillations and performing homodyne or heterodyne interferometric 
detection using lock-in techniques in order to demodulate the beat signal at higher 
harmonics of the tip vibration; therefore the far field part can be suppressed effectively 
to discriminate the small scattered signal from the tip apex against the diffraction lim- 
ited background, leading to a few nanometers resolution in the infrared spectral range 
[69-71]. 


Functionalized probes 


In order to manipulate and confine optical radiation at subwavelength scales a key 
improvement has been developed by integrating single light sources or nanostructured 
devices on the apex of standard probes. In this way a further nanoantenna is created 
on top of the near-field probe and a tremendous advancement in terms of spatial res- 
olution, throughput and sensitivity is achieved. 

On one hand, the modification of the spontaneous emission rate of an emitter 
placed in a nano-optical resonator device (Purcell effect [72]) can be monitored by 
grafting the light source at the end of a scanning probe, see Fig. (1.9 a), and then 
performing fluorescence lifetime measurements as a function of position [73, 74]. 
However, the disadvantage is that single nano-emitters may suffer of bleaching and/or 
blinking [75, 76]. 

On the other hand, thanks to high reliable nano-fabrication procedures which in 
most of the cases exploit FIB or electron-beam lithography, metal-based nano-anten- 
nas milled at the apex of standard probes are fabricated in order to achieve outstanding 
spatial resolution and bright near-field probes, thus overcoming the throughput con- 
strain imposed by the aperture dimension, by squeezing the near-field by means of 
localized plasmons coupling [77]. 

It is a well known principle of microwave engineering that the most efficient way 
to couple far-field radiation to a device of dimension much smaller than the light 
wavelength is to connect the device to an antenna [50]. If we follow the parallelism 
with antenna theory, in the optical regime we need to introduce metal materials losses, 
however such losses may be compensated by collective electron resonances, i.e. sur- 
face plasmons [78]. In fact, the recent enormous advance in nano-imaging techniques 
coincided with the introduction of concepts from plasmonics, such as surface plasmon 
polaritons (SPPs) confinement along metal-dielectric interfaces, into the design of 
nearfield probes [79, 80]. SPPs are waves that propagate along the surface of a con- 
ductor. By altering the structure of a metal surface the SPP interaction with light can 
be tailored on demand, offering the potential for developing new types of photonic 
devices, making optical antennas particularly efficient at selected frequencies and 
broadening the perspectives of near-field optics. With respect to metallic coated 
probes, plasmonic functionalized tips exhibit either a better spatial resolution (at con- 
stant throughput) or larger throughput (at constant resolution). Several high-perfor- 
mance near-field probe geometries, in which plasmonic nanoantennas are included, 
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have been engineered and optimized both experimentally and theoretically with re- 
spect to size, shape and material in order to efficiently transform light from the far- 
field to the near-field and the other way round, and for increasing the radiation effi- 
ciency of nano-optical emitters [81-83]. 

Although the design of optical antennas is likely to be inspired by radio-frequency 
or microwave regime working geometries, it can be expected that the presence of new 
physical phenomena will demand the exploration of new structures and material com- 
positions. The first developed plasmonic geometry was the renowned bow-tie nano- 
antenna that is constituted of two metal triangles facing each, as shown in Fig. (1.9 b), 
and optimizes the impedance matching, thus delivering high throughput and local 
field enhancement [84-86]. The practical realization on top of a scanning probe is 
performed by FIB, creating a bow-tie nano-aperture at the apex of a coated probe [87]. 
Other realizations are, for instance, mini-tip on aperture, single metal nanoparticle, 
see Fig. (1.9 c), and diabolo (complementary to the bow-tie) antenna structures have 
been grown on the end faces of aperture-type near-field probes [88—92]. However, all 
the reported functionalized probes rely on resonant structures with a limited spectral 
bandwidth that impose the spectral matching between the field enhancement effect 
and wavelength operation of the investigated sample. 

In Chapter (3) we will discuss SNOM experiments performed with a recently de- 
veloped next generation plasmonic nanoantenna, called “Campanile tip”. This near- 
field probe milled by FIB on top of a gold coated tip shows a broadband operation 
and allows us for mapping with unprecedented spatial resolution the influence on car- 
rier recombination in local trap states within individual solid state emitters (quantum 
wires) as well as the simultaneous electric and magnetic field distributions confined 
in nanoresonators. 
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Figure 1.1: a) Schematics concerning the image formation on the plane Z;of a point- 
like object (S) illuminated by a monochromatic plane wave with wavevector K per- 
formed by the optical system L (lens or objective). @ is the half-angle of the maxi- 
mum light cone accepted by L and Èp is the back focal plane. Red curves describe 
the behavior of radiation field wave front. The blue curve represents a one-dimen- 
sional profile cut of the Airy disk in 27. b) Typical three-dimensional Airy disk in- 
tensity pattern representing the image of the point-like source in Z,. 


Figure 1.2: a) Total internal reflection at the interface between two media withn, < 
nı if 0, > arctan (nz /n,). The displacement d represents the Goos-Hanchen 
shift. The evanescent wave in nz medium is schematically reported exhibiting a 
Poyniting vector S parallel to the interface. b) Frustrated total internal reflection oc- 
curs in a successive nj medium placed in close proximity to n, if the evanescent 
wave is not completely decayed. 
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Figure 1.3: In the angular spectrum representation the e.m. fields are evaluated in a 
plane perpendicular to an arbitrarily chosen z axis. In this case we evaluate the total 
field Aine + Ascatt- 
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Figure 1.4: a) Principle of confocal microscopy: the pinhole in front of detector fil- 
ters out the light that does not come from the focal spot. b) Comparison between 
normalized point spread function (PSF) of STED fluorescence microscope, confocal 
microscope and conventional scanning microscope reported by a, p and y curves, 
respectively, as evaluated in [21]. c) Pendry’s superlens principle of operation. Neg- 
ative refraction occurs inside the negative refractive index media (—n,) leading to a 
double focusing inside and outside the slab. The red curve describes the evanescent 
wave decay in n, and its amplification inside the superlens, so that in the image 
plane the initial evanescent wave amplitude is retrieved. 
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Figure 1.5: a) Principle of AFM microscopy: the tip attached to the cantilever scan 
the sample surface following the blue curve B reported in b), while the laser deflec- 
tions are measured by a four-quadrant photodiode. The electronic feedback control 
assures a constant tip-sample distance. b) Zoom on the tip-sample gap that high- 
lights the path scanned by the tip. Only the tip front atoms inyeract with the speci- 
men [39]. 
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Figure 1.6: Schematics of the four more widespread SNOM configurations. a) In- 
coming light illuminates the sample in the near-field by means of an aperture probe 
and the scattered light is detected in far-field by standard microscopes. b) Comple- 
mentary regime with respect of a): far-field illumination and near-field detection. c) 
Illumination/collection geometry in which both the near-field illumination and col- 
lection occur through the tip. d) Apertureless probe enhances the scattering in the 
near-field, whose information are analyzed by far-field detection. 
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Figure 1.7: Dielectric tapered fiber etching process in aqueous hydrofluoric acid 
[56]. a) First etching step: diffusive motion of the HF solution. b) Second etching 
step: convective motion of the HF solution. c) SEM image of an etched dielectric ta- 
pered fiber. 


Figure 1.8: a) SEM image of an aluminum coated probe with a nominal aperture of 
200 nm crated by shadowing evaporation, which is commercially available by the 
company LovaLite. b) Aperture probe lateral view schematics. It is highlighted the 
cutoff of the propagating modes when they approach the tip apex with an aperture 
size smaller than 4 that allows only evanescent waves to escape. 
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a) b) c) 


Figure 1.9: Schematics of functionalized near-field probes. a) Single molecule light 
source grafted at the end of the tip. b) Bow-tie aperture at the apex of a metal coated 
probe. c) Gold nanoparticle at the tip apex acts like a plasmonic nanoantenna. 
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Chapter 2 


Nanophotonics localized fields 


The nano-photonics research field has been quickly growing since the introduc- 
tion of artificial photonic crystal media engineered at the nanoscale. In particular, 
photonic crystal cavities fabricated on semiconductor slab are the state of the art de- 
vices to strongly localize light fields in volumes comparable to (A/n)? acting as effi- 
cient nano-resonators. They provide an efficient mechanism for controlling and inhib- 
iting spontaneous light emission, for reducing almost to zero the light group velocity 
and they can also be easily included in all-optical integrated circuits. Here, we de- 
scribe the basic properties of photonic crystal nanocavities. In particular, how, on one 
hand, SNOM techniques can experimentally probe the e.m. optical fields localized at 
the nanoscale and, on the other hand, how finite elements numerical calculations can 
simulate the behavior of light in such nano-resonators. In addiction, we address also 
the description of optical nano-resonators based on disordered arrangements of light 
scatterers. In these systems the interference between all the possible photon random- 
walks, due to multiple scattering, gives rise to random spatially and spectrally placed 
finite quality factor nanocavities. 


2.1 Photonic crystal nanoresonators 


2.1.1 Photonic crystals 


Nowadays, the most widespread and versatile devices conceived to mold the flow 
of light at the nanoscale are represented by artificial photonic crystals, which were 
introduced to the scientific community in the 1987 [1, 2]; although a century before 
Lord Rayleigh already studied a simplified analogous structure [3]. Photonic crystals 
are low-loss dielectrics systems in which a periodic modulation of the refractive in- 
dex, with a period comparable to A, is present. Dielectric modulation can be created 
along one, two or even three spatial dimensions thus defining 1D, 2D or 3D photonic 
crystals, respectively. As the ions in solid state materials occupy defined lattice posi- 
tions inducing an effective periodic potential to electron wavefunctions, the refractive 
index arrangement defines a lattice that gives the name to this engineered periodic 
crystal, in which the e.m. fields experience an effective periodic potential [4]. 

Due to the multiple scattering of light at the dielectric interfaces, interference phe- 
nomena deeply affect the the e.m. wave propagation inside photonic crystals. By en- 
gineering different periodic modulation geometry and also by altering the amplitude 
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of the refractive index contrast it is possible, for instance, to prevent light from prop- 
agating inside the photonic crystal in certain directions and with specified wavelength. 
Therefore, being fully back-reflected when it shine on the photonic crystal external 
surface [4]. Even thought, the behavior of electrons in semiconductors is different 
with respect the light propagation in photonic crystals the analogy can include Bloch 
formalism since photonic band-gaps resemble the electrons energy band-gap in sem- 
iconductors [5]. If the lattice potential is strong enough the photonic energy gap can 
extend over all possible directions, resulting in a complete photonic band-gap, as it 
happens in semiconductors for electron energy between the valence and the conduc- 
tion band [6-8]. 

Since fabricating photonic crystals in the optical spectral range requires a dielec- 
tric modulation of the order of (A/n), it is mandatory a huge technical effort, espe- 
cially in three-dimensional geometries. Due to the enormous progress in lithographic 
techniques established over the past decades for engineering semiconductor systems, 
fabrication resolutions down to few nanometers size are today attainable, allowing the 
manufacturing of high quality and high refractive index contrast photonic crystals op- 
erating at optical frequencies [9]. 

In this context, it is not surprising that nature has beaten us to the draw, since 3D 
photonic crystals can be largely found. In fact, they are responsible for the structural 
color that, as opposed to colors resulting from absorption by chemical pigments, ena- 
bles the bright and iridescent long-lasting colors owned, for instance, by a large vari- 
ety of butterflies, peacocks, jellyfishes and other animals (profitably due to millions 
of years lasting evolution) as well as by some opals and berry fruits [10-13]. 

The description of light behavior in photonic crystals must involve the solution of 
Maxwell’s equations in a periodic dielectric medium. The usual simplified hypothesis 
are the absence of e.m. sources, homogeneous and neutral linear dielectric media (with 
a magnetic permeability (4, = 1) so that the dielectric constant (e = n?) is a scalar, 
which depends only on the spatial coordinate. After factorizing the temporal evolu- 
tion, it follows that the magnetic field has to satisfy the Master equation (in c.g.s 
units) [4]: 


1 wy? 
{VA 5 A} H(r) = (È) H(r) (2.1) 
while the corresponding electric field can be evaluated as E(r) = (iw/c) A H(r). 
Solving Eq. (2.1) is equivalent to an eigenvalues problem for each field component, 
like the stationary Schrödinger equation. Therefore, if e(r) is perfectly periodic the 
solutions of Eq. (2.1) are characterized by a wavevector k and a band index n: 


H(r,t) = Hkn (r)exp(ik -r — iwynt) (2.2) 


where Hyn(r) has the same periodicity of the dielectric modulation e(r). 
Eq. (2.2) does not differ for any translation in the k-space equal to an integer multiple 
of the reciprocal lattice wavevector. 

A reliable fully-vectorial computational technique that allows to solve the Master 
equation in complex periodic structures is the plane wave expansion method (PWE). 
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By exploiting this method is can be easily calculated band structure (dispersion rela- 
tion) of any specific photonic crystal geometries [14, 15]. The algorithm expands in 
terms of Fourier series along the reciprocal lattice vector the periodic Bloch functions 
that represent fields and dielectric permittivity. Of course, this Fourier expansion is 
truncated to a fixed number of terms, limiting the accuracy of the calculation. A no- 
table example is the band diagram of a multilayer periodic in one dimension (1D pho- 
tonic crystal) calculated by PWE, as reported in Fig. (2.1). The lower energy band is 
denominated dielectric-band since it represents modes whose electric field maxima 
mainly occur in the higher refractive index layer. On the other hand, the higher energy 
band is called air-band, since it describes modes showing electric field maxima in the 
lower refractive index layer (that is usually air). PWE predicts also the positions and 
the widths of photonic band-gaps. 

The more convenient parameter to understand and quantify how photonic crystals 
mold the light behavior is the local optical density of e.m. states (LDOS) or g(r). This 
quantity, in analogy with solid-state physics, accounts for the number of available 
optical modes per unit of frequency and volume. Usually it is derived from the Green 
function of the e.m. system by taking the imaginary part of the trace of the Green 
dyadic summed over all the possible states [16]. Furthermore, it is a common strategy 
to discern between the electric and magnetic contribution and, especially when con- 
cerning light-matter interaction, to define partial LDOS for particular dipole direc- 
tions. For instance, the electric LDOS average over all solid angle results: 


olr, w) = Xn J, py dk 8(0 — xn)|Exnl (2.3) 


Eq. (2.3) highlights the fact that electric (magnetic) LDOS is related to the inten- 
sity of the electric (magnetic) field associated to all the modes that exist at fre- 
quency w. Furthermore, it can be proved that LDOS is inversely proportional to the 
group velocity of light [5, 17]. In order to understand the differences in light propa- 
gation induced by the presence of a periodic dielectric modulation, we report the 
LDOS in vacuum: 


w2 


Qvac(T, @) = (2.4) 


that is homogeneous and isotropic. Otherwise, inside a photonic crystal the LDOS 
experiences a strong inhomogeneity. Remarkably, in correspondence of the photonic 
band-gap LDOS vanishes. Since the spontaneous emission rate of a light emitter, ac- 
cording to Fermi’s golden rule, is proportional to the e.m. LDOS at the emitter’s po- 
sition, when the light source is located inside the photonic crystal and its energy emis- 
sion lies in the photonic band-gap, the spontaneous emission can be fully suppressed. 
This was also the first Yablonovitch’s idea, which led to the introduction of photonic 
crystals [1]. The LDOS is a crucial parameter in photonics because it is the only quan- 
titative way to describe the continuous part of the spectrum of any system inde- 
pendently of the excitation mode. This means that the LDOS provides spectroscopic 
information that are intrinsically independent on the illumination mode. 
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Photonic crystals based devices have been a strong impact in a large variety of 
optical applications. Concerning optoelectronic, a noteworthy result is that in 1D pho- 
tonic crystal Bragg reflection will occur and light with wavelength 4 incident along 
the periodicity direction will be reflected back with nearly 100% efficiency if the 
layered media is a quarter-wave stack, that is the period of the medium formed by nj 
and n, is equal to (A/2n, + 4/2n,) [18]. So they are used in modern devices such 
as the highly efficient and selective dielectric mirrors, dielectric Fabry-Pérot filters, 
in distributed feedback lasers, low-loss waveguides and in photovoltaic applications 
[4, 19-23]. It is possible to engineer photonic crystals in order to control photons dis- 
persion relation for achieving almost zero group velocity (“slow light” effect) to max- 
imize radiation/matter interaction, leading also to novel thin-film solar cells [24, 25]; 
or for creating an environment where temporally narrow pulses can travel great dis- 
tances without being extended or even being compressed, thus enhancing non-linear 
effects [26, 27]; or even for realizing optical metamaterials that exhibit a negative 
refractive index [28]. 


2D Photonic crystals on slab 


Although fabrication breakthroughs in layer-by-layer deposition and in self as- 
sembly structures drastically improved the performances of 3D photonic crystals, 
there are still large difficulties in obtaining efficient, reproducible, cost-effective de- 
vices exhibiting a three-dimensional modulation at the nanoscale [29]. In order to 
overcome these limitations the well-established manufacturing techniques that em- 
ploy semiconductor nanostructures has been emerged as a powerful tool to manipulate 
light in the three dimensions without the need to create 3D photonic crystals [30]. 
Fabricating 2D photonic crystals on a semiconductor slab suspended in air, with a 
thickness comparable to the in plane lattice constant, represents a smart strategy to 
solve the issue. In the out-of-plane direction, the slab structure provides an index guid- 
ing by total internal reflection, as in standard optical fibers. The amount of light that 
escapes from the 2D photonic crystal plane radiates outside the “light cone” (a/A = 
ak ;/2mn). In Fig. (2.2 a) is reported the design of a 2D photonic crystal of air holes 
arranged on a hexagonal lattice with a filling fraction (f, defined as the ratio between 
the area occupied by voids and the area occupied by dielectric medium in the unit cell) 
equal to 35%. The relative PWE calculated energy band diagram is reported in 
Fig. (2.2 b) in the first Brillouin zone, where the high symmetry points of reciprocal 
space are I =(0, 0), M=(7/a, v37/3a) and K = (47/3a; 0). It shows a photonic 
band-gap beneath the light line for the TE-like modes only, which exhibit a magnetic 
field mostly orthogonal to the slab plane. Actually, in a real slab only TE-like modes 
can be defined. While pure TE modes, with the magnetic field exactly orthogonal to 
the slab plane, exist only for pure 2D systems. In fact, the slab symmetry causes the 
modes to be purely transverse only in the mirror plane of the vertical direction, 1.e. in 
the middle of the slab. In conclusion, the realization of 2D photonic crystal hybrid 
slab structures allows to discard the quest for a periodicity in three dimensions and for 
a complete band-gap, since it conveys all the advantages of photonic crystal into 
standard light-guiding devices. 
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A further advantage of 2D photonic crystals is that nanoscale fabrication can eas- 
ily induce in them intentional and reproducible defects in the periodicity that can lo- 
cally sustain strong light confinement. In the thesis we will deal only with 2D photonic 
crystal on slab based devices, in particular we will focus the attention on hexagonal 
lattice of air holes just like the one reported in Fig. (2.2). 


2.1.2 Photonic crystal nanocavity 


The e.m. fields, solutions of Maxwell’s equations Eq. (2.1), with energies inside 
the bang-gap do not propagate in the photonic crystal. Nevertheless they are present 
as oscillating evanescent waves (whose intensity decay exponentially with superim- 
posed an oscillatory behavior) [31]. The scheme of permitted and prohibited modes 
can be altered by introducing a narrow bandwidth and spatially localized mode inside 
the photonic band-gap if the discrete translational symmetry, typical of photonic crys- 
tals, is broken by perturbing even a single site lattice, that corresponds to a defect in 
the periodic lattice. In fact, the LDOS becomes strongly modified in the close prox- 
imity of such defects and can be enhanced by many orders of magnitude. Conse- 
quently, the local symmetry breaking can trap light of a given J and behaves just like 
an optical cavity, being surrounded by perfectly reflecting walls due to the photonic 
band-gap. The analysis of the LDOS in a photonic crystal modified by the presence 
of defects tighten the analogy with semiconductor physics, where the doping is a com- 
mon technique able to introduce discrete energy modes inside the forbidden band-gap. 
In photonics, if the defect involves removal of dielectric then the localized mode in 
the cavity evolves from the dielectric-band (lower-energy) and can be made sweeping 
across the band-gap by adjusting the amount of dielectric removed. Similarly, if the 
defect involves an addiction of extra dielectric material then the cavity mode drops 
down from the air-band (higher energy) [4]. Therefore, acting on the defect geometry 
and on the refractive index it is possible to tune the localized mode wavelength to a 
target level inside the band-gap. There have been developed many types of photonic 
nanoresonators, most of them rely on planar photonic crystal on slab and a complete 
review can be found in [32]. We will deal with photonic crystal nanocavities (PCN) 
formed by the removal of one or more air holes in a 2D photonic crystal slab analo- 
gous of the one reported in Fig. (2.2). The typical volume occupied by such cavity 
modes can be very small, of the order of V = (A/n)?, and the quality factor Q = 4/44 
(where AA is the resonant mode linewidth) can become sufficiently large, up to 
10° [33]. The most widespread type of PCN consists on the absence of three adjacent 
holes along the unit cell direction, denominated L3 cavity [34]. 

In order to achieve ultra-high Q, it has been proved that the PCN realization 
should confine light gently to be effective, meaning the defect must not be an abrupt 
hole removal or hole filling. For instance, by inducing a displacement on the air holes 
that surround the L3 cavity the out-of-plane leakages are minimized since the field 
inside the light cone is reduced [33]. Also by carefully shifting more holes by super- 
imposing one or more additional periodic lattices to the photonic crystals, it is possible 
to control the angular out-of-plane emission pattern by making the band-diagram fold 
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in the Brillouin zone and enhancing the Q [35]. Recently, exploiting a genetic optimi- 
zation of the positions of few nearby cavity holes, a quality factor of the conventional 
L3 nanocavity that exceeds 10° was experimentally measured [36]. 

The realization of ultra-high Q and small mode volume PCNs is of the utmost 
relevance for a large research field that studies the interaction between a quantum 
emitter and a single radiation-field mode, called solid-state cavity quantum electrody- 
namics (CQED) [37, 38]. Experiments on atom-like emitters embedded in cavities can 
be described by elementary models, but reveal intriguing subtleties about the interplay 
of coherent dynamics with external couplings. 

In the noteworthy case of interaction between localized cavity fields and atom- 
like emitters, schematized as dipole light sources, we can discern two different re- 
gimes: the weak-coupling, characterized by an interaction that is lower than both the 
emitter and cavity mode decay rates, and the strong-coupling where the light/matter 
interaction is strong enough that the two initial states cannot be seen as separate enti- 
ties [39, 40]. In the weak-coupling regime Purcell predicted that the spontaneous 
emission rate could be tailored by modifying the density of available photon modes 
(LDOS) of the bare cavity [41]. For an ideal electric dipole source, spectrally resonant 
and with the same polarization of the cavity mode, placed at the maximum of the 
localized electric field, the nanocavity figure of merit, which gives an estimation of 
the electric spontaneous emission enhancement with respect to the emission rate when 
the dipole is in vacuum, is defined as the Purcell factor: 


r= 38) 2s 


where the electric mode volume is defined as: 


_ SJeMlE@MIzar 
az max[e(r)|Eo(r)|2] (26) 


with Eo(r) representing the electric field of the resonant mode of the bare 
nanocavity. Therefore, the Purcell factor of Eq. (2.5) can be strongly enhanced in pho- 
tonic crystal nanocavities by exploiting the deterministic fabrication techniques that 
allow for accurately positioning light emitters like quantum dots (QDs) in precise 
PCN positions [42-44]. The contribution of photonic crystal nanocavities triggers also 
the evolution of a variety of applications ranging from efficient light-emitting-diodes 
to ultra-fast and low-threshold lasers, including illumination, optical communication, 
solar energy and even quantum-information systems [45, 46]. 

On the other hand, if the emitter and the localized field are in the strong-coupling 
regime, coherent entangled states of quantum superposition exhibiting characteristic 
Rabi temporal oscillations and Rabi spectral splitting, even in the zero-detuning con- 
dition (which implies a perfect spectral matching between the bare resonant mode and 
the isolated emitter), has been predicted and experimentally observed [47, 48]. This 
feature extends the application range of CQED to fundamental physic concerning 
basic light/matter interaction. As an example, in the strong coupling regime nonlinear 
optics occurs even with only two photons, giving rise to the photon-blockade effect 
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[49]. In fact, the transmission or reflection from the cavity can be controlled by ab- 
sorption of just a single photon, since it modifies the energy states distribution of the 
system and prevent a subsequent identical photon to interact with the system. 

Self-assembled QDs are preferred candidates for the observation of these intri- 
guing light/matter interactions since they can be easily embedded inside 2D photonic 
crystal cavity on slab by exploiting the well-established molecular beam epitaxy sem- 
iconductor growing technique [50, 51]. Moreover, they are space-fixed nanostructure 
with a stable atom-like emission spectrum that reaches few ueV linewidth at low tem- 
peratures [52]. Single QDs have been reported as reliable non-classical light sources, 
able to generate defined single photon states. It is therefore very useful for quantum 
information and quantum cryptography applications to include these kinds of emitters 
into photonic nanocavities, which allow an efficient information transfer between a 
solid-state qubit and a flying qubit, as the emitted photon coupled to the cavity 
mode [53]. 

It is evident that the development of integrated nano-optical resonators is a fun- 
damental step in the route towards ultrahigh density all-optical integrated circuits as 
well as to all-integrated quantum photonic networks, since they are suited to stand for 
the building blocks of these systems such as optical switches and add-and-drop filters 
[54-57]. For instance, a single photon transistor in which one stored gate photon con- 
trols the resonator transmission of subsequently applied source photons has been re- 
cently demonstrated [58]. Photonic crystal fibers involves a linear defect line and are 
an established 2D photonic crystals device for fiber-optic communications, fiber la- 
sers, nonlinear devices and high-power transmission [59]. Moreover, exploiting 
strong light-matter interaction between photons and excitons in solid state quantum 
wells confined in PCN, polariton quasi-particles have been observed [60]. Also 
AND/OR gate operation has been shown, thus validating the connectivity of multiple 
transistors in the microcavity plane and paving the way to the implementation of po- 
lariton integrated circuits. Recent activity in this area has pioneered powerful ap- 
proaches to the study of quantum coherence and has fueled the growth of quantum 
information science. 

In conclusion, and more importantly, the characterization of e.m. fields localized 
in photonic crystal nanocavities with a sub-wavelength spatial resolution is essential 
for proving the intimate proprieties of CQED effects and for efficiently coupling light 
emitters in photonic devices. 


2.1.3 Finite-difference time-domain calculations 


Accurate analytic evaluation of the e.m. LDOS in complex photonic structures 
such as PCNs it is a difficult task. Plane wave expansion calculation method, pre- 
sented in Chapter (1), is not suited for estimating the spectral properties and the field 
distributions of the localized modes. A useful common numerical technique devel- 
oped to address this quest is the finite-difference time-domain (FDTD) method [61, 
62]. It is used for modeling full-vector e.m. waves evolution in arbitrary structures by 
finding approximate solutions to the associated time-dependent Maxwell’s equations 
over time in any specific positions. FDTD algorithms involve both the discretization 
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of the wave equations, by using second-order central-difference approximations to the 
space and time partial derivatives of the fields, and the mashing of the dielectric en- 
vironment by a dense discrete grid. Each cell of the grid has its own electric permit- 
tivity together with the electric and magnetic field defined both in the center and at 
the edges of the unit-cell faces, as reported in Fig. (2.3 a). The resulting finite-differ- 
ence equations are solved in a leapfrog manner by applying the usual tangential field 
continuity condition at each cell interface. In particular, the value of a field vector 
component at any lattice point depends only on its previous value and on the previous 
value of the components of the other field vector at adjacent points. Therefore, the 
updated electric field in time is dependent on the stored electric field and on the spatial 
variation of the magnetic field, as highlighted by the simplified scheme reported in 
Fig. (2.3 a). The process is iterated for every field component until the desired transi- 
ent or steady-state e.m. field is fully evolved. When the stability criteria are fulfilled, 
the FDTD algorithm provides a reliable computing scheme [61]. 

FDTD algorithms finally involve a Fourier-transform of the calculated time-var- 
ying fields to obtain the system spectral response. The electric and magnetic field 
spatial distributions at a given wavelength are therefore easily achieved, thus giving a 
direct estimation of both the electric and magnetic LDOS of the investigated structure. 

Typically, in any FDTD calculation extraordinary care must be taken in designing 
the dielectric environment, the light sources properties (dipoles, plane waves or 
Gaussian beams with defined polarization, wavelength, bandwidth, intensity, position 
that have to be suited for exciting a given PCN eigenmode) as well as in defining the 
actual “detectors”, which represent the positions whose fields are stored permanently 
in the calculator memory as a function of time. Since it is a time-domain method, 
FDTD solutions can simultaneously cover a wide frequency range. Moreover, the di- 
rect benefit in representing the field vectors in every cell, is achieving LDOS distri- 
butions with a high spatial resolution. It is, in fact, defined by the dimensions of the 
single cell, whose lower limit is given by the maximum memory that can be allocated 
in the calculator. On the other hand, the spectral resolution is bounded to the number 
of time-steps for which the fields evolve, and then it is limited both by the calculator 
memory and by the overall calculation time. 

Nevertheless, due to the time domain discretization, dispersive materials such as 
metals are difficult to handle with FDTD. A calculation procedure that overcomes this 
limitation performing approximate solution of partial differential equations is the fi- 
nite element method (FEM) [63]. The approach is based on eliminating the time de- 
rivatives, thus investigating only steady state problems. It is a fast method to obtain 
the stationary spatially resolved field distributions at a fixed wavelength for modes at 
least partially localized in any resonators, since the result is equivalent to a very long 
FDTD calculation. 

In conclusion, finite-element calculation methods represent the main theoretical 
predictions with whom compare the experimental data. Nowadays, many commercial 
finite-element are available. In this thesis the FDTD theoretical predictions concern- 
ing the LDOS of modes localized in PCNs will be calculated by means of Crystal 
Wave (by Photon Design), while the FEM calculations to estimate the local field en- 
hancement in plasmonic functionalized near-field probes will be performed by Com- 
sol. 
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As a remarkable example, in Fig. (2.4) we report a FDTD calculation performed 
by Crystal Wave on a modified L3 PCN where two lateral holes are gently shifted and 
reduced. In Fig. (2.4 b) is reported the top-view design of the PCN on slab as it is 
stored in the FDTD solver calculator with a cubic grid of 22 nm size. The detector of 
the near-field is placed 22 nm on top of the slab. The photonic crystal is designed to 
have a lattice constant of 308 nm, f = 35% and a slab refractive index of n = 

3.484 (that is the same of GaAs for A = 1.3 um). The single time step is fixed at 
0.4756 fs and the overall time sampling in this case covers 131072 time steps. In the 
calculation design we embed e.m. dipole sources in the slab with different orientations 
and a spectral range that covers all the photonic band-gap, reported in Fig. (2.4 a). 
Following this approach we excite all the available photonic modes of the system, 
hence including the localized PCN modes but also the “leaky modes” (non-stationary 
states inside the photonic band-gap that are non completely inhibited contributing to 
the PCN leakages). Therefore, the spectrum averaged over the whole detector size of 
Fig. (2.4 b), represents the spatially integrated LDOS as a function of the wavelength, 
as reported in normalized arbitrary units in Fig. (2.4 a). The intensity is almost zero 
in a large spectral range, which corresponds to the photonic band-gap where the 
LDOS vanishes, while the narrow peaks indicate the allowed modes inside the band- 
gap. The contribution of leaky modes is negligible in this case and the electric LDOS 
is dominated by few very localized modes. Figures (2.4 c) reproduces the electric field 
intensity distribution calculated by FDTD at the peak spectral position that exhibits 
the higher intensity mode of the PCN (A = 1285 nm). Such spatially resolved calcu- 
lation can be directly identified with the electric LDOS at a given wavelength. From 
this map it is evident that the related experimental challenge concerns the probing of 
the nanoscale modulation of the electric LDOS localized in PCNs. 


2.2 Near-field probing 


In this section a state-of-the-art review of efficient and commonly employed 
methods to detect the light fields localized in photonic crystal nanocavities is pro- 
vided. We also describe some typical near-field measurements performed with the 
SNOM setup exploited in the following Chapters, which describe experimental re- 
sults. 


Although several developments have improved the SNOM sensitivity, versatility 
and resolution, a complete understanding of tip-light-sample interactions, and specif- 
ically which components of the electromagnetic field are collected by SNOM probes, 
is still lacking. In fact, this issue depends on many subtle aspects such as the exact tip 
geometry, the orientation of electric and magnetic dipoles with which the probes can 
be approximated, the polarization of the investigated light, and more in general on the 
performed experimental configuration. By exploiting the reciprocity theorem, it is 
proved that SNOM and STM can be handled using the same formalism [64, 65]. It 
follows that a tip is able to detect high spatial frequencies if the tip itself creates high 
spatial frequencies when it is illuminated from the far field. More importantly, the 
detected field is given by an overlapping integral between the experimental field and 
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the response function of the probe, describing the spatial localization of the detection 
process, the polarization effects, and the spectral response [66-68]. 

In the experimental framework, nanostructures such as micro-disks, nanoparti- 
cles, nanoplasmonic resonators, and especially photonic crystal nanocavities represent 
successful platforms for investigating the nature of the signal detected by any SNOM 
probe, since they are able to localize e.m. fields in very small modal volumes [69]. It 
has been shown that, on one hand, the signal directly collected by dielectric tips or by 
metal coated circular aperture probes is the signature of the electric field intensity of 
the localized optical mode [70, 71]. On the other hand, aperture metal coated tips 
fabricated with a wavelength dependent specific metal thickness or with a split ring 
can directly probe a signal that is proportional to a specific spatial component of the 
magnetic field intensity of light [72-74]. Recently, unexpected results has been 
proved to overcome this dichotomy by showing that aluminum coated circular aper- 
ture probes directly detect a combination of all the in-plane electric and magnetic 
components, with relative weights that depend both on the probe geometry and on the 
probe-sample distance [75]. 


In order to improve the SNOM fidelity imaging of localized photonic modes and, 
at the same time, to circumvent the uncertainties related to the direct signal mapping, 
it has been demonstrated that an advantageous imaging technique relies on the tip- 
induced perturbation on the nanocavity environment. The strength of this method re- 
lies on the fact that, the presence of the probe itself is responsible of a spectral shift 
of the PCN localized modes, with respect to the unperturbed case, with a shift ampli- 
tude proportional to the mode e.m. field intensity [76-80]. As the presence of the tip 
above the cavity adiabatically changes the resonance wavelength without introducing 
significant additional scattering losses (assumptions confirmed by experimental ob- 
servations), the relative wavelength shift (64/4) is defined by the overlap integral 
between the unperturbed (Eo, Ho) and the effective (E, H) fields inside the perturba- 
tion volume Vp, as exactly derived from Maxwell’s equations [77, 81, 82]: 
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where the 3x3 matrix @@ and œin the coordinates basis (x, y, Z) are the electric 
and magnetic polarizability of the tip, respectively. The relation Eq. (2.7) predicts that 
the wavelength shift induced by a near-field tip is entirely set by the opto-geometric 
properties of the tip and the field distribution of the PCN mode. It also reveals that a 
tip exhibiting a dielectric polarizability perturbs only the electric field of the cavity 
while a tip with an effective magnetic polarizability interacts only with the magnetic 
field. Then, by mapping the spectral shift induced by a probe with an effective electric 
(magnetic) polarizability as a function of position it is possible to experimentally re- 
trieve the electric (magnetic) LDOS confined in the nanocavity for every attainable 
wavelength. 
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To describe the light-probe interaction and the effective tip polarizability the sub- 
wavelength tip apex can be approximated with an electric (p) or magnetic (m) dipole 
induced by the external fields [83]: 


p = a°E (2.8) 
m = a™H (2.9) 


Generally, in light-matter interaction the magnetic coupling is smaller than the 
electric coupling, in particular, at the macroscopic level the magnetic susceptibility of 
a material is generally about 104 the electric susceptibility [84]. Hence, light-matter 
coupling is mainly governed by the electric interaction mediated by the electric polar- 
izability. Therefore, for a tip apex composed by homogeneous media at optical wave- 
lengths the only non-negligible matrix is @®. However, it is important to stress that 
this behavior can change as a function of the probe geometry. 

For a dielectric probe the apex can be approximated as a sub-wavelength dielectric 
sphere. Due to the symmetry only the diagonal terms of the electric polarizability are 
non-vanishing and since all directions are equivalent we find af, = ayy = 
ag, meaning that an electric field oriented along a given axis cannot induce a dipole 
along a different one. Moreover, concerning the direct imaging only the induced di- 
poles that lie in the sample plane (xy plane) are coupled to the modes propagating in 
the fiber and therefore only TE-like modes can be directly detected, while the induced 
out-of-plane dipole does not emit radiation along z. Since the size of the dielectric tip 
apex can be very small, by investigating the perturbation effect induced by a° it is 
possible to retrieve the localized electric field distribution with an improved spatial 
resolution with respect to the direct signal detection that is affected by far-field radi- 
ation collected from lateral sidewalls. 

In particular, the insertion of a dielectric (uncoated) near-field probe provides a 
local increase of the dielectric constant de due to the tip positive effective electric 
polarizability. This induces a spectral red-shift (towards higher wavelengths) of the 
nanocavity resonances that is proportional to electric field intensity component, thus 
ensuring a sensitive imaging of the electric field intensity of the optical modes local- 
ized in PCNs [76, 78]. In the approximation of weakly perturbation (64/4 « 1) in- 
duced by a pure dielectric nanometric object placed at distance z from the PCN sur- 
face, the resulting spectral shift of a given mode is proportional to the electric field 
intensity. It follows by specifying Eq. (2.7) in the case where the in plane coordinates 
rjj are separated to the z axis and the unperturbed field Eg(r,,) is assumed constant 
over the small tip cross section [76]: 


2 
OME 12) ae 62/7) A -z/a (2.10) 
A 2V max|e(r)|Eo(r//)| | 
where the modal volume V is defined by Eq. (2.6) and d represents the exponen- 
tial decay length of the mode outside the slab in the z direction. The effective tip 
polarizability can be evaluated by assuming that only the very tip apex contributes to 
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the local dielectric perturbation, and that it can be approximated with a cylinder of 
volume Rd and dielectric constant £p, where Rp is the tip apex radius: 
ae =3 2 mRpd (2.11) 
Therefore, Eq. (2.10) points out that the spectral shift has a straight correspond- 
ence with the unperturbed mode distribution while it is inversely proportional to the 
unperturbed mode volume. Being for a dielectric material ep > 1, from Eq. (2.11) 
results that æ? is positive and hence induces a spectral shift towards lower wave- 
lengths. The dependence of V allows to intuitively explain the dielectric induced spec- 
tral shift towards higher wavelength, as a result of mode volume increase in the system 
composed by the PCN and the near-field tip, with respect to the PCN alone. In fact, 
the mode that is resonant to the tip-cavity system has an electric field intensity delo- 
calized also in the dielectric tip apex and experiences a higher effective refractive 
index that induces the mode wavelength to be red-shifted [4]. Moreover, Eq. (2.10) 
accounts for the exponential decrease of the spectral shift amplitude as a function of 
the tip-slab separation. By performing an experiment in which the dielectric probe is 
raster scanned at constant height above the PCN surface, Eq. (2.10) points out that the 
unperturbed electric field intensity distribution can be retrieved by mapping the tip 
induced spectral shift as a function of position. This result was for the first time theo- 
retically presented in ref. [76], as highlight in Fig. (2.5). 


In contrast, metallic coated tips with a circular symmetric aperture act as diamag- 
netic systems, inducing a mode blue-shift, which is found to be proportional to the 
intensity of the magnetic field component that is orthogonal to the PCN surface (TE- 
like mode) [79, 80]. In fact, due to the metallic probe geometry the aperture coated 
probe exhibits a negative magnetic polarizability that induces a magnetic response to 
the e.m. modes. By approximating the coated probe apex with a metallic hollow cyl- 
inder, a time-varying magnetic field oriented along the cylinder axis induces a current 
flowing in the metal, due to the Faraday’s law [80, 83]. This current, in turn, produces 
an induced magnetic field opposed to the total field inside the cylinder, hence letting 
the effective tip magnetic polarizability to be negative. In short, the probe geometry 
exhibits a non-vanishing magnetic dipole, given by Eq. (2.9), and its symmetry im- 
poses that the non-diagonal terms of a" are vanishing, while 4%} and ayy can be 
considered negligible with respect to q7} since the flux of the in-plane components of 
the magnetic field is negligible with respect to the flux of the out-of-plane component. 
The resulting perturbation of the PCN magnetic environment produces a sizeable 
blue-shift of the localized eigenmodes, which according to Eq. (2.7) accurately map 
the spatial distribution of the resonant magnetic field intensity oriented along z. For 
instance, the experiment described in ref. [79] shows that the spectral shift map per- 
formed with an Al-coated aperture probe on top of an L3 nanocavity and reported in 
Fig. (2.6 b), results in excellent agreement with the intensity of the magnetic field 
component orthogonal to the slab surface, shown in Fig. (2.6 c). 

In a simple and intuitive picture, the coated probe apex can be seen as a local 
metal mirror on top of the sample surface that reduces the PCN effective mode vol- 
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ume, thus shifting the resonant modes towards lower wavelengths [85]. It is also im- 
portant to note that a magnetic response of the constituent materials, namely glass and 
aluminum, is virtually absent and the effective magnetic response is only induced by 
the probe geometry. 


2.3 Disorder based nanocavities 


A small amount of disorder in a periodic photonic crystals would not destroy a 
photonic band-gap or localization light behavior. However, it would strongly change 
the properties of the localized modes, such as the resonant wavelength and the quality 
factor [86, 87]. Even in the limit of highly disordered dielectric structures light prop- 
agation can still be prevented and strong mode localization be achieved through the 
mechanism of Anderson localization [88-90]. 

Wave localization in random media occurs as a consequence of coherent construc- 
tive/destructive interference in multiple scattering and it plays a central role in the 
description of transport of quantum and classical waves. The physical mechanism be- 
hind localization, was introduced by Anderson in quantum electron transport [88]. 
The explanation is so general that the study of localization has not been restricted to 
electrons but light, sound waves, microwaves, and ultracold atoms have been consid- 
ered in theoretical and experimental investigations. Photons exhibit the advantage that 
no inter-particle interaction is present at moderate intensity (neglecting non-linear 
mechanisms), thus making photonic systems more suitable for studying Anderson lo- 
calization than electronic ones [91]. 

Disordered photonic systems are material in which the refractive index varies ran- 
domly in space. Different light mesoscopic regimes can be found as a function of the 
light wavelength (A), the sample size (L), and the photon scattering mean free path 
(lir), which quantifies the average distance travelled by the particle between succes- 
sive scattering events), as highlighted in Fig. (2.7 a-c). Neglecting absorption, the sin- 
gle scattering regime describes the case where the photon has a high probability to 
perform up to one scattering event before leaving the sample, exhibiting le ~ L, see 
Fig. (2.7 a). In general, the magnitude of the refractive index mismatch determines 
scattering strength. The ratio between 4 and the scatterer size gives the angular de- 
pendence of the scattering as well as the scattering cross-section in Mie’s theory [92, 
93]. For example, if the scatterer size is much smaller than A, the photon will be scat- 
tered nearly isotropically. On the other hand, if the size of the scatterer is comparable 
to A, the photon will be mainly scattered in the forward direction. When a disordered 
sample is illuminated by a coherent light source, the outgoing scattered light projected 
into a screen gives rise to a complex bright-dark intensity pattern, which is called laser 
speckle pattern. The analysis of the speckle pattern has developed a powerful tool for 
the study of the dynamic behavior of particles suspensions or macromolecules [94]. 

In another limit, for A « ler « L, the photon may scatter many times before leav- 
ing the disordered sample. This is the diffusion regime, which concerns a large variety 
of phenomena in opaque media, where multiple light scattering takes place, as re- 
ported in Fig. (2.7 b). Diffused photons undergo multiple scattering events and quickly 
loose the initial coherence before exiting the sample in a randomized direction. The 
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phase of the scattered waves after each scattering event is random, and thus one can 
model light propagation as a random-walk which leads to a multiple scattering pro- 
cess. But, the diffusive picture is not complete as some interference effects can survive 
the multiple scattering, as highlight by the two closed loop interference trajectories in 
Fig. (2.7 c). In fact, as the scatterers density increases and the condition (known as the 
modified Ioffe-Regel criterion [89]) lẹ < 4 holds, a transition from diffusion to lo- 
calization occurs for light, driven by multiple elastic scattering interference. In this 
case the interference effects are so strong that light transport comes to a halt and the 
waves becomes trapped in randomly distributed modes inside the sample. This phe- 
nomenon is called Anderson localization regime. Light behaves as being trapped in 
an optical nanoresonator. Therefore, disordered structure can sustain optical modes 
localized at the nanoscale. 

Of course, the main difficulty in the search for light localization relies on the re- 
alization of sufficiently strong scattering sample, in which the electric field can not 
even perform one oscillation before the wave is scattered for a second time. Consid- 
ering a semiconductor slab analogous to the one of Fig. (2.2 a), in which photonic 
crystals are designed, it is possible to achieve strong light localization by randomly 
engineering the air holes pattern distribution, see Fig. (2.8 a), as reported in the inten- 
sity distribution obtained in ref. [95]. This is also shown in Fig. (2.8 b), which exper- 
imentally demonstrates the presence of cavity-like confined optical modes. The two- 
dimensional random nanotextured slab has been proved to be a fertile system to study 
the localization regime, which is interesting also for other kinds of particles, because 
of the universality of this phenomenon in every random media wave transport. More- 
over, it accounts for necklace states occurrence and even for improving the absorption 
efficiency of thin film solar cells [95-98]. 


In Chapter (6) experimental near-field analysis of a random system on slab anal- 
ogous to the schematics reported in Fig. (2.8 a) will be presented. In particular, the 
light confinement will be engineered as a function of the design structural parameters 
and as a function of the post-fabrication oxidation technique that locally modifies the 
dielectric environment. 
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Figure 2.1: a) Schematics of the 1D photonic crystal along z direction with lattice 
constant a, for which the band diagram calculation reported in b) is performed. The 
details of the two media are: €, = 13 and layer thickness 0.2a, £, = 1 and layer thick- 
ness 0.8a, respectively. b) Energy band diagram calculated by PWE in the first Bril- 
louin zone for the infinite stack structure reported in a). In yellow are reported the 
photonic band gaps. The group velocity vg = dw/dk vanishes at the border of the 
Brillouin zone and also for k, = 0, with the exception of the lower energy band. 
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Figure 2.2: a) Schematics of the 2D photonic crystal of air holes arranged on a hex- 
agonal lattice with lattice constant a and filling fraction f = 35%, in a € = 12 slab 
with thickness 1.03a, whose band diagram calculated by PWE is reported in b). The 
Energy band diagram is reported as a function of the in plane wavevector in the first 
irreducible Brillouin zone (reported in the inset with the high symmetry points). Only 
the five lower energy band for the TE-like and TM-like modes are reported. The yel- 
low stripe represents the photonic band-gap that occurs for the TE-like modes only. 
The blue regions correspond to modes above the light cone that are not confined in 
the slab. 
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Figure 2.3: a) Representation of the unit-cell mesh used for FDTD calculations. It 
highlights the positions where local electric and magnetic field components are eval- 
uated. b) Simplified scheme involved in the FDTD calculation loop to evaluate the 
e.m. fields by solving the approximated Maxwell’s curl equations in which the deriv- 
atives become finite-differences. 
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Figure 2.4: a) FDTD calculated normalized spectrum (blue line) of the modified L3 
PCN on-slab reported in b) together the normalized emission of the electric dipoles 
placed inside the slab (black line). b) Top-view zoom of the 3D environment where 
the fields evolve according to the FDTD algorithm. The unit cell size of the discrete 
grid is 22 nm. Black holes represent air (n = 1) and orange environment simulates 
GaAs (n = 3.484). The withe dashed rectangle highlights the edges of the FDTD de- 
tector. c) Electric field intensity distribution calculated at the wavelength of the most 
intense resonant mode (A = 1285 nm). The scale bar is 600 nm. 
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Figure 2.5: a) Schematic of the PCN exhibiting a resonance at 2 = 1500 nm, whose 
surface is raster scanned by cylindrical silicon tip with a diameter of 125 nm that is 
placed at 30 nm above the slab. b) Electric field intensity distribution of the unper- 
turbed mode calculated by FDTD ina plane placed at 30 nm above the slab. The black 
circumferences outline the air pores. c) Calculated resonance frequency shift induced 
by the tip as a function of position by exploiting Eq. (2.10), Eq. (2.11) and FDTD 
calculations. Figures adapted from [76]. 
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Figure 2.6: a) Schematics of the Al-coated aperture probe. b) Measured spectral shift 
as a function of position in a PCN exhibiting a resonant mode at A = 1530 nm. The 
withe circles outline the air pores. c) FDTD calculation of the intensity ofthe magnetic 
field component orthogonal to the slab surface. Both maps are adapted from [79]. 
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Figure 2.7: Schematics of the three main light regimes in a disordered structure made 
of randomly placed dielectric spheres, for increasing scatterers density. The yellow 
spot is the light source and the red arrows describe the light random-walk trajectories. 
a) Single scattering regime takes place if the scattering mean free path is of the order 
of the sample size. b) Diffusive regime is found when multiple scattering occurs along 
with the condition À « ler. c) Light localization can take place if lẹ- < A due to strong 
interference. For instance, two trajectories propagating in opposite directions along a 
closed loop exhibit constructive interference and maximize the probability of locali- 
zation. 


Figure 2.8: a) Schematics of the randomly nanostructured semiconductor slab. b) Ex- 
perimental three-dimensional view of the near-field intensity distribution of a local- 
ized mode corresponding to the peak wavelength shown in the inset. Figure b) is 
adapted from [95]. 
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Chapter 3 


Deep sub-wavelength imaging 


Photonic and plasmonic devices rely on the nanoscale control of the spatial and 
energy distributions of the local density of optical states in dielectric and metallic 
environments. The tremendous progresses in designing and tailoring the electric 
LDOS of nano-resonators demand for a deep sub-wavelength imaging tool able to 
access the optical mode features, triggering the development of different nanoscale 
imaging techniques. By exploiting the near-field probe perturbation imaging tech- 
nique, on one hand, SNOM equipped with dielectric tapered probes provides a sub- 
wavelength electric LDOS detection; on the other hand, SNOM equipped with aper- 
ture metal coated probes provides the localized magnetic field component of light. 
The advent of artificial metamaterials operating at optical frequencies, in which the 
magnetic interaction with light can be as relevant as the electric response, makes it 
straightforward to simultaneously detect both field components. Here, we develop a 
plasmonic near-field probe to strongly enhance the collection efficiency from light- 
emitters and, more interestingly, to achieve a simultaneous electric and magnetic 
fields deep subwavelength detection. 


3.1 Experimental setup 


The scanning near-field optical microscope we used is a room-temperature com- 
mercial apparatus (TwinSNOM by Omicron), that is assembled on an Zeiss micro- 
scope conveniently modified to house the tip holder together with the scanning stages. 
The standard optical microscope is used for correctly positioning the probe on the 
sample surface but also to measure the far-field emission. The sample is mounted on 
piezoelectric stages (P-500 by Physik Instrumente) that enable a nominal spatial ac- 
curacy along x,y and z directions (where the xy plane is parallel to the sample sur- 
face) equal to (1, 1, 0.1) nm in the range (100, 100, 10) um. The scanning speed is 
kept sufficiently low to prevent artefacts caused by the limited photo detector response 
time. As we have discussed in Chapter (1), in order to achieve a sub-wavelength im- 
aging resolution two requirements have to be fulfilled: i) the probe has to be placed in 
the near-field of the sample; ii) the probe must have a sufficiently high throughput. 

In the experimental setup the optical fiber terminating in the sharp-tip is glued on 
an insulating module equipped with a piezoelectric bimorph, as reported in Fig. (3.1 a- 
b) following the scheme developed in [1]. Bimorph is made of two plates of piezoe- 
lectric material with opposite polarization vectors as shown in Fig. (3.1 c). If a voltage 
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is applied to bimorph electrodes, one of the plates will extend and the other will be 
compressed, resulting in a bend of the whole element. In order to control the probe- 
sample distance the excitation piezoelectric, layer 2 in Fig. (3.1 b), induces a small 
lateral tip oscillation parallel to the sample surface, with a frequency resonant to the 
system own mechanical frequency (about 60 KHz). While the detection piezoelectric, 
layer 1 in Fig. (3.1 b), measures the effective tip oscillation signal with amplitude and 
phase, in analogy with the tuning fork mechanism reported in [2]. As the probe ap- 
proaches the surface, due to probe-sample interacting shear forces that are effective 
only at few nanometers distances, the detected amplitude decreases and the relative 
phase between excitation and detection signals changes. In order to keep constant this 
phase difference, an electronic feedback drives the sample stage along z, thus ensuring 
a constant probe-sample distance throughout the whole scan. Finally, by reporting the 
signal that controls the z piezo as a function of position, we are able to map the surface 
topography resembling the AFM detection. Furthermore, condition i) is satisfied since 
by a proper calibration of the measured phase difference we can fix the tip-sample 
distance in a range between 5 nm and 20 nm, which corresponds to the near-field 
regime, if working with visible or infrared light. On the other hand, condition ii) re- 
quires a detailed characterization of the apex of the near-field probe. This SNOM 
setup allows to employ different probes as well as different configurations, such as 
illumination-collection or near-field collection with far-field illumination. Moreover, 
many laser sources are available in the visible and near-infrared spectral range in order 
to perform photoluminescence (PL) of optically active samples or resonant scattering 
measurements in samples in which light sources are not present. 

The schematic of the typical experimental setup performed to detect the PL signal 
that is coupled to the localized PCN modes, thus allowing us to address the optical 
properties of the nanocavities, is reported in Fig. (3.2). The SNOM is used in illumi- 
nation-collection geometry. The PL of the QDs is excited by a laser coupled into a 
given near-field probe. In the excitation path laser collimated light transmitted through 
an optical Faraday isolator (F), to prevent back-reflection that alters the laser stability, 
then through a simple telescope (T) formed by two lens, with f,/f; < 1 to reduce the 
beam spot size. Before being coupled to the optical fiber that ends with the near-field 
probe, it passes through a dichroic mirror (DM) that with high efficiency (up to 90%) 
selectively transmits visible light (up to 1.0 um) and reflects infrared light. The laser 
power emission is possibly adjusted by neutral density filters and the optical fiber 
used is a single mode fiber operating at 1300 nm. The PL is collected by the same 
probe, reflected by the dichroic hot mirror, dispersed by a spectrometer (SP, with a 
30 cm focal length equipped with a 600 lines/mm grating that exhibits 1.0 um blaze). 
Finally, it is detected by a liquid nitrogen-cooled InGaAs array detector (DTC, 
equipped with an array of 1024 photodiodes, each with a 25 um x 500 um size). The 
overall spectral resolution is 0.11 nm, as obtained from by the FWHM of the experi- 
mental spectral response to a delta-like CW laser emitting at A =1300 nm. At every 
tip position a PL spectrum is collected, thus obtaining a hyperspectral image of the 
PCN. 


62 


Niccolò Caselli 
3.2 Test-bed sample: the D2 photonic nanocavity 


The investigated 2D photonic crystal, together with the inserted nanocavity, is 
realized on a air-suspended membrane of Gallium Arsenide (GaAs) and it is com- 
posed by hexagonal lattice of air holes with f= 35% that has been presented in Chapter 
(2) [3, 4]. The fabrication process involves four steps schematically reported in 
Fig. (3.3 a). In the middle of the 320 nm thick membrane three layers of high-density 
(10!! cm?) Indium Arsenide (InAs) QDs are grown by molecular beam epitaxy, by 
exploiting the Stranski-Krastanov self-assembly technique [5]. The membrane is 
grown on top of a 1.5 um thick Alo.7Gao.3As sacrificial layer by patterning a 150 nm 
thick SiO2 mask with a 100 kV electron beam lithography and subsequent CHF3 
plasma etching. Then, the negative of the mask is transfer on the GaAs layer by 
SiCl4/O2/Ar reactive ion etching. The suspended membrane is finally released by a 
selective etching of the sacrificial layer in a Hydrofluoric acid (HF) solution. To ob- 
tain PL emission we illuminate the sample by the SNOM probe with a laser that en- 
sures sufficient light energy to excite the QDs exciton recombination, in particular a 
A =780 nm cw diode laser or a A =514 nm cw Argon laser are available. At room- 
temperature the PL of the large ensemble of quantum dots acts as a broadband internal 
light source that covers the spectral region in which the PCN modes are present 
(1200 nm and 1350 nm). The near-field PL spectrum averaged on a membrane region 
outside the photonic crystal (where the QDs emission is not reduced by the photonic 
band-gap) is reported in Fig. (3.3 b) and shows two broad peaks corresponding to the 
single QDs ground state and first excited state, respectively. When the photonic crys- 
tal structure is superimposed to the QDs layers, only the allowed narrow are clearly 
observed since a large amount of QDs emission is inhibited. The investigated photonic 
nanocavity is formed by removing 4 adjacent holes from the crystal lattice in dia- 
mond-like geometry, as highlighted by the SEM image of Fig. (3.4 a). This kind of 
PCN is called D2 nanocavity and has already been widely studied [6-10]. Due to the 
peculiar geometry, the local defect lowers the symmetry of the two-dimensional lat- 
tice and the PCN belongs to the Cov spatial symmetry point group and each PCN 
localized mode is attributed to one of the group irreducible representations [11-13]. 
Therefore, the D2 cavity supports non-degenerate resonant modes that exhibit differ- 
ent TE-like polarization with defined spatial symmetries in the localized e.m. field 
distribution [7, 8]. 

The typical SNOM PL spectrum averaged over the entire PCN and collected by a 
dielectric probe is reported in Fig. (3.4 b). It shows two narrow and spectrally sepa- 
rated peaks with Q = 3000 corresponding to non-degenerate resonant modes. It has 
been shown that the lower energy mode (M1) has a linear polarization along x while 
the first excited energy mode (M2) shows elliptical polarization with the dominant 
contribution parallel to y axis [7]. Different structure, where the photonic crystal lat- 
tice constant varies in the range 0.22 < a/A < 0.29 but with the same filling fraction, 
are realized in order to overcame the unavoidable fabrication induced disorder in the 
photonic crystal, which can largely shift the PCN resonant modes. These scaled struc- 
tures show the same optical properties apart from different resonant wavelengths, and 
we choose to investigate the one with the best wavelength matching between the QDs 
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emission and the PCN modes. Moreover, the same PCN is replicated several times to 
verify the fabrication reproducibility. 


3.2.1 Electric field intensity imaging 


By performing SNOM measurements it is a straightforward result to achieve the 
intensity map of every resonant mode on the sample surface. In Fig. (3.6 a-b) are 
plotted the experimental intensity distributions for the M1 and M2 mode, respectively, 
centered on the D2 nanocavity and evaluated by reporting as a function of position 
(i.e. for each spectrum collected) the intensity corresponding to the maximum wave- 
length of each mode. In particular, the M1 mode is found to be elongated along the y 
axis with three circular equally spaced high intensity regions, while M2 shows two 
lateral crescent shapes distributed along x direction and a dark central region. Even- 
tually, both intensity distributions have even parity along both axes. The choice of 
plotting the maximum wavelength maximizes the map intensity contrast. 

As explained in Chapter (2) a sophisticated method to improve the SNOM spatial 
resolution and mode fidelity imaging, without performing further measurements, re- 
lies in the study of the tip-induced perturbation imaging and the accurate analysis of 
the localized modes spectral shift (64) as a function of position. In the case of dielec- 
tric near-field probes this spectral shift corresponds to the electric field intensity, i.e. 
the electric LDOS, as stated in Eq. (2.10). The principle underlying the spectral shift 
evaluation is briefly explained in Fig. (3.5). In Fig. (3.6 c-d) are shown, for both res- 
onant modes, the spectral shift amplitude maps evaluated by performing a Lorentzian 
function fit of each spectrum in the peak wavelength range and then by extrapolating 
the spectral shift as a function of position following the method reported in Fig. (3.5). 
Hence, the blue regions of Fig. (3.6 c-d) represent a wavelength difference equal to 
Zero. 

The high fidelity with which the spectral shift maps reproduce the electric LDOS 
of the nanocavity localized modes is demonstrated by the very good agreement with 
the three-dimensional FDTD calculated electric LDOS reported in Fig. (3.6 e-f) for 
both M1 and M2 with unit-cell edge length equal to 10 nm and a spectral resolution 
of 0.03 nm. The spectral shift distributions clearly show more and smaller details with 
respect to the intensity maps. In the case of MI, Fig. (3.6 c), are present five intensity 
lobes along y that exhibit the maximum spectral shift in correspondence of four sym- 
metric hot-spots. In order to evaluate the SNOM imaging spatial resolution we com- 
pare the normalized vertical profiles centered on the experimental and theoretical 
maps of the mode M1, as reported in Fig. (3.7). In particular, in Fig. (3.7 b) on one 
hand, the row PL intensity data do not reproduce every spatial oscillation of the elec- 
tric field intensity as calculated through FDTD, resulting in a spatial resolution insuf- 
ficient to detect all the relevant mode features. On the other hand, the spectral shift 
profile nicely shows the same modulation of the theoretical profile. However, to figure 
out the actual experimental spatial resolution it is not correct to evaluate the smallest 
spatial feature of the map. The reason is that the investigated PCN localized modes 
do not show delta-like spatial features, but the smallest map lobe calculated by FDTD 
is of the order of 100 nm. Therefore, to estimate the spatial resolution we compare the 
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experimental spectral shift profile to the profile obtained from the convolution be- 
tween the calculate LDOS map and two-dimensional Gaussian point-spread functions, 
characterized by different full widths at half maximum (FWHM). In Fig. (3.7 c) are 
reported the experimental data and the theoretical profile obtained by a convolution 
with a Gaussian function of 70 nm FWHM, that minimizes the y? distribution result- 
ing to a reduced-y?of 1.1. Furthermore, it is important to notice that since the resonant 
wavelength are extrapolated from Lorentzian fits, the resulting uncertain is lower than 
the experimental spectral resolution and it is directly provided by the fitting algorithm 
for each spectrum, as reported in the error bars of Fig. (3.7 c). Therefore, the SNOM 
PL, combined with the tip-induced spectral shift analysis, allows us to obtain high 
fidelity electric LDOS mode imaging with an ultra-subwavelength spatial resolution 
of 70 nm, comparable to the so far best result obtained on PCNs [14]. 

By the way, the reported spectral shift measurements are in principle analogues 
to the result obtained in ref. [6] but with an improved spatial resolution, down to 4/19 
where A is the vacuum wavelength of the localized mode. Moreover, the reported 
near-field PL experiments are performed at sufficiently low power laser excitation 
(60 uW) to avoid the presence of a spectral red-shift induced by local high-power 
laser heating, that may be superimposed or even comparable to the tip-induced spec- 
tral shift [9]. 

The dielectric tip induced perturbation is referred to as “gentle” perturbation 
meaning that an observable spectral shift is achievable without a large degradation of 
the cavity Q factor, thus preserving the resonant mode features such as the electric 
field spatial distribution. This assumption is confirmed by far-field confocal experi- 
ments where the SNOM probe in Fig. (3.2) is replaced by a standard 50X objective 
and a single mode optical fiber acts as a pinhole, that gives the estimation of the actual 
unperturbed nanocavity Q [8, 15]. In particular, for the same nanocavity investigated 
in Fig. (3.4) the bare PCN possesses Q = 3100 and the dielectric probe locally in- 
duces a Q decrease that ranges in between 5% (when the tip is positioned on the mode 
tail) and 20% (for tip positioned on the mode electric field maximum) with respect to 
the unperturbed value. 

As in any real experiment, the reported SNOM perturbation imaging technique 
has also some limitations. To address the issue of reproducibility the same measure- 
ments are repeated several times also with different dielectric tips. We find that not 
every probe gives the same result and a selection of good tips is needed. Only about 
25% of the total number of dielectric probes that we fabricate by tube chemical etch- 
ing gives reproducible results to map the electric LDOS by inducing a gentle and de- 
tectable spectral shift. Concerning the imaging performed by a given tip, we obtain 
similar spectral shift maps and analogue values of maximum spectral shift inde- 
pendently on the investigated nanocavity. However, every near-field probe tends to 
be damaged by the use, exhibiting a lifetime of the order of a couple of weeks of 
continuous work. Moreover, since the spectral shift is inversely proportional to the 
unperturbed resonant mode volume, the perturbation imaging technique could be ef- 
ficiently applied to any optical nanoresonator in which the perturbing volume is not 
negligible with respect to the localized cavity mode volume, as stated in [16]. For 
instance, in the case of M1 mode in the D2 nanocavity we find by accurate FDTD 
calculation (by employing many large detectors parallel to the slab plane at different 
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height both inside and outside the slab) and exploiting Eq. (2.6) a mode volume V = 
0.033 um}, which is comparable to the L3 nanocavities mode volume [17]. In addi- 
tion, we are able to give an estimation of the effective tip apex size that is responsible 
of the dielectric perturbation, approximated with a glass sphere of radius Rp and 1.5 
refractive index. We use Eq. (2.10) and Eq. (2.11) to link the measured maximum 
spectral shift to FDTD parameters such as the nanocavity mode volume and the out 
of slab electric field intensity decay rate (d=40 nm, calculated employing multiple 
field detectors outside the slab). Then, we obtain the tip effective electric polarizabil- 
ity a a° =1.1-107! mie and the radius size Rp = 100 nm. This vale of Rp is compa- 
rable to the estimation evaluated by SEM image of the tip. Therefore, the assumption 
of a tiny but non negligible perturbation is correct, since the perturbing volume Vp = 
tR8d is about 30% of the mode volume V. Finally, it is worth stressing that the esti- 
mation of the experimental electric LDOS imaging spatial resolution, 70 nm reported 
in Fig. (3.7 c), is in agreement (being of the same order of magnitude) with the esti- 
mation of the effective tip apex radius that induces the observed spectral shift. 


3.2.2 Magnetic field intensity imaging 


The capability of the SNOM imaging can be enlarged towards the detection of the 
magnetic component of light localized in a PCN if the dielectric probe is replaced 
with a metal (aluminum) coated aperture probe. In fact, as discussed in Chapter (2), 
the aperture probe placed on top of a nanocavity acts as a conductive ring in which 
the PCN localized magnetic field induces an opposed magnetic dipole, see Fig. (3.8 a- 
b). The perturbation effect induces a resonant mode blue-shift (towards lower wave- 
lengths) that is proportional to the intensity of the magnetic field component orthog- 
onal to the PCN plane (H2), as stated in Eq. (2.7). The magnetic intensity can be 
retrieved with a sub-wavelength spatial resolution by performing the PL measure- 
ments and spectral shift analysis employing metal coated tips. Therefore, we have the 
capability to image a quantity that in general interacts very weakly with ordinary mat- 
ter. The first reported observations has been performed simultaneously in ref.s [10, 
18]. The former work, investigating a L3 nanocavity, found an effective negative mag- 
netic polarizability for the aperture tip, equal to @?} =12-107! m3/y9. The main re- 
sults of ref. [10] are reported in Fig. (3.8 c-h), since they concern both M1 and M2 
modes of the D2 nanocavity. The collected PL intensity maps, Fig. (3.8 c-d), repro- 
duce the basic features of the localized electric field intensity for both modes, with a 
100 nm spatial resolution limited by the probe aperture size. However, normalizing 
for the incident laser power used in experiments with dielectric probes, a huge de- 
crease of the collected signal, down to 10°, is found. This requires time-consuming 
measurements with coated probes. In addiction, by reconstructing the spectral shift 
amplitude maps, Fig. (3.8 e-f), are obtained the images of H? for both modes in agree- 
ment with the FDTD calculations reported in Fig. (3.8 g-h) and, with a sub-wave- 
length spatial resolution. However, it is clear that the experimental spectral blue-shift 
maps are affected by strong noise. This behavior is due to the fact that in the regions 
of maximum blue-shift the PL signal is minimum. In fact, for a spatially confined 
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mode the magnetic field has a local maximum (proportional to the spectral shift am- 
plitude) in the positions where the electric field (proportional to the signal intensity) 
is minimum and vice-versa. 

In conclusion, the perturbation imaging allows, to map the PCN localized electric 
field intensity, by employing dielectric near-field probes, and also to retrieve the in- 
tensity of the localized magnetic field component that is orthogonal to the photonic 
slab, by using metal coated aperture probes. 


3.3 Campanile tips for near-field probing 


Given the electric and the magnetic field imaging, it is tempting to combine the 
two spectrally opposite dielectric and metallic perturbations by employing a probe 
featuring both an electric and a magnetic effective polarizability. A promising ap- 
proach can be found in near-field probes functionalized with plasmonic based nano- 
antenna, as described in Chapter (1). In fact, by designing an advanced metal-dielec- 
tric structure on top of a dielectric tip we will show how it is possible to induce a 
simultaneous electric and magnetic interaction that are comparable in strength. More- 
over, this plasmonic probe can be easily employed in the SNOM illumination /collec- 
tion configuration, since it allows for an efficient far-field to near-field bidirectional 
coupling over a wide spectral range with an outstanding sub-wavelength spatial reso- 
lution. 

The geometry of the presented plasmonic SNOM probe consists of a square-based 
three-dimensional tapered structure terminating in a metal-insulator-metal (MIM) 
nanogap, which resembles the shape of a campanile bell tower, hereafter referred to 
as “campanile tip”. The campanile tip design is reported in Fig. (3.9 a). This structure 
merges high e.m. field enhancement and confinement at the nanoscale with a broad- 
band operation, since it exploits adiabatic plasmonic compression [19-22]. The fab- 
rication process is performed at the apex of a glass-core optical fiber (single mode for 
A =633 nm, with core diameter size equal to 2.5 um ) tapered by the tube etching 
technique described in Chapter (1). Using a Gallium focused ion beam (FIB) milling 
(Zeiss Crossbeam 1540), the campanile geometry is carved into the etched glass fiber 
tip. Using shadow evaporation, 50 nm thick Au layer is evaporated on two opposing 
sides of the campanile structure. Typically, the resolution of the angle resolved evap- 
oration is not sufficient to create a well-defined gap between the two triangular plates 
at the tip apex. Hence, a gap is subsequently cut by FIB at the tip apex, resulting in a 
40-50 nm wide and 30 nm long gap. The SEM images of the campanile tip apex that 
exhibits a clear nanogap formation are reported in Fig. (3.9 b-c). Of course, the fabri- 
cation resolution can be even increased by using a Helium FIB. The campanile tip 
geometry allows strong electric field enhancement confined at the nanogap, as shown 
by the finite element calculation performed by COMSOL Multiphysics 4.2 software 
reported in Fig. (3.10). The calculation are evaluated at 2 =666 nm for a nanogap size 
of 10 nm. The other design parameters are reported in the caption of Fig. (3.9), and 
the frequency dependent relative permittivity of gold is taken from Palik’s handbook 
[23] while the refractive index of the SiO2 layer is set as n = 1.5. The electric field 
enhancement distribution inside the campanile tip, that is the ratio of the electric field 
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strength, |E|?, in the nanogap and the incoming field strength, |Eo|7, is reported in 
Fig. (3.10 a-b). It shows that the probe provides efficient coupling between far-field 
and near-field. In fact, the fundamental mode in an MIM structure, polarized along 
the direction perpendicular to the metal surface (y), is supported without any cutoff 
frequency, no matter how thin the insulating layer [20]. The plasmonic mode in the 
campanile is confined to the gap region, which thus defines the improved spatial res- 
olution with respect to conventional metal coated aperture probes. The field enhance- 
ment is plotted in Fig. (3.10 c) for an aperture size of 10 nm. Typically, in conven- 
tional aperture probe, a 30-100 nm thick metallic cladding is employed to confine 
light inside the fiber and funnel it to the nanoscale aperture, as described in Fig. 
(1.8 b). However, this metal coating implies boundary conditions for the e.m. wave 
travelling inside the fiber leading to a mode cutoff, including the SPP mode, before 
reaching the aperture. The strong field confinement in the campanile tip is due to the 
smooth reduction of the radius of the metal tapered waveguide that produces an adia- 
batic transformation of the SPPs. Phase and group velocity of the SPPs tend to zero 
while approaching the tip, and at the same time the SPPs e.m. fields are largely en- 
hanced, thus producing a nano-focusing effect. Moreover, the high throughput is ef- 
fective over a broad spectral range, as highlighted in Fig. (3.10 d). The operating 
bandwidth is limited only by metal absorption at short wavelengths and extends into 
the infra-red region. The calculations of Fig. (3.10 d) also point out that the electric 
field enhancement is larger than the one of a bow-tie antenna with same-sized gap. 
Here, we evaluate only the electric enhancement, which corresponds to the tip ability 
to detect a signal in PL measurements, since we are interested in the campanile tip 
throughput. While the tip induced perturbation effect will be afterwards analyzed. An- 
other advantage of the campanile tip relies on the highly confined field in the gap 
region and results in almost background free acquisition. By operating with the ex- 
perimental setup presented in Fig. (3.2), the PL signal collected back through the an- 
tenna gap is affected by the background arising from the edge-scattered light, which 
is found to be far below the noise threshold. 


3.3.1 Mapping local charge recombination heterogeneity 


In order to demonstrate the notable performances of the plasmonic nanoantenna, 
we employ campanile tips with 40 nm gap for mapping the inhomogeneous radiative 
recombination in individual InP quantum-wires. These nanostructures are chosen be- 
cause of the intense and linearly polarized PL emission in the near-infrared and the 
potential application as light-harvesting nanostructures, due to the 1.4 eV bandgap 
and low surface recombination rates [24, 25]. 

The indium phosphide nano-wires are prepared by chemical vapor deposition us- 
ing trimethylindium and tertbutylphosphine as In and P sources, respectively. The 
nanowires are grown at 415 C for 10 minutes and under a flow of H2, on a Si [111] 
substrate with a 2 nm thick gold film acting as catalyst for the nano-wire formation. 
The process gives 1-5 m long and 30-100 nm thick nanowires. TEM analysis shows 
that the wires have a zinc-blend structure and grow in the [111] direction with cylin- 
drical shapes with clean, often atomically flat, cleaved ends. The nano-wires are then 
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dried and transferred to 100 nm thick SiO2 support with alignment marks to allow us 
to investigate the same wire with multiple techniques. A SEM image of one of the 
examined wires is reported in Fig. (3.11 e). 

The single nano-wire PL is collected by campanile tips employing the experi- 
mental setup of Fig. (3.2) with a He-Ne laser excitation, with a 50/50 beam splitter 
for replacing the dichroic mirror, with a 695 nm long pass filter in front of the spec- 
trometer and a Si-based CCD camera for detection at wavelength around 850 nm, with 
a spectral resolution of 0.10 nm. The light polarization is controlled during the prop- 
agation along the optical fiber by using mechanical compensator based on a Babinet- 
Soleil compensator, that permits to apply a controlled pressure and rotation to the 
fiber. Therefore, acting on the compensator, by maximizing the PL signal ofthe nano- 
wire, the excitation light is forced to be polarized along y when it approaches the 
campanile tip. At each pixel during the scan a spectrum is recorded, thus hyperspectral 
imaging is obtained as shown in Fig. (3.11 a-c). In Fig. (3.11 d and e) the same InP 
quantum-wire is investigated by a confocal setup, with a 100x and 0.7 NA objective, 
to highlight a comparison with the results obtained by campanile tip. 

In Fig. (3.11 c) normalized spectra acquired along the wire length by campanile 
tip are reported. They show evident local spectral variations. In particular, the hot 
spots of the intensity map in Fig. (3.11 a) shows spectra with a blue-shift emission of 
the band-edge, as well as stronger contributions from trap-related spectral components 
above the bandgap that broadens the spectra considerably. In fact, the emission peak 
at 839 nm (1.47 eV) corresponds to the expected 100-meV blue-shift relative to bulk 
InP wire, independent on quantum confinement [26]. 

In contrast, confocal intensity map of the same wire evaluated at the same wave- 
length, reported in Fig. (3.11 d), displays two maxima but without any spectral varia- 
tion along the wire, see Fig. (3.11 f). It was previously observed that strong PL en- 
hancements and spectral blue-shifts result from passivated InP nano-wire sur- 
faces [26]. This was attributed to Coulombic interactions between excitons and posi- 
tively charged trap states on the wire surfaces, that influence the optoelectronic prop- 
erties of nanocrystals and nano-wires [27]. The exciton diffusion length in these ma- 
terials is hundreds of nanometers, and therefore individual trap states within the dif- 
fusion volume should strongly influence the local absorption energy and charge re- 
combination rate. 

We believe that the heterogeneity observed along the wires, Fig. (3.11 a), on 
length scales well below the exciton diffusion length, are direct maps of trap-state 
modifications of the local exciton properties. The observed hot spots are probably due 
to an increase of trap-state densities (and changes in the native oxide layer) at the wire 
ends. Their spectral characteristics are consistent with a trap-induced Stark shift, pre- 
dicted to be about 60 to 70 nm above the band-edge for positive trap states [26, 27]. 
The absence of spectral variations in the confocal measurements is attributed to the 
low spatial resolution obtainable in the far-field, especially in the z direction that force 
to probe the entire wire thickness. Then surface-specific effects are obscured by the 
bulk behaviour. A different experimental approach based on cathodoluminesence 
could achieve a comparable spatial resolution. However, the large number of incident 
electrons fill the trap states and do not detect any spatial variation in emission from 
InP nano-wires [28]. 
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In order to test the imaging reproducibility, we investigated a second quantum- 
wire, whose SEM image is reported in Fig. (3.12 a), by a different campanile tip. The 
PL intensity distribution evaluated at every wavelength in the peak emission shows 
the same heterogeneity of the previous Fig. (3.11 a). The campanile tip provides an 
optical resolution much higher than the confocal approach, as can be seen in Fig. (3.12 
d) by comparing the intensity profile along the line orthogonal to the wire length ac- 
quired by campanile tip, by confocal setup and by SEM topography, highlighted by 
the red lines in Fig. (3.12 a, b and c) respectively. A line scan along the wire length, 
reported in Fig. (3.12 e), reveals a spatial resolution of 40 nm approximately equal to 
the nanogap size, whereas the topographic line scan is almost flat and shows negligi- 
ble variations. Here, the spatial resolution is evaluated as the spatial difference be- 
tween 90% to 10% of the PL signal in the measured sharpest variation, since, differ- 
ently from PCNs, we cannot model or calculate the emission properties of the inves- 
tigated quantum-wires. 

In conclusion, we prove that campanile-style far-field to near-field transformers 
provide a pathway for understanding energy conversion processes at critical length 
scales, yielding insights into the role of local trap states in radiative charge recombi- 
nation in InP nano-wires. More generally, our study demonstrates the impact of the 
campanile geometry on a wide range of nano-optical measurements, because it can 
virtually be implemented in many imaging and spectroscopy methods, including Ra- 
man and infrared/Fourier transform, as well as white-light nano-ellipsometry/interfer- 
ometric mapping of dielectric functions. We expect that the large bandwidth enhance- 
ment makes them ideal for ultrafast pump-probe and nonlinear experiments down to 
molecular length scale [29, 30]. They also could be used for ultrasensitive medical 
detection, photo-catalysis and quantum-optics investigations, as plasmonic optome- 
chanics and circuitry elements. Finally, concerning the measurements presented in 
this section, seven campanile tips on different optical fibers were produced. We tested 
all of them and they gave reproducible results. 


3.3.2 Simultaneous imaging of electric and magnetic fields 


In order to demonstrate that campanile tips are able to simultaneously image both 
the electric and the magnetic field intensity of light localized in optical nanoresona- 
tors, we perform PL measurements on the D2 photonic crystal nanocavity presented 
in Fig. (3.4). We employ the experimental setup of Fig. (3.2). As reported in the sche- 
matics of Fig. (3.13 a), the measurements are performed exciting the QDs embedded 
inside the sample with the cw diode laser (780 nm) and a power of 120 uW delivered 
to the near-field probe. The typical SNOM spectrum acquired in this configuration in 
the middle of the PCN with lattice constant a =311 nm is shown in Fig. (3.13 b), in 
the spectral range of M1 and M2 modes. It highlights a high signal/noise ratio, signa- 
ture of an excellent tip throughput with respect to standard coated aperture probes. In 
fact, if we compare the signal detected in Fig. (3.4 b) with the one reported in Fig. 
(3.13 b) on the same PCN, we find only one order of magnitude reduction in the 
throughput achieved by campanile tips with respect to dielectric tapered probes. 
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The high throughput is unfolded by the strong electric field enhancement confined 
at the nanogap, together with a magnetic field enhancement inside the tip in close 
proximity to the apex. These results are shown by the FEM calculation, evaluated at 
A = 1.3um for a 40 nm nanogap, reported in Fig. (3.13 c-d), respectively, in close 
analogy with the behavior reported in Fig. (3.10 b), which is evaluated at smaller 
wavelength and for a smaller nanogap. Moreover, we demonstrate that the high tip 
throughput is effective in a broad spectral range as highlighted in Fig. (3.13 e). Con- 
cerning the induced perturbation on the PCN, the metal-insulator-metal (MIM) apex 
of the campanile tip can be modelled by two metal dipole nanoantennas spaced by 
40 nm of dielectric material. Mimicking this kind of system by an equivalent RLC 
circuit is a widely used approach [31]. It is well-known that RLC circuits show an 
effective electric interaction associated with the capacity C together with an effective 
magnetic interaction that is driven by the inductance L. In fact, the electric polariza- 
bility describes dielectric probes as effective capacitors, while the diamagnetic re- 
sponse allows schematizing metallic coated aperture tips as inductances. The interac- 
tion of the campanile tip apex with AC electromagnetic fields can be mimicked by an 
equivalent RLC circuit. This analogy helps us to understand the simultaneous pres- 
ence of both electric and magnetic polarizability. In particular, in the campanile tip 
the inductance is driven by a combination of electric and displacement currents flow- 
ing at the MIM apex structure. Moreover, around the resonant wavelength the electric 
and magnetic energies stored in the circuit are of the same order of magnitude [32]. 
The key point we exploit is that, while the induced electric dipole is parallel to the 
external electric field, the L response is diamagnetic, meaning that the associated mag- 
netic dipole is antiparallel to the inducing magnetic field. This has strong effects on 
the tip perturbation of the nanocavity mode wavelengths, since the electric induced 
shift is expected to be toward red spectral regions, while the magnetic perturbation 
induces a blue-shift. Therefore, the role of the campanile tip is twofold. On one hand, 
it combines the dielectric perturbation effect with the perturbation induced by a metal 
coated tip, resulting in a concomitant electric and magnetic polarizability. On the other 
hand, the electromagnetic field enhancement accounts for PL spectra with a sig- 
nal/noise ratio that allows monitoring with high accuracy the spectral position of the 
resonant mode maximum. We evaluate the PCN localized modes spectral shift that is 
related to the electric and/or magnetic field component of light, depending on the ef- 
fective electric and/or magnetic interaction. 

By exploiting the perturbation approach, we experimentally succeed, during the 
same scan, in extracting high quality maps of the electric and the magnetic field in- 
tensities of the localized modes with a sub-wavelength spatial resolution. This allows 
us to overcome the long-standing problem, discussed in Sec. (2.2), that is related to 
the critical interpretation of the signal directly collected by SNOM probes. This issue 
normally requires an accurate model of the experimental configuration and compari- 
son of the collected signal with the calculated spatial distribution of the electric and 
magnetic fields, thus differentiating our work with respect to ref. [33]. 

We find that the peak wavelength of each mode is strongly dependent on the tip 
position on the sample surface. The effect is highlighted for M1 in Fig. (3.14 a-b), 
which report the schematic of the sample where three different spots (A, B and C) are 
indicated, together with the normalized PL spectra collected at these positions by the 
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campanile tip, respectively. The PL spectrum recorded at position A, where M1 mode 
has maximum electric field as stated in Fig. (3.6 e), shows a red-shift with respect to 
the spectrum detected outside the PCN (position C, where the intensity of the mode 
MI is quite negligible).At the same time the PL peak observed when the tip is posi- 
tioned in B is blue-shifted with respect to the spectrum detected in C. In fact, in B the 
cavity resonant magnetic field intensity has a local maximum, as highlighted by the 
dominant out-of-plane intensity component calculated in Fig. (3.8 g). Since the tip 
perturbation can be considered negligible a few microns outside the nanocavity, the 
observed red and blue spectral shift demonstrates the simultaneous electric and mag- 
netic susceptibility of the campanile near-field probe. In order to establish the precise 
resonant wavelength of both modes, Lorentzian fits of each spectrum are performed. 
To highlight the variation of the shift, we report in a two-dimensional map the wave- 
length difference between the peak wavelength at every position and the peak wave- 
length outside the PCN. The results are shown in Fig. (3.14 c-d), and correspond to 
the spectral shift maps of the M1 and M2 modes, respectively. Notably, a significant 
modulation of the spectral shift is observed in both maps. Of course, a useful charac- 
teristic of the investigated PCN in discriminating the tip induced shift of both modes 
is that they are more spectrally separated than the maximum observed spectral shift. 
The experimental ratio between the magnetic and electric response of the campanile 
tip can be evaluated as the ratio p between the maximum blue-shift and the maximum 
red-shift. For both modes it is found to be of the order of 0.5 and in particular: 


B(M1) = |--0.28 nm] / 0.60 nm = 0.47 (3.1) 
B(M2) = |—0.18 nm] / 0.33 nm = 0.55 (3.2) 


which demonstrate a well-balanced dielectric and diamagnetic response of the 
campanile probe. This characteristic gives simultaneous access, in a single measure- 
ment, to both fields with similar sensitivity, which is intriguing if one considers that 
the diamagnetic response of normally occurring materials in the optical spectral range 
is usually many orders of magnitude smaller than the dielectric response [34]. 

The experimental data show that the campanile probe has a concomitant positive 
electric and negative magnetic polarizability, as initially supposed. The different sign 
of the dielectric and diamagnetic response, together with the fact that for any standing 
wave the nodes of the electric field correspond to the maxima of the magnetic field 
and vice-versa, allow us to discriminate with high fidelity between the two field com- 
ponents. In principle, similar properties could also apply to other plasmonic nanoan- 
tennas. The use of a campanile probe has a double benefit, as it enhances the near field 
coupling and also makes the magnetic and electric perturbation interaction of compa- 
rable magnitude. 

To corroborate this interpretation we perform FDTD calculations of the field dis- 
tributions of the localized photonic modes. The maps associated to M1 and M2 are 
reported in Fig. (3.15 a-b) for the electric field intensity (red color scale) and in 
Fig. (3.15 c-d) for the magnetic field intensity in negative units (blue color scale), 
respectively. In order to compare the numerical calculations and the experimental data 
reported in Fig. (3.14 c-d), the experimental parameters p between the magnetic and 
the electric response of the campanile tip evaluated in Eq. (3.1) and Eq. (3.2) are used 
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to reproduce the expected spectral shift map. In particular, we construct a linear com- 
bination of the FDTD calculated electric and magnetic normalized intensities, assum- 
ing a negligible electric field at the magnetic field maxima and conversely, by using 
a simple relation formula deduced from Eq. (2.7): 


5A EMI lHo(r) 
A di Eo) lFnax Li p Ho) lFnax GS) 


The maps calculated with 8(M1) = 0.47 and B(M2) = 0.55 for the modes M1 
and M2 are reported in Fig. (3.15 e-f), respectively, and closely match the experimen- 
tally measured field patterns of Fig. (3.14 c-d) in all the relevant details. 

The experimental spatial resolution of the perturbation imaging method obtained 
with the campanile tip is estimated by analyzing the map along one-dimensional pro- 
file cut. For instance, in Fig. (3.16 a-b) are shown the M2 experimental spectral shift 
map and the FDTD calculation using Eq. (3.3), respectively, where the central hori- 
zontal dashed lines represent the two spatial profiles to be compared. The nanogap 
size at the apex of the campanile probe is a priori expected to give a spatial resolution 
of 40 nm. This assumption is examined by comparing the spectral shift data profile 
and the theoretical profiles obtained by the convolution between the FDTD map of 
Fig. (3.16 b) and two-dimensional Gaussian point spread functions, characterized by 
different FWHM and whose intensities are normalized to the experimental spectral 
shift amplitude range. The profile obtained with a Gaussian point spread function of 
40 nm FWHM nicely reproduces the measured values within the experimental uncer- 
tainties, by minimizing the y? distribution. In Fig. (3.16 c) we report the spectral shift 
data profile along with the FDTD results convoluted with Gaussian functions of 40 nm 
and 80 nm FWHM. The former convolution better reproduces the experimental data. 
Therefore, the spatial resolution achieved by the campanile probe imaging technique 
reaches the expected 40 nm limit (down toA/30), determining a huge improvement 
with respect to standard dielectric or metal coated probes. 

There are many advantages with respect to direct near-field signal collection im- 
aging performed with standard coated tip. First of all, the technique is sensitive not 
only to the in-plane e.m. fields, thus neglecting the dominant contribution to the mag- 
netic component of the TE-like mode, but it can access the total intensity of both 
fields, which correspond to the electric and magnetic LDOS of the localized modes. 
Secondly, the mode perturbation induced by the tip polarizabilities is directly propor- 
tional to the LDOS and does not require a detailed model of the near-field probe signal 
collection. This not only means a high fidelity imaging but also that after calibrating 
the mode shift on a well characterized sample it is possible to obtain a quantitative 
measurement of the electric and magnetic LDOS for any nanoresonators. 

Finally, and importantly for many applications, the campanile tip has an enhanced 
throughput, which leads to the possibility of investigating also weak resonances and 
even individual quantum emitters. For instance in the reported measurements the col- 
lected PL signal with campanile tip is three orders of magnitude higher than the PL 
collected in nominally identical experiments performed with Al-coated aperture 
probes with a nominal aperture of 200 nm, normalized at the same excitation power 
and integration time. The excellent agreement between the experimental data and the 
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near-field numerical calculations clearly indicates that the campanile probe induced- 
perturbation imaging is able to catch the details of the electric and the magnetic inten- 
sity distributions of optical modes, with an ultra sub-wavelength spatial resolution. 

Combined access to the near-field electric and magnetic components of optical 
modes, could improve innovative devices based on the interplay between the electric 
and the magnetic optical response, such as nano-metamaterials at optical frequencies. 
These nanoscale architectures, where the magnetic interaction with light can be as 
relevant as the electric response, exhibit unprecedented light behaviors, such as neg- 
ative or zero refractive index, optical cloaking and superluminal phenomena [35, 36]. 
Furthermore, they are promising both for fundamental science and for cutting edge 
technologies, triggering the quest for novel sensing applications and field-enhance- 
ment effects at the nanoscale [37, 38]. One possible actual implementation could be 
coupling light nano-emitters (like molecules or quantum dots) to metamaterial or 
nanoresonator for tailoring the spontaneous emission rate through electric and/or mag- 
netic dipoles interaction (electric or magnetic Purcell effect [39, 40]), by exploiting 
the knowledge of electric and magnetic local density of states with a unprecedented 
spatial resolution. 
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Figure 3.1: a) Image of the optical fiber that ends with the SNOM tip. The fiber is 
glued to the module that possesses piezoelectric controller whose electric contacts are 
clearly visible. b) Lateral view of the module schematics to highlight the bimorph 
piezoelectric (orange) whose detection and the excitation components are labelled as 
1 and 2, respectively. Black arrows show the tip controlled oscillation plane. c) Sche- 
matics of the piezoelectric bimorph operation, where E and P represent the electric 
field and polarization, respectively. 
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Figure 3.2: Schematics of the illumination-collection SNOM setup able to detect the 
photoluminescence signal emitted in the localized PCN modes. F corresponds to a 
Faraday optical isolator, T is a telescope formed by two lens, HM an hot dichroic 
mirror, P a long-wavelength pass filter, L a lens that focuses the PL signal to a spec- 
trometer (SP) and DTC a detector composed by an array of InGaAs photodiodes. The 
piezoelectric bimorph and scanning system are both driven by the same electronic 
central unit. 
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Figure 3.3: a) Four subsequent fabrication steps (1-4) performed to create the GaAs 
(light-gray) photonic crystal slab with InAs QDs embedded (red layer). The dark-gray 
region represents the Alo.7Gao3As sacrificial layer and the black layer is the SiO2 
mask. b) Typical spectrum of the quantum dots PL detected by SNOM equipped with 
dielectric probe in a region outside the photonic crystal, as highlighted by the inset. 
The excitation light is provided by the A = 780 nm cw diode laser with 60 uW power 
delivered to the probe. 
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Figure 3.4: a) SEM image of the investigated D2 PCN with lattice constant a = 
331 nm. b) Typical spectrum averaged over the D2 cavity detected by SNOM die- 
lectric probe as highlighted by the inset. The excitation light is provided by the A = 
780 nm cw diode laser with 60 uW power delivered to the probe. Both resonant 
modes M1 and M2 exhibit clear PL peaks. 
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Spectral Shift map 

Figure 3.5: Schematic explanation of the spectral shift evaluation technique. Two 
spectra, blue and red lines, are reported as a function of the wavelength and corre- 
spond to the signal collected by the SNOM dielectric tip in position A and B on the 
D2 nanocavity, respectively. The spectrum acquired inside the PCN (position B) has 
higher counts and a higher wavelength peak position with respect to the spectrum 
collected outside the cavity. By analyzing the number of counts collected as a function 
of tip position, we retrieve the intensity map of the mode. On the other hand, by re- 
porting the peak wavelength as a function of tip position we obtain the spectral shift 
map. 
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Figure 3.6: a)-b) Map of the intensity collected in the xy plane few nanometers above 
the slab at M1 and M2 peak wavelength, respectively. c)-d) Map of the induced spec- 
tral shift amplitude for M1 and M2, respectively. e)-f) FDTD electric field intensity 
distribution calculated at M1 and M2 unperturbed resonant wavelength, respectively. 
The scale bar in all the maps is 300 nm. 
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Figure 3.7: a) Map of the electric LDOS intensity calculated by FDTD at MI resonant 
wavelength. The dotted orange line highlights the position where the profiles reported 
in b) and c) are considered, both in the experimental and in the theoretical map. 
b) Comparison between the normalized vertical profiles of M1 obtained by the FDTD 
calculation (black line), by the PL intensity (red dots) and by the spectral shift analysis 
(blue triangles), respectively. c) Comparison between the normalized vertical profiles 
of M1 obtained by a convolution of the FDTD map with 70 nm FWHM two-dimen- 
sional Gaussian point-spread function (red line) and the experimental spectral shift 
profile (blue triangles), where the error bars are provided by the Lorentzian fit. 
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Figure 3.8: a) Schematics of the metal coated probe interaction with the magnetic field 
of light (blue arrows) localized in the PCN. b) The near-field probe is represented as 
a conducting ring that generates a magnetic dipole (green arrow) opposite to the out- 
of-plane magnetic field. c)-d) Map of the intensity collected in the xy plane few na- 
nometers above the slab for the M1 and M2 peak wavelength, respectively. e)-f) Map 
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of the induced spectral shift amplitude towards lower wavelength for M1 and M2, 
respectively. The excitation light is provided by the A = 780 nm cw diode laser with 
2.0 mW power delivered to the probe. g)-h) FDTD intensity distribution of the mag- 
netic field out-of-plane component calculated for M1 and M2, respectively. The scale 
bars in all the maps are 300 nm. Adapted from [10]. 


Q 


di 


Figure 3.9: a) Schematics of the campanile tip geometry composed of a gold tapered 
metal-insulator-metal waveguide fabricated at the end of a tapered glass fiber with the 
identification of the characteristic parameters: d=50 nm, 0=30°, Wo=20 nm, 
L=200 nm, Wi=200 nm, Do=10 nm, D=200 nm. b)-c) SEM images of the campanile 
tip with a 40 nm nanogap, which resembles the shape of a bell tower, as shown in the 
photo of the Berkeley campanile reported in the inset of b). 
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Figure 3.10: a) Finite element method calculation of the electric field enhancement 
distribution at the campanile tip apex for a nanogap size of 10 nm and A = 666 nm. 
The mode propagating in the fiber is excited by a plane wave linearly polarized 
along y. b) Same result of a), but restricted to the nanogap region and in a saturated 
color-scale contrast to show the nearly background-free enhancement at the tip apex. 
c) Electric field enhancement distribution at the apex of a conventional gold coated 
probe for an aperture size of 10 nm. d) FEM electric field enhancement calculation as 
a function of the incident wavelength for a campanile tip with a 2 nm gap (red curve) 
that extends over larger bandwidth than does the response of a bow-tie antenna (blue 
curve) with the same gap or a sharp gold tip with a 20 nm size aperture (gray curve). 
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Figure 3.11: a)-c) SNOM PL measurements performed by campanile tip on a InP 
quantum-wire, using He-Ne cw laser emitting at A = 633 nm and employing an ex- 
citation power of 100 uW inside the tip. a) PL intensity distribution evaluated at 2 = 
839 nm; b) topography image; c) waterfall plot of normalized spectra acquired at 
positions 1 to 11 reported in a) shows strong local spectral variations. d) and f) PL 
measurements performed on the same quantum-wire by a confocal setup. d) and f) are 
the confocal analogous of a) and c), respectively. e) SEM image of the wire. 
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Figure 3.12: a) SEM image of a different quantum wire with respect to Fig. (3.11). b)- 
c) PL intensity distribution evaluated at A = 802 nm by campanile probe and by con- 
focal setup, respectively. d) Comparison between intensity acquired in a)-c) along the 
central profile orthogonal to the wire length (red lines). e) Intensity profile (red curve) 
and topography profile (dotted curve) collected by SNOM along the wire length reveal 
a spatial resolution of about 40 nm. 
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Figure 3.13: a) Three-dimensional schematics of the experimental setup. The excita- 
tion light is provided by the cw diode laser (780 nm) with 120 uW power. DM is the 
dichroic mirror and DTC represents the detection unit. b) Typical PL spectrum col- 
lected by the campanile tip in the middle of the PCN with lattice constant a = 311 nm 
reported in the inset. c)-d) FEM calculations of the electric and magnetic field en- 
hancement distribution, respectively, at the campanile tip apex for a nanogap size of 
40 nm and 7 = 1.3 um. The mode propagating in the fiber is excited by a plane wave 
with electric field linearly polarized along y, it follows that the electric (magnetic) 
field at the tip apex is mainly y (x) polarized. e) FEM electric field enhancement 
calculated over the near-infrared spectral range for a campanile tip with a 40 nm gap. 
The dashed vertical lines highlight the experimental spectral region. 
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Figure 3.14: a) Schematic of the investigated PCN where the three yellow-blue rings 
A, B and C represent the positions where the spectra reported in b) are collected by 
the campanile tip. b) Normalized PL signal of the M1 mode collected in the three 
positions reported in a). In the same scan the evidence of both a red (A) and a blue 
(B) tip induced spectral shift, with respect to the peak detected outside the cavity (C), 
is clear. c)-d) Experimental spectral shift amplitude distributions for the M1 and M2 
mode, respectively. They are obtained reporting the maximum peak wavelength, eval- 
uated by a Lorentzian fit, as a function of the tip position. The white color, for each 
map, corresponds to the smaller-perturbed wavelength detected in the position about 
1 um away from the center of the PCN. 
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Figure 3.15: Finite-Difference Time-Domain calculations. a)-b) (E*)distribution nor- 
malized to the maximum for M1 and M2 mode, respectively. c)-d) (~H?) distribution 
normalized to the maximum for M1 and M2 mode, respectively. e)- f) Distributions 
evaluated using Eq. (3.3) with p (M1) = 0.47 and B(M2) = 0.55, respectively. These 
maps nicely reproduce the experimental spectral shift maps reported in Fig. (3.14 c- 
d). The spatial scale bar is 300 nm in every map. 
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Figure 3.16: Spatial resolution of the campanile tip perturbation imaging method. a)- 
b) Comparison between the spectral shift data and the FDTD calculation using 
Eq. (3.3) with (M2) = 0.55. The horizontal dashed lines indicate where the spatial 
profiles reported in c) are evaluated. The scale bars are 300 nm. c) Spectral shift data 
(red dots) collected, as a function of position, along the dashed line in a), and com- 
pared to the theoretical profiles obtained by the convolution between the FDTD map 
of b) and two-dimensional Gaussian point spread functions, characterized by FWHM 
of 40 nm (black curve) and 80 nm (green curve), whose intensity are normalized to 
the experimental spectral shift amplitude range. The error bars are provided by the 
Lorentzian fit. 
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Chapter 4 


Phase sensitive Fano imaging 


Photoluminescence based imaging methods require optically active samples with 
constrain of spectral and spatial matching between the photonic mode and the light 
sources, which may suffer of bleaching or blinking. Therefore, a pure optical method 
that can be applied on any kind of high Q nano-resonators to retrieve the confined 
mode distribution is actually missing. Here, for the first time, a resonant scattering 
setup is implemented on a scanning near-field optical microscope. 

This scheme allows us to perform hyperspectral imaging of the light states local- 
ized in optical nano-resonators, independently on the material they are made of. This 
novel imaging technique deeply exploits the nature of Fano resonances, which emerge 
from the interference between a resonant and a non-resonant signal. Actually, the de- 
tailed analysis of the Fano spectral lineshapes not only allows us to retrieve the elec- 
tric-LDOS of nanoscale modes, but also provides unprecedented insight into the con- 
fined electric field phase spatial modulation. In fact, Fano imaging intrinsically in- 
cludes a near-field integrated interferometer for a direct phase retrieval of any inves- 
tigated polarization components. Hence, Fano imaging constitutes a substantial ad- 
vance for mapping localized optical states. In particular, we prove that it stimulates a 
deep sub-wavelength imaging of both the electric LDOS and the electric field phase 
of the modes localized in photonic crystal nanocavities. 


4.1 Fano resonances 


In this section we present a brief review of the previously developed resonant 
scattering experiments on PCNs and the theory that accounts for the observed asym- 
metric spectral resonances, that significantly differ from the symmetric well-known 
Lorentzian shape of a conventional resonance. 

The beginning of the story traces back in 1935, when unexpected experimental 
evidences of strong asymmetric profiles in absorption spectra were observed in Ry- 
dberg atoms by Beutler [1]. In the same year U. Fano, Fig. (4.1), working with 
E. Fermi proposed the theoretical explanation that predicts the shape of the observed 
spectral lines, based on a superposition principle from quantum mechanics calcula- 
tion [2]. According to Fano: “The Beutler spectra showed unusual intensity profiles 
which struck me as reflecting interference between alternative mechanisms of excita- 
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tion”. The interpretation provided by Fano is based on the interference between a dis- 
crete excited state of an atom and states described by a continuum energy band, if the 
continuum bandwidth covers the energy of the discrete state. 

In 1961 he published a much more elaborated theoretical work, today considered 
his masterpiece, where are presented and analyzed peculiar shapes of the inelastic 
scattering spectra of electrons on He atoms [3]. The phenomenon described by Fano 
turned out to be universal, because the manifestation of interference effect does not 
depend on the peculiar system considered and it was called “Fano resonance”. There- 
fore, it is not surprising that, over the years, Fano resonances have been observed 
across many different branches of physics, including nuclear, atomic, molecular and 
condensed matter [4-11]. Recently, also a classical analogy oftwo coupled oscillators 
has been proposed [12]. 


Fano resonance takes place when a narrow band of virtually any origin (corre- 
sponding to a discrete atomic level) is superimposed on a spectrally flat background 
spectrum (corresponding to a continuous band of atomic states). In particular, they 
arises when a particle may reach the same final state via two different scattering path- 
ways, as highlight in Fig. (4.2 a): (i) scattering mediated by the interaction with a 
discrete energy state, where the field phase changes by z as a function of the wave- 
length (resonant path), or (ii) through direct scattering with a continuum of energy 
states where the field phase and amplitude are nearly constant (non-resonant path). 
Therefore, the analysis of Fano resonances provides fundamental information on the 
geometric configuration or internal potential of the system. 

Fano derived a simple expression for the cross section of an e.m. transition de- 
scribed by a given quantum operator (for instance the absorption of a photon) from an 
initial bound state (we shall consider the ground state) to a final continuum state owing 
the same energy of an excited discrete level, and that can be reached through the path- 


ways (i) and (ii) [3]: 


gee (4.1) 


1+? 


where the reduced energy e = 2(E — E,)/T is a function of the particle energy E, 
the energy and the linewidth of the discrete state Ey and T, respectively, while the 
phenomenological shape-parameter q is introduced by Fano so that Tq? /2 represents 
the ratio between the transition probability to the mixed state (discrete state modified 
by the interference with the continuum) and to the unperturbed continuum states in an 
energy bandwidth T [3]. 

Eq. (4.1) nicely reproduces the asymmetric observed lineshapes, as schematically 
reported in Fig. (4.2 b). Moreover, in order to highlight the basic properties of the 
Fano relationship, in Fig. (4.3) are plotted various lineshapes as a function of the re- 
duced energy for different values of the parameter q. The shape parameter governs 
the lineshape asymmetry. In fact, as reported in Fig. (4.3 a), the case of q = 0 de- 
scribes a symmetrical dip, which corresponds to the case where the non-resonant scat- 
tering path (ii) is the dominating contribution. When q = 1 the continuum and discrete 
transitions are of the same strength, resulting in a perfectly odd-profile with respect 
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to the energy axis. For increasing q the lineshape is mainly determined by the transi- 
tion through the discrete state and, in the limit q — +09, it approaches the standard 
Lorentzian profile of a Breit-Wigner resonance. For any q # 0 the Fano profile shows 
the maximum amplitude value at energy Ey +T/(2q) and the minimum at Eo — 
qI'/2. The resonant energy of the discrete level may lie somewhere between the max- 
imum and the minimum of the asymmetric profile, and the parameter q defines the 
relative deviation. For q — +00, the resonant energy coincides with the lineshape 
maximum, while in the case of q = 0 it coincides with the minimum. Notably, for 
q = 1 it is located exactly at half the distance between the minimum and maximum. 
The effect of the shape parameter sign inversion is evaluated in Fig. (4.3 b), where 
two profiles calculated with opposite q values are reported, showing that for a nega- 
tive shape parameters the lineshape undergoes a reflection with respect to the vertical 
axis. Furthermore, every lineshape shows an energy for which the amplitude vanishes, 
implying that for that energy destructive interference occurs. The basic properties that 
distinguishes Fano resonance from conventional resonances, in which an external ex- 
citation at a proper energy enhances the system response, is therefore the presence of 
a possible destructive interference that manifest itself as an antiresonance, that means 
o = 0 for a particular energy. In Fano resonances the character of the interference 
changes rapidly as a function of the energy, giving rise to asymmetric lineshapes. 


4.1.1 Fano resonances in photonic crystal cavities 


In optics, Fano resonances have been widely observed especially in plasmonic 
based nanostructures, metamaterials and photonic crystals; large reviews can be found 
in [11, 13]. Since we will focus on the Fano lineshape variations in scattering near- 
field measurements performed on PCNs, it is important to introduce previous experi- 
ments and theoretical works that deal with Fano resonance in photonic nanocavities. 
In fact, it has already been demonstrated that the q parameter can be tuned. For in- 
stance, the absorption spectra of a single quantum dot vary from a symmetrical profile 
to an asymmetrical one with increasing laser power, indicating an enhancement of the 
continuum transition [14]. Even reflectance spectra on single quantum-well structure 
performed at low temperature show periodic Fano lineshape variation, from disper- 
sive-like to absorptive-like and inverse absorptive-like, as a function of the barrier 
thickness layer [15]. This behavior is due to the interference between the reflected 
waves from the sample surface and from the quantum-well, and their analysis carries 
information about both the absorption coefficient and the exciton transition. 

A wealth class of studies concerning Fano resonances have been performed on 
photonic crystals nanocavity platform. In a disordered photonic crystal Fano reso- 
nance may occurs via the interaction of the narrow Bragg band with the extra scatter- 
ing induced by disorder [16, 17]. In photonic crystals slab structures, a three-dimen- 
sional frequency and time-domain theoretical analysis proved the presence of spectral 
sharp asymmetric lineshapes superimposed upon a smooth background [18]. The 
sharp features come from the guided resonances confined in the slab that can be cou- 
pled to radiation modes, while the slowly-varying background corresponds to Fabry- 
Pérot spectral oscillations of a uniform slab (without holes). The Fano lineshapes 
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emerges as the interference between a non-resonant pathway, which corresponds to 
the direct transmission or reflection, and a resonant pathway, in which the remaining 
portion of the incident light interacts with the exponential decaying amplitudes of the 
slab guided resonances. Notably, the resulting spectral sharp profiles show different 
shapes and symmetry, depending on the relative spectral position (phase shift) be- 
tween the guided mode and the background. 

For analytically investigating the Fano resonances in optical cavities a temporal 
coupled-mode theory has been developed [20]. This theory models the spectral re- 
sponse in term of the transmission and reflection coefficients of the medium together 
with the resonant energy and the quality factor of a localized mode. It has been applied 
to the interesting case of an optical system, consisting of a waveguide side coupled to 
a single-mode cavity, that offers potential applications in optical modulations and 
switching [19, 21]. In particular, concerning ref. [19], as reported in Fig. (4.4), they 
studied a system consisting of a square based photonic crystal of high-index dielectric 
rods, in which a line defect forms the waveguide that is coupled to a point defect 
nanocavity and in the middle near the nanocavity presents two small dielectric cylin- 
ders separated by l, to provide partial reflection for waveguide modes. The FDTD 
calculations of Fig. (4.4 a-c) prove how the presence of partially reflecting elements 
induces the spectral lineshape evaluated at the output of the waveguide to be asym- 
metric, with a fast frequency switch from zero to maximum output intensity. Differ- 
ently, by removing the two rods inside the waveguide the spectrum is found to be an 
inverted Lorentzian dip. However, both the lineshapes reported in Fig. (4.4 c) can be 
described by Eq. (4.1), ad discussed in Fig. (4.3 a). In fact, at the waveguide output 
we find the interference between the field directly transmitted through the line defect 
and the field resonant to the cavity mode, that is trapped inside the cavity for a char- 
acteristic time (inversely proportional to cavity Q) and then after a delay is coupled 
back to the waveguide mode. The role of the reflecting elements lies only in introduc- 
ing a given retard, that is a phase shift, between the two pathways, and then changing 
the shape of the output spectrum. The coupled-mode theory calculations, reported in 
Fig. (4.4 d-g), faithfully reproduce the FDTD results and clearly show that, by arbi- 
trary varying the phase shift between the cavity mode and the background, the output 
spectrum shape drastically changes. Therefore, the method allows to catch the funda- 
mental behavior of Fano resonances. 

Profound insight in high-Q PCNs has been pursued by resonant scattering far- 
field experiments which display Fano-like resonances [22-25]. These alternative tech- 
niques exploit the light resonant scattering using a large bandwidth laser or a tunable 
laser able to span the spectral range of the investigated nanocavity illuminated by 
standard microscope objective. In this geometry Fano resonances arise from the inter- 
ference of (i) light resonantly coupled to the cavity mode and then released in free 
space, and (ii) light directly diffused by the surface of the photonic crystal pattern, 
through a continuum of extended states, that is not coupled to the cavity mode. How- 
ever, it is difficult to observe the cavity spectrum directly, because only a small frac- 
tion of the incident light couples to the PCN owing to poor mode matching between 
the gaussian probe beam and the cavity mode. For this reason, the signal is monitored 
in crossed-polarization detection scheme, thus eliminating the huge direct reflection 
signal and allowing the detection of the scattering amplitude. In particular, if the PCN 
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mode is expected to be linearly polarized along a given in-plane axis, in order to max- 
imize the detection of the resonant mode over the background, the incoming light is 
set to be linearly polarized at 45° with respect to the resonant mode axis. Then, a 
second linear polarizer, whose axis is orthogonal to the polarization of the impinging 
light, filters the backward scattered light. This configuration, combined with the cross- 
polarization detection, suppresses the non-resonant contribution, as desired, while re- 
ducing the resonant contribution only by a factor of 2. It turns out that the two contri- 
butions are almost balanced and the Fano lineshapes can be nicely observed. The ob- 
tained Fano spectral lineshape intensity F(E) are commonly reproduced by the fol- 
lowing relationship derived from Eq. (4.1) [23]: 


[q+2(E-E0)/T]? 


F(E) = Ao + Po T Tar 


(4.2) 


where Ao and Fo are amplitude factors. The former corresponds the minimum 
intensity (Fmin = Ao); the latter to is related to the amplitude of the Fano lineshape, 
defined as the difference between the maximum and the minimum intensity, that cor- 
responds to Enay — Fmin = Fo(1 + q°). 

By fitting the experimental data with Eq. (4.2) a characterization of the cavity 
mode in non-optically active samples is then achieved, concerning both the resonant 
wavelength (Eo) and the quality factor (Q = E,)/T). But a resonant scattering tech- 
nique that, by means of Fano resonances detection in the near-field, would map the 
localized field distribution with a sub-wavelength resolution, is actually missing. In 
ref. [23] are also found clear lineshape variations, from strongly asymmetric profiles 
to symmetric Lorentzians, as the impinging beam waist is increased over the whole 
cavity area. This demonstrates an enhancement of the scattering in the continuum 
pathway and a decrease of the asymmetry parameter q. 

Moreover, Eq. (4.2) can be understood more in details. In fact, considering ener- 
gies far away from the cavity resonance, |E/E | > 1, the amplitude F tends to the 
value B = (Ao + Fo), which has to be equal to the background scattering signal non- 
resonant the cavity mode. Therefore, Eq. (4.2) can be written as [25]: 


_ (ata) _ 
FE) =B + Fo eee aE 1} (ta) 


Notably, Eq. (4.3) remains invariant under the transformations (q’ = —1/q; Fo = 
—F)q7). Hence, two Fy with opposite signs give the same lineshape, and without any 
loss of generality is commonly considered that Fọ>0 [25]. 


4.2 Deep sub-wavelength Fano imaging 
The resonant scattering technique overcomes the drawback of the most common 


experimental approaches to achieve PCN mode imaging characterization. They 
mainly rely on either internal light sources (PL experiments), or on evanescent cou- 
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pling between the investigated PCN and an optical waveguide [26, 27]. The first ap- 
proach, which is applicable only to the small class of optically active media, such as 
the PCNs analyzed in Chapter (3), cannot be extended to silicon, glass or metal based 
materials. It also requires a spectral and spatial matching between the light source and 
the photonic modes. Moreover, it is limited by pump-induced losses of the active me- 
dium and to the presence of efficient light emitters, which may lower the cavity qual- 
ity factor. Although the second approach is well suited for optically passive PCN, it 
requires the fabrication of an additional coupling waveguide, which complicate the 
experimental geometry and also may induce a loading effect on the PCN modes that 
can result drastically perturbed. 

An alternative approach, based on cathodoluminescence, provides a deep sub- 
wavelength (about 4/20) imaging of a two-dimensional SizN4 based PCN in a broad 
spectral range [28]. However, this technique is limited to the collection of the electric 
LDOS in dielectric regions and it is prevalently sensitive to the out-of-plane electric 
field component, which is the less important one in two-dimensional photonic crystal 
nanocavity on slab. 

Pure photonic investigation techniques are potentially more powerful due to the 
covering of a broad spectral range, the possibility to address the distributions of local- 
ized modes also in air regions or in metal devices, the absence of detection polariza- 
tion constraints and, in combination with heterodyne detection in an external interfer- 
ometer, the possibility of phase retrieval as they preserve coherence [29, 30]. This 
approach relies on near-field surface enhanced light scattering performed with aper- 
turless probes (a-SNOM), and have recently been applied to various contexts [29, 30]. 
The great sensitivity of the method is largely due to the use of high refractive index 
material (metal or silicon) probes in order to enhance the scattering. But, in the case 
of optical micro and nano resonators (such as ring resonators, microdisk cavities, 
PCNs) investigation, the a-SNOM probe would perturb too much the local dielectric 
environment, with respect to glass tapered probes, thus deteriorating the high quality 
factor and also detuning/mixing modes localized in coupled PCNs or in random cav- 
ities [31, 32]. Moreover this technique, like cathodoluminescence, suffers from being 
mainly sensitive to the out-of-plane component of the electric LDOS, and it is not 
suited for TE-like localized modes. 

Therefore, it is clear that pure optical imaging techniques based on resonant scat- 
tering, that can be applied on high Q resonators in order to retrieve the intensity and 
the phase of the confined electric field, is highly desirable. In order to meet this crucial 
need for nanophotonics, we modify the scanning near-field optical microscopy setup 
presented in Fig. (3.2) to perform resonant scattering spectroscopy that, equipped with 
dielectric tapered probes, enables the hyperspectral imaging of the electric LDOS with 
outstanding spatial resolution, also for non optically active samples. We investigate 
the D2 photonic crystal nanocavity, exploiting a detailed spectral analysis of the ob- 
served Fano resonances. In addition, we use aperture aluminum coated probes and 
campanile tips to prove that the method intrinsically includes a near-field integrated 
interferometer for phase retrieval without the need of external heterodyne detection, 
thus revealing the sub-wavelength spatial modulation of both the intensity and the 
phase of the localized electric fields. 
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4.2.1 Resonant back-scattering 


The schematics of the SNOM setup in the illumination/collection geometry 
equipped with a dielectric probe to perform resonant back-scattering (RBS) measure- 
ments is shown in Fig. (4.5 a). In particular, to cover the interesting spectral range in 
a single scan only we use a pulsed supercontinuum (SC) laser (Leukos-STM with 
pulse temporal width of 1 ns, repetition rate of 5 KHz), whose emission at the output 
of a photonic crystal fiber is reported in Fig. (4.5 b). Light coming from the laser is 
filtered by a long-wavelength pass filter (LP) that allows only 2 > 1.2 um to be trans- 
mitted, and it is linearly polarized by the transmission through a polarizing beam split- 
ter cube (PBS), whose extinction ratio is 2 - 104. Then it is coupled to the dielectric 
near-field probe to illuminate the sample in the near-field. We convey inside the near- 
field probe an almost constant power spectral density of about 5 uW per 10 nm. 

Crossed polarization detection is performed to overcome the huge reflection sig- 
nal, which has the same polarization of the incident light. In fact, only balancing the 
resonant and the diffused signals we are able to observe clear Fano lineshapes. The 
implementation of SNOM-RBS therefore requires the polarization control in both in- 
put and detection channel paths, which share the near-field probe. Therefore, a Babi- 
net- Soleil compensator, reported in Fig. (4.5 c-d), is mounted on the optical fiber in 
order to control the light polarization at the end of the tip by twisting and squeezing 
the fiber. This provides a linear polarization, with an extinction ratio of 120 in the 
(1320 + 20) nm wavelength range, for the light incident on the sample with the elec- 
tric field oriented at 45° with respect to the x axis of the photonic crystal. The polari- 
zation control is performed by collecting the light transmitted through the sample by 
a 50x 0.4 NA objective (Ob.) and by measuring with a photodiode detector (PD, 
Femto-Watt series by New Focus), in a 40 nm spectral band, the ratio between the 
light polarized as the light impinging the sample, and the light orthogonally polarized 
with respect to the incoming beam. Moreover, we test that the investigated sample is 
not a birefringent material by measuring the polarization maintaining properties of 
linearly polarized light at 1300 nm. In fact, the incoming collimated free-standing 
laser light shows the same polarization degree after the transmission through the sam- 
ple. The backward scattered light collected by the probe is filtered in crossed polari- 
zation configuration by the polarizing beam-splitter cube, and spectrally analyzed by 
the same detection scheme employed in PL measurements, see Fig. (3.2), with typical 
integration time of 100 ms, leading to a spectral resolution of 0.11 nm, that in photon 
energy units corresponds to 0.08 meV. Spectra are collected at every tip position on 
the sample surface. A typical raw RBS spectrum collected in the middle of D2 
nanocavity with lattice constant 321 nm, is reported in Fig. (4.6 a). The spectrum dis- 
plays two asymmetric resonances corresponding to the two lower energy modes (M1 
and M2) superimposed to a large wavelength dependent background, whose intensity 
is generally comparable to the maximum observed amplitude (i.e. difference in counts 
between the maximum and minimum of a single sharp resonance). This background 
consists only in a non-resonant scattering signal that can be experimentally addressed 
by the averaged signal detected in a region located about 4 um away from the center 
of the cavity, where the photonic crystal is periodic and the PCN localized modes are 
no more detectable, as shown in Fig. (4.6 b). Of course, in order to observe clear Fano 
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lineshapes the resonant and non-resonant scattering signals must be almost balanced. 
This is the scope of the cross-polarization setup where the non-resonant contribution, 
which is mainly polarized as the incoming beam and larger than resonant channel, is 
strongly suppressed. Nevertheless, after the cross-polarization procedure, a non reso- 
nant contribution is still left and it is used to produce interference with the resonant 
channel. This remaining non-resonant contribution is polarized as the resonant mode, 
therefore for mode MI is linearly polarized along the x axis of the PCN. 

In order to directly observe the spectral region where the negative or positive in- 
terference between the two scattering signals occurs, the background has to be sub- 
tracted, as we reported in Fig. (4.6 c). In the resulting spectrum both resonances show 
a dispersive-like lineshape, fingerprint ofthe Fano resonance, and each mode includes 
both positive and negative values, evidence of a constructive and destructive interfer- 
ence, respectively. When the background-subtracted spectrum is averaged over adja- 
cent positions during the SNOM scan, we achieve high signal/noise ratio lineshapes, 
see Fig. (4.6 d). Only after the background subtraction of each acquired spectrum we 
perform the Fano function fit of both modes intensity, as a function of the photon 
energy, by employing Eq. (4.2). The experimental data are well reproduced by the 
Fano formula, as shown in Fig. (4.7 a-b), and an accurate evaluation of the relevant 
parameters is pursued. 

The power of the near-field microscope, configured for resonant scattering spec- 
troscopy, fully emerges when optical imaging is achieved. By mapping, as a function 
of position, the parameters deriving from the Fano lineshape analysis, it is possible to 
retrieve the electric-LDOS imaging of both localized modes. In order to map the RBS 
amplitude for M1 and M2 we plot in Fig. (4.8 b and e), respectively, the value 
Fo(1 + q?) provided by the fit, that corresponds to the difference in counts between 
the maximum and the minimum of each Fano lineshape. These maps reproduce the 
main features of the electric LDOS at the energy of the resonant modes obtained by 
finite difference time domain (FDTD) calculations reported in Fig. (4.8 d and g), re- 
spectively. However, they show a limited spatial resolution (about 250 nm) due to the 
fact that the amplitude map is directly proportional to the light collected by the die- 
lectric uncoated near-field probe, whose resolution is bound to the few hundred na- 
nometers large collecting area. 

An improved spatial resolution is achievable from the same set of data making 
the best use ofthe tip local perturbation on the PCN dielectric environment. The slight 
perturbation induced by the near-field tip gently shifts the resonant modes towards 
lower energies proportionally to the strength of the local unperturbed electric field 
intensity, as demonstrated in Chapter (3) [26, 34]. This phenomenon is underlined by 
the two normalized intensity spectra reported in Fig. (4.8 a). The energy spectrum 
reported with black dots shows the largest tip perturbation and it is collected in the 
proximity of the maximum of the M1 intensity, corresponding to the red circle in 
Fig. (4.8 c). On the other hand, if the tip is located in the green circle position in 
Fig. (4.8 c), the PCN experiences the slightest disturbance resulting in a higher energy 
resonance, as highlighted by the spectrum reported with withe triangles in Fig. (4.8 a). 

Figures (4.8 c and f) show the spatial distributions of the tip induced spectral shift 
for M1 and M2, respectively, as obtained by reporting the fitting parameter 
(Eo — Eo min) as a function of position for each mode. The results are in excellent 
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agreement with the FDTD based electric LDOS calculations, shown in Fig. (4.8 d and 
g). Since the smallest LDOS spatial feature is about 100 nm, in order to estimate the 
spatial resolution of the Fano imaging technique, we compare the experimental M1 
spectral shift map to the FDTD map convoluted with two-dimensional Gaussian point 
spread functions, characterized by different FWHM. In Fig. (4.9 a-b) are reported 
vertical and horizontal profiles, as highlighted in Fig. (4.8 d), of the experimental 
spectral shift data compared to theoretical profiles obtained by convolution with 
Gaussian of FWHM equal to 50 nm, 70 nm and 100 nm. The function with 70 nm 
FWHM minimizes the y? distribution. Therefore, the combination of SNOM and RBS 
spectroscopy allows us to obtain high fidelity electric LDOS mode imaging with an 
ultra sub-wavelength spatial resolution of 70 nm, down to 4/19, comparable to the 
PL technique presented in Chapter (3) and to the so far best result obtained on 
PCNs [28]. 

Overall, the spatial resolution achieved, as well as the spectral shift amplitude 
obtained through Fano fits for both resonant modes are comparable to the spatial res- 
olution and amplitude evaluated by Lorentzian fits of PL spectra measured with sim- 
ilar dielectric probes on equivalent D2 nanocavities. This feature confirms that the 
dielectric perturbation technique does not depend on the employed imaging method 
but it is a property of the cavity-probe system only. It is straightforward to notice that, 
although the investigated PCN has light source embedded (InAs quantum dots), in 
resonant scattering experiments we do not excite the QDs, since the illumination is 
performed with wavelengths grater than 1.2 um and at very low power. In fact, outside 
the photonic crystal, where only QDs are present in the slab, we do not collect any 
signal resembling the QDs emission, see Fig. (3.3 b). 

The presented Fano imaging can be potentially applied to any nano-photonic and 
nanopasmonic optical resonator that exhibits sharp spectral resonances, by confining 
light, for example, in silicon based materials as well as glass, polymer or metallic 
materials. By mapping the Fano lineshapes as a function of position in the near-field 
of the D2 nanocavity not only spectral shift and amplitude variations are evident. In 
fact, if for sake of simplicity we focus on the behavior of the M1 mode, which is 
known to possess a liner polarization component parallel to the x axis of the PCN, we 
find strong lineshape modifications during a single large scan over the PCN, as re- 
ported in the two spectra of Fig. (4.10 a). The Fano lineshape collected in position A 
exhibits positive signal with respect to the background, while the spectrum detected 
in B has almost exclusively negative values, as if the interference abruptly reverses 
from constructive to destructive when moving from point A to point B. 

Moreover, in resonant scattering SNOM experiments a clear signature of the lo- 
calized modes is detectable over a larger area with respect to PL measurements, until 
3 um away from the cavity center, and also in the regions where LDOS is negligible. 
But this is not the only advantage, since the Fano lineshape analysis allows for map- 
ping resonance amplitude, spectral shift, quality factor (as in PL measurements) and 
also, more strikingly, shape modification, which is a unique property of the resonant 
scattering technique. The idea is that the additional information contained in Fano 
shape modifications as a function of position allows to retrieve the spatial phase mod- 
ulation of the localized modes. In fact, we already discussed that the shape of the Fano 
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resonance is given by the interference between non-resonant light diffused with a con- 
stant phase over the whole PCN and resonant signal, which has a defined phase spatial 
dependence. In our imaging setup, by employing a coherent light source and looking 
at the changes of the Fano lineshape as a function of position, we should get an equiv- 
alent integrated interferometer to locally probe the relative phase distribution of the 
PCN localized modes. 

In order to obtain a distribution of the lineshape modifications it is a straightfor- 
ward result to map the RBS intensity distribution at the energy of the most intense M1 
peak, as reported in Fig. (4.10 b), which maximizes the counts difference. From this 
map the MI distribution is recognizable in its main features and lineshapes inversions, 
as the one reported in Fig. (4.10 a), only occur between the blue and red regions of 
Fig. (4.10 b). But these regions are 3 um apart and do not correspond to adjacent 
intensity lobes. So, they do not correspond to the actual mode phase distribution be- 
havior, as it is predicted by FDTD calculation of the real part of the MI electric field 
x component evaluated in a plane placed at 22 nm from the slab surface, reported in 
Fig. (4.10 d). But if the intensity distribution is evaluated at an energy that maximizes 
the color contrast between regions with opposite intensity sign, we obtain a nice agree- 
ment with the calculated real part of the electric field x polarized evaluated at 200 nm 
height and shown in Fig. (4.10 e). The dielectric tapered probe collects light also from 
the lateral sidewall at longer distance from the PCN than the probe apex. This mixing 
of near-field and intermediate-field detection fairly gives a probe effective collection 
placed at about 200 nm distance from the sample surface, causing a depression of the 
obtained spatial resolution [24]. 

For an accurate sub-wavelength phase-sensitive imaging of the electric field com- 
ponent of light localized in the PCN a probe able to collect only the near-field contri- 
bution is therefore recommended. However, the RBS setup does not allow the use of 
Al-coated probes because of the extremely low throughput. On the other hand, the 
campanile near-field tips presented in Sec. (3.3), by exploiting the surface plasmon 
assisted enhancement, show high throughout and can be employed in the illumina- 
tion/collection geometry. In Fig. (4.11 a) are reported three RBS background-sub- 
tracted spectra collected by campanile tip on different positions corresponding to ad- 
jacent electric LDOS lobes of the mode M1. A clear lineshape variation is present 
between positions C and D, while a similar spectrum of C is recovered in position E. 
In order to highlight these lineshape modifications as a function of tip position the 
intensity map evaluated at E = 931.28 meV is reported in Fig. (4.11 b). It shows a 
reversal of the interference condition between adjacent LDOS modes, that is in good 
agreement with the electric field phase variations reported in Fig. (4.10 d). Therefore, 
the local changes of the Fano lineshape give direct access to the spatial modulation of 
the localized mode phase by means of an integrated interferometer related to the na- 
ture of the Fano resonances. It must be stressed that here we are not interested to the 
actual lineshape, (that can be a peak, a dip ora a dispersive-live curve), but only care 
about the lineshape inversions in different positions, since the peculiar Fano lineshape 
depends also on the daily setup alignment that modifies the non-resonant signal and 
on the electric LDOS intensity. We estimate the spatial resolution of the phase sensi- 
tive Fano imaging by fitting the experimental data obtained by a vertical profile cut 
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of Fig. (4.11 b) with the error function, which comes out by convolving the experi- 
mental Gaussian point spread function and the supposed step-like variation of the 
phase modulation. We obtain the best fit using a Gaussian point spread function with 
FWHM equal to 85 nm, that corresponds to a deep sub-wavelength spatial resolution. 
The observed spatial resolution reduction is likely due to the prolonged use across 
many scan experiments, that induces the opening of the metallic nanogap. 

For completeness, by fitting every spectrum acquired with the campanile tip with 
Eq. (4.2), in Fig. (4.11 d and e) are also reported the spectral shift map and the Fano 
lineshape amplitude map. These images fairly well reproduce the M1 behavior when 
the PCN is observed and perturbed by the campanile tip, as discussed in Sec. (3.3.2). 
In fact, the simultaneous blue and red spectral shift are present in the regions of max- 
imum magnetic and electric field intensity, respectively, with a strength comparable 
to that observed in Fig. (3.14 d). 


4.2.2 Resonant forward-scattering 


In order to ensure a homogeneous excitation of the whole PCN with fixed relative 
phase of the incoming field between adjacent points we also implement a near-field 
collection setup in transmission geometry, whose schematics is reported in Fig. 
(4.12 a). This SNOM setup allows us to detect resonant forward scattering (RFS) of 
light from the investigated nanocavity. The linearly polarized light of the supercon- 
tinuum laser (transmitted through a Glan- Thompson polarizing prism, with extinction 
ratio of about 3 - 105, and oriented at 45° with respect to the PCN x axis) illuminates 
the sample from the far-field by means of a 50x objective (NA = 0.4) that ensures an 
excitation spot diameter of about 2 um. The forward scattered light collected by an 
Al-coated aperture near-field probe (produced by LovaLite) is then filtered in crossed- 
polarization configuration and analyzed by the same polarizing beam-splitter and de- 
tection scheme of the RBS setup. By means of the Babinet-Soleil compensator tool 
the crossed-polarization condition in checked (providing a linear polarization re- 
sponse of the system fiber plus tip with an extinction ratio of 120) by minimizing the 
signal transmitted by PBS and detected by PD when the incoming light shines outside 
the photonic crystal on the bulk GaAs. 

Figures (4.12 b) shows three RFS background-subtracted spectra spectra collected 
by the Al-coated aperture tip in three different cavity positions, corresponding the J, 
K and L defined in the inset of Fig. (4.12 a) and lie into adjacent M1 intensity lobes. 
All the reported spectra show both positive and negative values, depending on the 
interference condition between the resonant and non-resonant scattering channels. In 
particular, the spectra collected in J and L positions show a constructive interference 
for lower energies, while for higher energies the interference becomes destructive. On 
the contrary, the spectrum collected in K has an opposite lineshape, and this inversion 
is a clear indication that the resonant mode in K is characterized by a 7 phase variation 
with respect to J and L regions. Furthermore, the recovery of the J lineshape in the L 
position, points out that the photonic mode has the same phase in these regions. 

In order to confirm the ability to retrieve the mode phase modulation through Fano 
lineshape changes, we perform three-dimensional FDTD calculations, as reported in 
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Fig. (4.13 a). A linearly x-polarized plane wave impinges the bare nanocavity and the 
forward signal is collected by three detectors placed 10 nm above the sample surface 
at the positions indicated by a, £ and y. Figure (4.13 b) shows that the spectra calcu- 
lated in a and y are characterized by a similar lineshape, which differs only because 
the electric LDOS has distinct values in the two positions (i.e. the ratio between reso- 
nant and non-resonant scattering is slightly changed). The spectrum collected by the 
detector in f shows a lineshape inversion along the energy axis with respect to the 
other two calculated spectra, without any effective change of the resonant energy: 
Eo(a) = Eo(Y) = 946.639 meV, Ep(B) = 946.636 meV. The three investigated po- 
sitions are shown in the near-field calculation of the real part of the electric field of 
MI x component, reported in Fig. (4.13 c), to highlight that the localized field phase 
is reversed by m with respect to æ and y positions. These calculations provide the 
interference reversal evidence as a function of position. 

By reporting the RFS intensity map at a proper energy, as shown in Fig. (4.13 d), 
we can image the lineshape modifications also in resonant forward scattering. The 
spatial intensity sign alternation along the vertical direction nicely agrees with the 
phase modulation of the calculated electric field of Fig. (4.13 c). The good fidelity 
with which the map of RFS intensity reproduces the theoretical phase alternation map, 
proves the actual phase modulation imaging of the localized mode. As was done in 
the case of the campanile tip, we estimate the spatial resolution of the phase sensitive 
Fano imaging performed in RFS by a coated probe, by fitting the vertical profile of 
the experimental data in Fig. (4.13 d) with the error-function. Data and fitted curve 
are reported in Fig. (4.13 e), showing a mapping the phase modulations with a spatial 
resolution of 110 nm, as given by the FWHM of the Gaussian point spread function. 
The obtained resolution is limited by the probe aperture, which is nominally about 
200 nm. Therefore, we also demonstrate that the phase modulation retrieval is inde- 
pendent on the tip nature and therefore the possibly tip induced phase distortions do 
not dramatically alter the phase imaging. In fact, with the Campanile probe in RBS 
configuration we obtain qualitatively analogue spatial lineshape inversions, with re- 
spect to Al-coated probes, which result in a similar intensity map, even if the spatial 
resolution are different. Finally, to address the issue of reproducibility the same meas- 
urement are repeated more times, giving the same results for every probe but, demon- 
strating that a selection of good tips is needed, both dielectric and metallic ones. In 
particular, concerning the electric LDOS imaging performed by a selected tip, we ob- 
tain every time a similar spectral shift map and analogue value of maximum spectral 
shift, independently on the alignment or on the investigated nanocavity. Concerning 
the phase modulation imaging by the Fano lineshape analysis, every selected probe is 
able to collect the resonant scattering signal and to show clear Fano lineshapes. Alt- 
hough the same Fano lineshape is not detected again in the same position on the PCN, 
if the setup is re-aligned or if the near-field probe is replaced, the same mode lineshape 
inversions occur between the same spatial intensity lobes, thus giving equivalent 
phase imaging results. We believe that spurious reflections may play a role in chang- 
ing the amplitude and phase of the non-resonant scattering and that the method limi- 
tations are related to the pre-selection of good tips and to the aging of the tip after 
several measurements. 
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So far, we investigated the linearly polarized M1 mode. For this mode, to equili- 
brate the detection of resonant and non-resonant contributions, we took advantage of 
the crossed-polarization technique obtaining well resolved Fano resonance spectra. 
On the other hand, addressing the phase modulation imaging of the PCN excited mode 
M2 requires a more challenging experiment, since M2 shows an elliptical polariza- 
tion, where both x and y polarization components are present (with the y component 
larger than the x one) [35]. Obviously, the phase distributions of the two perpendicular 
components are quite different, as reported in the FDTD near-field calculations of Fig. 
(4.14 d and h), and we must develop an experimental method to separate them. 

Whenever the polarization of a given mode is unknown the cross-polarization de- 
tection apparatus can not be used. We then apply some modification to the RFS setup 
illuminating the sample alternatively with x or y linearly polarized light to retrieve 
separately the x and the y M2 phase spatial distributions. For each independent polar- 
ization we succeed in balancing the cavity resonant contribution and the transmitted 
light not coupled to the cavity, by exploiting the small size aperture of the Al-coated 
probe and eliminating the cross-polarization setup. In fact, differently from dielectric 
probes, metalized probes collect signal only from the small aperture placed few na- 
nometers on top of the sample, spatially filtering out very efficiently the non-resonant 
component, thus increasing the ratio between the resonant and the non-resonant chan- 
nels in detection. This polarization sensitive apparatus cannot be implemented in 
RBS, since releasing the crossed-polarization scheme by using dielectric probes leads 
to a huge background signal, which hides the resonant channel. On the other hand, by 
using metal coated probes it is not possible to collect any signal because of the very 
low tip throughput. The spectra reported in Fig. (4.14 a and e) show that the signal to 
noise ratio is still quite high and it allows us to image the phase modulation distribu- 
tion of M2 in both polarization components by mapping the mode intensity, as shown 
in Fig. (4.14 c and g), respectively. The image fidelity of the experimental phase mod- 
ulation is higher for the stronger polarization component (parallel to y axis) with re- 
spect to orthogonal polarization component (parallel to x axis), as inferred by the 
comparison between the experimental intensity maps in Fig. (4.14 c and g) and the 
real part of the x and y component of the electric field calculated by FDTD, shown in 
Fig. (4.14 d and h), respectively. Moreover, it is important to stress that the described 
experimental configuration for Fano phase modulation imaging does not need any 
pre-knowledge of the localized mode polarization. 

In conclusion, we achieve an ultra-subwavelength Fano imaging of the electric 
LDOS and of the phase modulation of localized modes in a two-dimensional PCN by 
using a pure optical effect that bridges together resonant light scattering and near-field 
spectroscopy. This method has the considerable advantage of being applicable to res- 
onators based on any kind of material and in a wide spectral range. In addition, for the 
first time, we directly probe the relative phase modulation of the localized mode with- 
out the use of an external interferometer, by exploiting the interferential nature of the 
Fano resonances. Near-field Fano imaging opens up new strategies to investigate with 
a deep-subwavelength spatial resolution the electric LDOS and the phase distribution 
of modes localized in a wealth class of nanophotonic and nanoplasmonic resonators, 
also based on materials without optical emission, and it can image different polariza- 
tion components. 
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Figure 4.1: Prof. U. Fano (1912-2001), Italian theoretical physicist pioneer in devel- 
oping the theory of resonant configurations interaction to describe asymmetric spec- 
tral lineshapes. 


=> 


i 
i 
® i 
1 
| 
© 
i 
i 
i 


Autoionization\ 


2 2 2 
Ig) d 1 poe qe &. E = (i q) 
€e +1 e +1 e +1 


Figure 4.2: a) Fano resonance as a quantum interference between (i) autoionization of 
two excited electrons followed by the Auger effect and (ii) direct ionization of a deep 
inner-shell electron. The two pathways are represented as transitions from the atomic 
ground state |g) either to a discrete excited autoionizing state |d) or to a continuum 
of states |c), where the discrete energy of |d) lies in the continuous range of values of 
|c). Dashed lines indicate double excitations and ionization potentials. b) Illustration 
of the Fano formula of Eq. (4.1) as a superposition of the Lorentzian line shape of the 
discrete level with a flat continuous background. Figures adapted from [11]. 
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Figure 4.3: a) Fano lineshape as a function of the reduced energy calculated from 
Eq. (4.1) in arbitrary units for different values of the shape parameter: q = 0 blue 
curve, q = 1 red curve, q = 2 green curve. b) Fano lineshapes with opposite values: 
of q = 1 red curve, q = —1 black curve 
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Figure 4.4: The system studied in [19] consists of a square lattice of dielectric rods 
with a line defect, in which two small cylinders separated by l, highlighted by yellow 
dots, are introduced to provide partial reflection. A nanocavity is side-coupled to the 
waveguide. a)-b) FDTD calculated electric field steady state distributions relative to 
the frequency of minimum and maximum output amplitude of the solid-line spectrum 
reported in c), respectively. c) Spectrum obtained by Fourier transforming the fields 
as a function of time, calculated at the output of the waveguide with (solid line) and 
without (dashed line) the two rods. d)-g) Transmission spectra calculated by coupled- 
mode theory (solid lines) for increasing cavity resonant frequency. The dashed lines 
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represent the transmission spectrum through the two reflecting elements without the 
presence of the cavity. Figures adapted from [19]. 
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Figure 4.5: a) Schematics of the experimental SNOM setup employed in illumina- 
tion/collection geometry, with dielectric probes for resonant back-scattering experi- 
ments. The part of the setup below the sample measures the degree of polarization 
after the near-field probe, controlled by means of the Babinet-Soleil compensator, 
whose schematics and image are reported in c)-d), respectively. SC is the supercon- 
tinuum laser source, whose spectral emission is reported in b) [33]. LP is a long-wave- 
length pass filter. P1 is a linear polarizer. PBS is a polarizing beam splitter cube. Ob. 
is a 50x, 0.4 NA objective. PB is a narrow band pass filter. PD is a photodiode light 
detector, L a lens that focuses the signal to the entrance slit of the spectrometer (SP) 
and to the detector (DTC) composed by an array of InGaAs photodiodes. The black 
arrow in the incoming light beam and the black dot in the collected light beam high- 
light the cross-polarization between the two branches. The electronic central unit and 
the piezoelectric movements are not shown, but are the same of Fig. (3.2). 
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Figure 4.6: a) RBS raw spectra collected by dielectric probe in the middle of the D2 
nanocavity with lattice constant 321 nm whose SEM image is reported in the inset 
with a scale bar of 300 nm. b) Background spectrum collected about 4 um away from 
the center of the cavity, as highlighted in the inset. c) Spectrum obtained by subtract- 
ing the background b) from the raw data a). d) Background subtracted spectrum aver- 
aged over 6 adjacent positions at the center of the cavity during the SNOM scan. 
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Figure 4.7: a) RBS spectrum of MI acquired during a SNOM scan and background 
subtracted, is reported as a function of the photon energy (black dots) and fitted by 
the Fano formula of Eq. (4.2) (red curve). b) Experimental RBS of M2 acquired during 
a SNOM scan and background subtracted, is reported as a function of the photon en- 
ergy (black dots) along with the Fano fitted lineshape (red curve). The calculated pa- 
rameters and relative uncertain are shown in both insets. 
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Figure 4.8: a) Dots and triangles represent normalized RBS spectra of the mode M1 
collected at two different positions, corresponding to red and green circles in c). Both 
curves are reproduced by Eq. (4.2), which gives the red and the green normalized 
profiles with resonant energies Ey =930.72 meV and E =930.88 meV, highlighted 
by red green vertical dashed lines, respectively. b) Fano amplitude distribution, eval- 
uated as Fy(1 + q?) for the M1 mode. c) Tip induced spectral shift of M1. d) Electric- 
LDOS calculated by FDTD at the M1 peak energy. The dotted orange lines in d) in- 
dicate the positions where the profiles reported in Fig. (4.9) are considered. e) Fano 
amplitude distribution of the M2 mode. f) Tip induced spectral shift of M2. g) Elec- 
tric-LDOS calculated by FDTD at the M2 peak energy. The scale bar is 250 nm. 
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Figure 4.9: a) Comparison between the vertical profiles obtained from the normalized 
map of the spectral shift experimental data (red dots) of Fig. (4.8 c) and the ones 
obtained from the normalized convolution of the calculated electric-LDOS intensity 
map of Fig. (4.8 d) with two-dimensional Gaussian point spread functions of FWHM 
equal to 50 nm, 70 nm and 100 nm. The error bars are provided by the Fano function 
fit of each spectrum acquired. b) Same as a) for the horizontal profile of Fig. (4.8 c 
and d). 
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Figure 4.10: a) RBS background-subtracted spectra of the mode M1 collected in two 
different points, labeled A and B in b) and c). b) RBS intensity distribution evaluated 
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at the peak maximum energy of the M1 Fano resonance, E = 930.93 meV. Withe 
color corresponds to zero intensity. c) RBS intensity distribution evaluated at the E = 
929.70 meV to maximize the color contrast between regions with opposite intensity 
sign. d) Real part of the x component of the electric field calculated by FDTD on a 
detector placed at 22 nm from the slab surface, evaluated at the energy of M1 and 
normalized in [-1; + 1]. e) FDTD calculation of the real part of the x electric field, 
as in d), but in the case of a detector placed at 200 nm distance from the cavity. 
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Figure 4.11: a) RBS background-subtracted spectra collected by the campanile tip on 
three different positions (C, D, and E) shown in the intensity map b), which is evalu- 
ated at E = 931.28 meV. All the spectra recorded in the red (blue) regions have line- 
shapes similar to the C (D) spectrum. In the inset is reported the schematics of the 
campanile tip. c) Intensity data acquired along the central vertical profile between 
adjacent lobes in b). The error bars are provided as the square root of the intensity 
counts and the red curve is the error-function which best reproduces the experimental 
data. d)-e) Spectral shift map and amplitude map, evaluated as F (1 + q?), for the 
same SNOM measurement shown in b). The scale bar is 250 nm. 


106 


Niccolò Caselli 


o 
© 


PBS 
Polarization Hem 


Compensator 


q 
o 


4 


A 


fad 
© 


RFS Intensity [Kc/s] 
ò o 


ob. | 
SC Laser P Pi 
= Energy [meV] 


Figure 4.12: a) Schematics of the experimental SNOM setup in transmission geometry 
and equipped with coated probes or campanile tips for resonant forward-scattering 
experiments. All the optics and the instruments employed here are also present in the 
illumination/collection setup of Fig. (4.5 a) for RBS measurements. SC laser light 
linearly polarized illuminates the sample from the far-field by the objective (Ob.); the 
forward-scattering signal is analyzed in cross-polarization scheme and acquired by 
the spectroscopy detection apparatus. SEM image of the investigated nanocavity is 
reported in the inset along with the positions (J, K, L) where the background-sub- 
tracted spectra in b) are detected by aperture coated probe. 
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Figure 4.13: a) Schematic of the FDTD calculations modelling the RFS imaging. The 
PCN is excited from the bottom by a linearly polarized plane wave and three detectors 
(green squares), centered on the lobes indicated by a, p and y in c), are placed 10 nm 
above the sample surface. b) Calculated RFS spectra collected by the three detectors 
in a); every spectrum is background subtracted and divided by its own maximum 
value. c) Real part of the x component of the M1 electric field calculated by FDTD 
10 nm above the sample surface. d) Spatial distribution of the RFS intensity evaluated 
at 932.0 meV. J, K and L correspond to the position highlighted in Fig. (4.12 a). The 
inset represents the schematic of the coated probe. e) Intensity data acquired along the 
central vertical profile between adjacent lobes in d). The error bars are provided as 
the square root of the intensity counts and the red curve is the error-function that best 
reproduces the experimental data. The scale bar is 250 nm. 


108 


Niccolò Caselli 


a 
) 8 (O) M2y 
4 
È 
O 
= 
z|... 0 
Val 
= 
g S @ M2y 
= 
n 
LL 
E å 
0 


944 946 948 950 952 944 946 948 950 952 
Energy [meV] r Energy [meV] 


Figure 4.14: a) RFS spectra background-subtracted of the excited mode M2 collected 
by the Al-coated probe in the positions M and N reported in c). As highlighted by the 
schematic in b) the incoming light has a linear polarization along the y axis, and in 
the detection path no cross-polarization scheme is mounted. c) Spatial distribution of 
the experimental intensity of the M2 y-component evaluated at 947.2 meV. d) Real 
part of the electric field relative to the M2 y-component calculated by FDTD and nor- 
malized between [—1; + 1]. e) RFS spectra background-subtracted of M2 collected 
by the Al-coated probe in the positions P and Q reported in g) for the case of x-polar- 
ized incoming light, as reported in f). g) Experimental spatial distribution of the x- 
component intensity evaluated at 947.2 meV. h) Real part of the electric field relative 
to the M2 x-component calculated by FDTD and normalized between [—1; + 1]. The 
scale bar in all the maps is 250 nm. 
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Chapter 5 


Photonic crystal molecules 


The coupling between adjacent optical nano-resonators is a topic of vast attraction 
both for fundamental physics and for high density photonic integrated circuits appli- 
cations. Interacting photonic nanocavities emerge as analogous of the basic inter- 
atomic coupling that gives rise to real molecules, described by wave-functions with 
defined symmetry character. They are called photonic molecules. The high spatial 
resolution imaging of the photonic molecules resonant modes it is a primary task to 
be achieved. Here we describe the basic principles of the interaction between coupled 
resonators and we study, both experimentally and theoretically, the light behavior in 
a photonic molecule formed by two and three interacting photonic crystal nanocavities 
aligned along different symmetry axes. We notably find that the photonic molecule 
ground state can be either bonding or antibonding and that the effective mode inter- 
action cannot be limited to the nearest-neighbour. Moreover, we develop a post-fab- 
rication method to tune the PCN localized modes by inducing a laser assisted nano- 
oxidation that produces a large blue-shift of the modes. This technique is also ex- 
ploited to engineer the coupling strength in photonic crystal molecules and to prove 
the symmetry character inversion between the molecule ground state and the first ex- 
cited state. 


5.1 Coupled photonic nanocavities 


Dielectric point defects in photonic crystals can generate localized light states in 
the photonic band gap, similarly to the case of electronic states in the band gap of 
semiconductor materials, as we discussed in Sec. (2.1.2) [1]. The light states localized 
in photonic crystal nanocavities play a relevant role in the development of light emit- 
ters and optical waveguides [2]. Extending the analogy between photonic modes res- 
onant in single PCNs and electrons described by quantum mechanics in atomic-like 
states, coupled PCNs are also denominated photonic crystal molecules, due to an ef- 
fective interaction between close nano-resonators [3-5]. Light behavior in systems 
based on aligned PCNs resembles the molecular interaction where the resulting nor- 
mal modes exhibit energy splitting. This condition is achieved by an evanescent pho- 
ton tunneling between each single PCN, whenever the resonant wavelength matching 
and the electric LDOS spatial overlap between them are fulfilled. Photonic molecules 
have been implemented on the basis of any kind of optical nanoresonator, ranging 
from the case of nominally identical micropillars, microdisks or dielectric spheres that 
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sustain whispering gallery modes, to the case of coupled identical or even non-iden- 
tical PCNs, which lead to the formation of homo-atomic or hetero-atomic photonic 
crystal molecules, respectively [3, 6-10]. In the simplest and remarkable case of only 
two coupled PCNs, following quantum mechanics denominations, the diatomic mo- 
lecular state that is symmetric in the exchange of the two atomic orbitals is named 
bonding (B) state and it corresponds to the ground state in the case of attractive (neg- 
ative) coupling. While the antisymmetric molecular state is denominated antibonding 
(A) and it is associated to the ground state only if the coupling is positive. Interest- 
ingly, the molecular-like modes exhibit different radiation patterns, thus engineering 
the coupling strength is proposed as an alternative method to tailor the angular distri- 
bution properties of the resonant modes. Electromagnetically induced transparency 
like effects can occur in coupled optical resonators due to the coherent interference, 
as well as slow-light phenomena [11-13]. In quantum-information application it has 
been proved that a photonic molecule can enhance the spontaneous emission rate of 
an embedded single quantum dot, which emits a single pair of entangled photons, by 
ensuring a deterministic coupling of the molecule normal modes to the biexciton and 
exciton energy transitions, thus improving the source of entangled photons [14]. Pho- 
tonic molecules would also represent the basis for developing two-qubit gates in inte- 
grated structures [15]. Also the possibility of using strongly coupled photonic mole- 
cules to efficiently increase the mutual interaction between atomic-like light emitters 
(quantum dots), even at large distance, has been theoretically addressed [16]. 

In a given resonator the strong coupling regime between a single emitter excita- 
tion and the cavity mode gives rise to a dressed-photon state (polariton). This behavior 
induces non-linear optics phenomena at extremely low light intensities, such as pho- 
ton blockade and photon induced tunneling [17]. Within this framework, a fascinating 
implementation has lead to the proposal of an optical analogue of the Josephson in- 
terferometer consisting of two coherently driven linear optical cavities connected 
through a central cavity with a single-photon non linearity [18]. The big advantage of 
the system composed by interacting resonators relies in the accessibility to local prop- 
erties such as the possibility to measure the output of a single cavity. Recently, com- 
plex architectures formed by large ensemble of strongly interacting dressed-photonic 
molecules have been proposed to address quantum many-body aspects of quantum 
gases, as the Mott insulator to superfluid transition, resembling the case of ultracold 
atoms in optical lattices [19, 20]. In order to create proper quantum optics devices a 
fundamental requirement is the design and control of adjacent nanocavity modes at 
the target wavelengths, within an accuracy which is not directly obtainable due to the 
fabrication tolerances. The compensation of the fluctuations related to the structural 
disorder and more generally the control of the resonance wavelength of each resonator 
within the array and also of the tunneling coefficient between adjacent nanocavities 
is of the utmost relevance in the roadmap toward quantum photonic devices. 


5.1.1 Coupled modes theory 


The analytical expression of the coupled PCNs mode features can be modeled by 
considering a system of two oscillators of resonant energy hw, and hwy, respectively, 
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whose interaction is mediated by a potential g that defines the coupling strength. The 
problem can be solved by calculating the eigenvalues and eignevectors of the matrix: 


(e) DI 


whose eigenvalues solutions, considering also the damping (where the resonant 
mode losses are yı and y, for the single oscillators and F for the supermodes), are 
given by: 


(5.2) 


where Q_ and Q, represent the molecular-like ground state and excited energy 
state, respectively. While the corresponding eigenvectors on the basis of the single 
oscillators are: 


w= (“eee tee 1) (5.3) 
Y, = (Motes he tee 1) (5.4) 


29 i 


where Y_ and Y, are the molecular states relative to Q_ and Q4, respectively. In 
the limit of negligible detuning, that is |A] = |hw, — hw2| « g, the lower energy 
molecular state, P_, can be either an even state (1; 1) if g < 0, or an odd state (-1; 1) 
if g > 0. While Y, has always an opposite parity with respect to the V_ state. 

Moreover, if we consider g < 0, as happens for attractive potentials, and a negli- 
gible detuning, the molecular states result: 


W_==(91 +p) (5.5) 
Y, == (91 - 92) (5.6) 


where 1 and @z are the wave-function of the two isolated oscillators, respec- 
tively. In particular, the molecular ground state eigenfunction (5.5) is given by a linear 
symmetric combination of the atomic-like states, which is called bonding state in anal- 
ogy with real molecules behavior. While the eigenfunction (5.6) highlights that the 
excited state is formed by an antisymmetric combination, giving rise to the antibond- 
ing state. However, if the effective coupling strength is positive, an inversion between 
the molecular modes eigenfunctions of (5.5) and (5.6) occurs. In fact, the ground state 
acquires an antibonding character and the excited state becomes bonding. 


In photonics the role of the wave-function is accomplished by the e.m. fields while 
the effective potential is related to the dielectric environment [1]. Since the e.m. fields 
are vectors, we need to separate them in each polarization component in order to get 
a fruitful comparison with the scalar theory that we presented. For instance, for a given 
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localized mode the symmetry character can be defined only for one single field polar- 
ization component and this information cannot be retrieved by means of PL measure- 
ments, independently on the spatial resolution, since they allow only electric and mag- 
netic field intensities detection. In order to probe the symmetry character of the cou- 
pled modes a revisited Young’s double slit experiment in the far-field has been tested 
with photonic molecules [21, 22]. The experimental results can be compared with 
accurate numerical calculations of the field, exhibiting both amplitude and phase in- 
formation. Moreover, resonant scattering phase-sensitive measurements on photonic 
molecules are still lacking and they will be implemented in the future, employing the 
Fano-imaging setup proposed in Sec. (4.2.2). 

The expression of the coupling strength in photonic molecules can be derived 
from the tight binding theory, in the approximation of effective electric interaction, 
see also the Appendix [23, 24]: 


__ 2rhc  f Ae(r)Eo(r)-Eo(r+D)dr 57 
A flmolHo(r)|? +e) |Eo(r) |? ]ar e 
where Eo and Ho are the electric and magnetic fields of the single PCN; D de- 

scribes the spatial separation between adjacent cavities; Ae is the difference between 
the dielectric constant of the coupled system and the dielectric constant of the single 
resonator. The relation (5.7) highlights how the photon tunneling is governed by the 
electric LDOS spatial overlap between neighboring resonators. It depends on both the 
system geometry and on the dielectric permittivity values. 


5.1.2 The D2 cavity-based photonic molecules 


Starting from the D2 photonic crystal nanocavity, considered as the atomic-like 
building block of complex coupled systems, two kinds of photonic molecules have 
been realized and investigated [6]. The one in which two D2 cavities are aligned along 
the M photonic crystal axis, and the other where the nano-resonators are aligned along 
the K symmetry direction, reported in the insets of Fig. (5.1 a) and Fig. (5.2 a), re- 
spectively. In the spectrum of M-aligned molecule, shown in Fig. (5.1 a), are present 
two pairs of peaks, P1-P2 exhibiting a large splitting and P3-P4 with a small spectral 
separation, comparable to the mode linewidth. All of them are roughly in the same 
spectral region of the M1 and M2 modes of the single D2 cavity. By analyzing 
Eq. (5.2), neglecting the losses, we find that the molecular mode splitting (Q, — Q_) 


reduces to yA? + 4g. Then, the splitting between molecular modes is proportional 
to the coupling strength but it is also affected by detuning. Therefore, the presence of 
a double peak in the spectrum is not enough to state the presence of an effective cou- 
pling. More insight is pursued by the electric field intensity distribution analysis, per- 
formed by mapping for each mode the wavelength spectral shift induced by the 
SNOM dielectric probe, as reported in Fig. (5.1 b-e). P1 and P2 intensity distributions 
covers the whole molecule. In both maps is clearly recognizable the superposition of 
M1 mode intensity distribution on the two PCNs that constitute the molecule. There- 
fore, the combined spectral and spatial information makes clear that P1 and P2 are 
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molecular-like modes resulting from a strong coupling between the M1 modes of the 
single M-aligned PCNs. For the investigated sample, also the maps concerning P3 and 
P4 are delocalized on the entire molecule and are formed by a superposition of the 
M2 modes of the single PCNs, as reported in Fig. (5.1 d-e). Thus they result as mo- 
lecular-like modes, whose spectral separation is ascribed to a small but effective cou- 
pling. The behavior of the M-aligned molecule is predicted by Eq. (5.7). In fact, a 
large overlap between the isolated PCN modes leads to a strong coupling, as happens 
for MI that is elongated parallel to the M-axis resulting in P1 and P2 molecular-like 
modes. While a small overlap results in a reduced coupling that can be in competition 
with the disorder induced detuning, as it happens for the single cavity mode M2, 
which is elongated parallel to the K-axis and gives rise to P3 and P4 molecular modes. 

Conversely, K-aligned molecule exhibits the opposite behavior with respect to the 
M-aligned molecule. In fact, due to the spatial properties of M1 and M2, the coupling 
between P1-P2 is smaller with respect to the interaction between P3-P4. In the inves- 
tigated sample, whose spectrum is reported in Fig. (5.2 a), PI and P2 modes show a 
small splitting. But the evidence of a negligible coupling, with respect to disorder 
induced detuning, is found only after the analysis of the electric field intensity maps 
of Fig. (5.2 b-c). Where it is clear that both P1 and P2 modes are localized only on 
one PCN and, therefore, instead of molecular modes they have to be considered as 
atomic-like modes. On the other hand, the widely spectrally separated modes P3 and 
P4 experience a strong coupling characterized by spatial distributions delocalized on 
the whole system, see Fig. (5.2 d-e). Therefore, in the investigated molecule the split- 
ting between P1 and P2 modes has to be mainly ascribed to disorder induced detuning 
(g < A), while the splitting between P3 and P4 modes is largely attributed to an ef- 
fective interaction (g > A). 

The straightforward result of the angular emission pattern analysis, performed in 
the far-field for each polarization, is the symmetry character probing of the strongly 
coupled modes in both M-aligned and K-aligned D2-based photonic molecules [21]. 
In particular, the M-alignment is found to sustain a bonding ground state (P1), conse- 
quently exhibiting g < 0. On the other hand, in the K-alignment the molecule lower 
energy state symmetry character cannot be retrieved since the P1-P2 weak coupling 
can be deeply affected by the detuning between the two PCNs and, therefore, P1 can 
result in an atomic-like state, as reported in Fig. (5.2 b). 


5.2 Post-fabrication control of evanescent coupling 


Quite remarkably, in Sec. (5.1.2) it is shown that in photonic molecules the sign 
of the interatomic-like coupling is not uniquely defined and so is the parity of the 
modes. It is also proved, by the analysis of the far-field emission, that it is possible to 
change the character of a given molecular mode from even to odd, depending on the 
lattice direction along which the two PCN are aligned [21, 25]. Naturally emerges the 
question if in photonic molecules a fine tuning of the coupling strength expressed by 
Eq. (5.7) is achievable both at the design stage and with post-fabrication techniques, 
changing from positive to negative values, thus leading to an antibonding ground 
state. Firstly, we exploit accurate numerical FDTD-based calculations to neglet the 
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unavoidable fabrication induced disorder that leads to a significant detuning between 
the two PCNs. Such analysis allows us to demonstrate the possibility to tune the 
ground state of a given photonic molecule from a bonding to an antibonding character, 
by operating a local modification of the dielectric environment in between the two 
coupled nano-cavities. This also leads to possible implementation of a fine control of 
the photon tunnelling in coupled systems. Secondly, as a proof of principle, we ex- 
perimentally probe the antibonding ground state nature along with the coupling 
strength induced variation of the K-aligned molecule. In order to achieve this goal we 
develop a local post-fabrication technique based on laser assisted non-thermal nano- 
oxidation [26]. 


5.2.1 Engineering the mode parity 


For the numerical characterization of the resonant modes in photonic crystal mol- 
ecules we employ the three-dimensional FDTD solver package presented in 
Sec. (2.1.3). We consider a 320 nm-thick semiconductor membrane with a refractive 
index n = 3.484, the same of GaAs material, but we could have chosen Si or SizN4 
as well. The photonic structure is a two dimensional hexagonal lattice of holes (n = 
1) with lattice constant of a = 308 nm, the pores have a diameter size d = 193.2 nm, 
leading to a 35% filling fraction. Since the unit-cell by whom the dielectric matrix is 
discretized has size (x = a/28,y = V3a/28) matched with the hexagonal lattice, 
every computed hole is exactly equal to the others and no detuning effect are present. 
The near-field detector is placed parallel to the PCN surface at a distance of 22 nm 
and covers the whole molecule. 

The simulated photonic molecule is formed by two D2 nanocavities aligned along 
the K axis, since this configuration allows a small but sizable coupling between the 
MI modes of the single PCN, resulting in the P1, P2 modes. In fact, the calculated 
spectrum of the nominal structure without detuning, reported in Fig. (5.3 a), shows a 
double peak structure with an energy splitting of the order of 1 meV. Neglecting the 
losses in Eq. (5.2), the molecular mode splitting hO, — AQ_ reduces to 2g for A = 0, 
thus the splitting between P1 and P2 in FDTD spectrum gives a direct estimation of 
the coupling strength absolute value. Also the FDTD calculated electric field spatial 
distributions, reported in Fig. (5.3 b-d), are in agreement with the existence of a small 
but effective coupling, since both modes are delocalized over the molecule. The 
modes arising from M1 of the single PCN are linearly polarized along the cavity align- 
ment direction, thus facilitating the numerical analysis since only one polarization 
component has to be investigated. Since M1 exhibits an even electric field distribution 
along both x and y axes, as reported in Fig. (4.10 c), the analysis of the maps in Fig. 
(5.3 b-c) points out that PI mode has an antibonding character showing an odd spatial 
distribution with respect to the alignment axis inversion. While P2 mode is a bonding 
molecular mode with an even spatial distribution along x. Therefore, as highlighted 
by the electric field distribution of Fig. (5.3 b) and differently from real molecules, 
the K-coupling supports an antibonding ground state. An important outcome of the 
antibonding mode character is that the far-field emission is exactly zero along k, = 0 
and in particular at normal direction, as shown in Fig. (5.3 c). This is a consequence 
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of destructive interference between the two out of phase coupled single PCNs [21]. 
For similar reasons, the bonding mode P1 has a strong emission along the normal 
direction, due to the constructive interference of the two in phase PCNs along k, = 
0, as reported in Fig. (5.3 c-e). Therefore, an inversion in the mode parity and sym- 
metry can determine a large variation in the angular distribution of the photonic mol- 
ecule modes. 

As it is pointed out in Eq. (5.7) the character (B or A) of the ground state depends 
on the overlap integral between the two atomic-like modes weighted over Ae. In par- 
ticular the interference of the two single cavity fields can result either in a positive or 
in a negative g, which induces the molecular ground state to be A or B, respectively. 
In fact, the single PCN localized mode Eo exponentially decays moving away from 
the cavity with typical oscillation between positive and negative values superimposed. 
Therefore, the constructive interference between the two Eo localized fields can pro- 
vide g > 0 and the antibonding ground state character. 

It was also experimentally observed, for both the investigated polarization, that 
the K-aligned molecule sustains an antibonding lower energy state (P3) arising from 
the strong coupling between the M2 mode of the single D2 cavity, that give rise to P3 
and P4 modes [21]. Analogous behavior can be also found in other kinds of nanocav- 
ity based photonic molecule [16, 22, 25, 27]. On the contrary, investigating the overall 
K-aligned molecule ground state (P1) is a difficult task. Since P1 and P2 modes are 
weakly coupled, a post-fabrication technique, able to compensate the disorder induced 
detuning, is required before measuring the coupled mode symmetry character. 

Nevertheless, finding an antibonging lower energy molecular state is a counterin- 
tuitive result with respect to the standard expectation on the basis of the quantum me- 
chanical analogy, for instance with standard coupled quantum wells. However, an an- 
tibonding ground state has been recently observed in coupled InAs quantum dots, 
where the heavy-hole light-hole mixing and the spin-orbit interaction causes this 
counterintuitive reversal for particular values of the barrier thickness [28]. An alter- 
native description of the physics involved in photonic molecules, can be derived by 
recalling that the eigenstates of an arbitrary photonic system has to satisfy the varia- 
tional principle, related to the minimization of the energy functional U [1]: 


2 
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The ways to minimize Eq. (5.8) are two: (i) by reducing the number of mode 
nodes and (ii) by increasing the field concentration in the high-dielectric regions. As 
in quantum mechanics, bonding states always satisfy the condition (i) and this is the 
reason why the ground states usually have bonding character. However, while the 
wave-function in atoms decays exponentially with distance from the nucleus, the lo- 
calized electric field in a PCN shows an exponentially damped oscillation with dis- 
tance from the local dielectric defect. The overlap between the oscillating evanescent 
electric field of the two single PCN modes produces constructive or destructive inter- 
ference. Obviously, in the case of constructive (destructive) interference, the electric 
field is enhanced (reduced), and this decreases (increases) U. Therefore, the presence 
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of an antibonding ground states in photonics is associated with constructive interfer- 
ence in the defect region for the odd state, which may enhance the electric field con- 
centration in the high-dielectric regions. 

An open question is whether it would be possible to control both the photonic 
coupling strength value and sign by local modifications of the structure design. This 
adjustment has to be symmetric with respect to the molecule geometry, otherwise it 
induce a sort of detuning between the two single cavities. The region of strong overlap 
between the single molecular modes E(r) and Eg(r + D) can be evaluated by map- 
ping the value: 


\Eg(r) |? — |Ea(r) |? (5.9) 


which for the considered K-coupled photonic molecule is reported in Fig. (5.4 a). 
The overlap is distributed on the whole system, but the maximum value is close to the 
four pores around the central hole between the two cavities. Then we reduce the di- 
ameter size of five central holes of the molecule, highlighted in Fig. (5.4 b), imposing 
a symmetric perturbation. Although such perturbation affects a region where Ae = 0, 
the coupling strength given by Eq. (5.7) is modified by a variation of the single cavity 
resonant mode Eo equally for both PCNs, as assured by the symmetric dielectric al- 
teration. This perturbation leads to a marked tuning of both molecular modes but with 
different strength. The peak energy values of the B and A modes are reported in 
Fig. (5.4 c) as a function of the modified holes diameter size d. Starting from the nom- 
inal d (193.2 nm), which corresponds to an A ground state with 1.0 meV energy split- 
ting, we see that the diameter size reduction produces a continuous decreasing of the 
photonic coupling, directly evaluated from the mode splitting. For d = 183.2 nm the 
molecular spectrum shows a single peak, meaning that the photonic coupling is de- 
creased below the mode broadening. This leads to a degeneracy of the two molecular 
modes that is purely accidental, since the C2, symmetry of the K-coupled D2 photonic 
cavities cannot support degenerate modes accordingly to the group theory [29]. With 
a further reduction of d the ground state changes into a bonding character, as a con- 
sequence of the change of the photonic coupling sign. For d = 173.2 nm the mode 
splitting is 1.1 meV, almost equal to the nominal design, but with modes exhibiting a 
reversed parity. Therefore, a smooth transition of the lower energy state from bonding 
to antibonding can be achieved. 

A more intuitive picture of the physics underlying the data reported in Fig. (5.4 c) 
can be obtained by considering the dielectric perturbation induced tuning in the basis 
of the molecular modes. In fact, decreasing d means adding dielectric material, ac- 
cordingly with the observed red shift of both modes. In addition, the perturbation be- 
ing very local, the mode shift is larger the higher is the mode electric field intensity at 
the dielectric modification position. Therefore, the presented local perturbation results 
in a more pronounced red-shift of the B mode with respect to the A mode. 

Further information about the photonic molecular modes parity transition is ob- 
tained by analyzing the mode quality factors as a function of d, as shown in Fig. 
(5.4 d), where the mode broadening of both P1 and P2 modes are reported as a func- 
tion of d. It is well known that the formation of the molecular modes with opposite 
parity also tailors the diffraction losses [30]. For the nominal structure with identical 
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holes the lower energy broadening is 0.38 meV, while a smaller broadening 
(0.22 meV) is found for P2, in agreement with the far-field patterns of Fig. (5.3 c and 
e), which show that P1 has larger diffractive losses. These values are almost constant 
when reducing d from the parity inversion, then they reach an intermediate value 
when the coupling vanishes and finally they are exchanged after the transition. In 
other words, we find that the Q of the A and B modes are almost unaltered. Moreover, 
the A mode is always broader than the B one, in contrast to the intuitive findings for 
coupled L3 cavities, associated to the strong (weak) emission of the B (A) mode in 
the vertical direction [30]. However, the discrepancy is only apparent, since always 
the broader mode has larger diffractive losses in the light cone. The different behavior 
is simply due to the differences in the angular patterns of the single modes, as L3 
cavity is compared to the D2 cavity. Therefore, the broadening exchange between P1 
and P2 is a consequence of the parity transition achieved with the manipulation of five 
central pores. 

The counterintuitive change of the photonic molecule ground state from a bonding 
to an antibonding orbital character is also confirmed by performing numerical guided 
mode expansion method calculations [31]. In Fig. (5.5 a-b) are reported the resonant 
energy of both bonding and antibonding modes in M-aligned and K-aligned photonic 
molecules, respectively, as a function of the cavity-to-cavity distance. For these cal- 
culations, we use a supercell dimensions of up to 16a X 7wy for the K-alignment 
(where wo = V3a), and 10a x 13w for the case of M-alignment. The photonic mol- 
ecules discussed so far are reported in the insets of Fig. (5.5 a-b) and correspond to a 
center-to-center distance equal to 3V3a and 3a for the M-aligned and K-aligned pho- 
tonic molecules, respectively. Figure (5.5 a) shows that the ground state for the M 
coupled molecule, for increasing center-to-center distances, is always bonding and 
that the modes energy splitting decreases monotonically, as expected for standard 
coupled systems. On the contrary, in the case of K-aligned molecules the ground state 
symmetry does not follow a trivial behavior, as results from Fig. (5.5 b). In fact, the 
lower energy mode has either an even or an odd character as a function of the cavity 
separation. Furthermore, the small energy mode splitting shows little variations for 
center-to-center distance greater than 3a. These results agree with experimental data 
and FDTD calculations thus confirming the subtle nature of the photonic coupling, 
driven by a complex interference between localized light states. 

It could be interesting to explore whether also a post-fabrication tuning method 
would be able to control the molecular mode parity, thus it could be used to used to 
realize optical switches exploiting the inversion of the peculiar far-field patterns of 
the A and B modes. In order to propose experimental methods to achieve a post-fab- 
rication coupling control, we consider a gentle modification of the refractive index 
within a circular central region of the photonic molecule with a diameter of the order 
of 0.8 um. To address the problem, a controlled reduction of the effective refractive 
index could be experimentally obtained either by a local photo-induced oxidation of 
the semiconductor slab that will be presented in Sec. (5.2.2); or dynamically by in- 
jecting free carriers, for instance with a fs pulsed laser [32-34]. 

We theoretically implement such techniques in the case of the modified photonic 
molecule with the five central hole diameter sizes reduced to 173.2 nm. Since starting 
from such geometry we can obtain the opposite ground state character transition to 
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the one reported in Fig. (5.4 c). The results are reported in Fig. (5.6), where the sche- 
matics of the modified photonic molecule are drawn as left panels and the correspond- 
ing calculated spectra are reported in the right panels. Starting from the nominal re- 
fractive index equal to 3.484, we reduce the refractive index of the circular central 
part of the underlying slab at 3.416. The corresponding calculated spectrum shows a 
single peak only, denoting a very small (if any) photonic coupling. Eventually, for a 
refractive index of 3.348, we find that the P1 character is changed into antibonding 
and that the spitting reaches 1.2 meV. 

In summary, it is possible to tune a photonic molecule mode from the bonding to 
the antibonding character and vice versa. The control can be experimentally achieved 
by combining a careful design of the central region between the two PCNs or by post- 
fabrication techniques such as local laser-assisted oxidation or dynamic free-carriers 
injection exploiting fs pulsed laser, for instance. These properties would open the road 
to possible applications in quantum optical information processing, including the op- 
tical switch of the far-field pattern of quantum emitters. 


5.2.2 Photoinduced non-thermal oxidation tuning 


Concerning most of the PCNs applications in optoelectronic and quantum optics, 
a fundamental requisite relies on the design and control of the localized modes at the 
target wavelengths, within an accuracy which is not directly obtainable due to the 
unavoidable fabrication tolerances [35-37]. Different post-fabrication processing 
methods, able to tune the nanocavity modes for compensating the fabrication imper- 
fections, have been demonstrated in the last few years, such as tip and thermal induced 
tuning, wet chemical digital etching, atomic layer deposition, AFM induced nano- 
oxidation, liquid nano-infiltration, and nano-electromechanical tuning induced by 
controlling the applied voltage to PCN based devices [32, 38-43, 25]. The combina- 
tion of these techniques allows large and fine tuning of the resonant mode, with re- 
spect to the mode linewidth, and enables both a reversible or irreversible mode spec- 
tral shift towards lower or higher wavelength. 

Recently, a very simple method based on the thermal oxidation by laser irradiation 
in air atmosphere has been realized in GaAs and Si structures [33, 44]. This method 
has the great advantage of being local without needing extra materials and processing 
tools. However, it requires a large thermal budget which can produce undesired an- 
nealing and modification of the underlying semiconductor layers. In order to over- 
come this drawback, we implement a method to achieve photo-induced fine-tuning of 
PCN resonant modes with very low thermal budget. It is based on light induced oxi- 
dation in air atmosphere which produces a local reduction of the GaAs membrane 
effective thickness and at the same time an enlargement of the holes diameter effective 
size resulting in a large and smooth blue-shift of the nanocavity modes. 


The experimental setup employed to analyze the photo-induced oxidation is a 
confocal apparatus where a 50X microscope objective (NA = 0.7) replaces the near- 
field probe in Fig. (3.2). The sample is excited either with light from a diode laser 
(780 nm) or by an Ar ion laser (514 nm), and the emitted PL signal is coupled to a 
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single-mode optical fibre of 6 um core diameter, acting as a confocal pinhole con- 
nected to a spectrometer with a 300 lines/mm grating that provides a spectral resolu- 
tion of 0.5 nm. Finally, the sample is mounted on a xy stage so that it can be scanned 
in respect to the objective, allowing the collection of two-dimensional PL maps. 
Therefore, both the nanocavity mode observation and the laser-driven nano-oxidation 
are performed in the far-field to validate a proof-of-principle that can be extended, as 
we will show, to SNOM configuration for improving the induced oxidation resolution. 
The spatial resolution obtained by the confocal setup is evaluated performing a PL 
intensity map of the fundamental PCN mode M1 at low excitation power (P = 
0.1 mW) with the 2 = 514 nm laser, where heating can be neglected, as reported in 
Fig. (5.7 a). The laser spot FWHM size results 1.1 um thus defining the spatial reso- 
lution. The local temperature of the samples under laser exposure is directly measured 
by using the PCN modes as a local thermometer and taking advantage of the thermal 
red-shift of the resonant modes that, around room temperature, has been calibrated for 
the M1 mode of the D2 nanocavity to be 0.12 nm/K [45]. 

By increasing the laser density power illuminating and collecting PL signal in the 
center of the PCN we observe a sizeable red-shift, as reported in Fig. (5.7 b) for 
780 nm laser excitation. To get more insight into the temperature induced distribution 
the thermal induced spectral shift map and the average radial profile are evaluated in 
Fig. (5.7 cd), respectively. The map shows that the thermal gradient extends over a 
2 um region, that is slightly larger than the excitation spot and allows us to retrieve 
the temperature distribution. However, during laser exposure, the PCN modes may 
shift due to oxidation, optical Kerr effect, free carrier injection, and heating. Oxidation 
leads to an irreversible blue-shift which is negligible for the experimental time expo- 
sure (100 ms) corresponding to the integration time of the spectra acquired in 
Fig. (5.7). Optical Kerr effect occurs at much higher excitation power, due to the small 
non-linearity of GaAs [46]. The tuning due to free carrier injection gives a blue-shift 
estimated to be of the order of 1 nm for the data reported in Fig. (5.7 c) by using the 
parameter given in ref. [34]. Heating gives a red-shift of the modes and, therefore, in 
the reported experimental condition proves to be the dominant effect. Still the correc- 
tion for the blue-shift (of the order of 20% of the experimental red shift) due to free 
carrier injection has to be accounted for an accurate estimation of the sample temper- 
ature. The experimental spectral shift of M1 at the center of the spot in Fig. (5.7 c) is 
4.7 nm, corresponding (including the estimated blue-shift for carrier injection) to a 
temperature rise as low as AT < 47 °C. The temperature increase is large enough to 
be carefully measured but extremely low for thermal effects such as annealing or ther- 
mal oxidation to be activated. Data for green light from the Ar laser are very similar: 
at the maximum power used (P = 0.7 mW), we observe a temperature increase as 
AT = 22°C. 

In order to verify that our data reflect the real temperature of the sample, we use 
FEM calculations to estimate the laser heating for the realistic structure in the exper- 
imental conditions. The temperature profile during laser irradiation is calculated by 
solving the heat conduction equation, taking into account the temperature-dependent 
thermal conductivity of bulk GaAs and the power lost due to reflection at the top 
surface of the membrane. We note that the actual thermal conductivity in the PCN 
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region may be lower than the bulk values due to phonon scattering at the etched in- 
terfaces. Heat losses due to convection and conduction in air are neglected. Multiple 
reflections inside the membrane are also not considered; for green light, the absorption 
in a single pass through the membrane is more than 90%, and therefore, the approxi- 
mation is valid. 

The spatial temperature gradient calculated by FEM at the top surface of the PCN 
at a laser power P = 0.7 mW and 4 = 516 nm (because of the availability of the pa- 
rameters required for the calculation) is reported in Fig. (5.8 a-b). It nicely agrees with 
the experimental data, giving a thermal distribution extending over 2 um. In particu- 
lar, the temperature reaches the maximum value (AT = 16 °C) at the center of the 
laser spot and drops while moving away from it; thus, experimental data and theoret- 
ical simulations nicely fit together. On the other hand, in the case of infrared light 
assuming an excitation power of 2.7 mW and A = 780 nm, the FEM calculation gives 
a maximum heating equal to 21 °C. This value is underestimated because the mem- 
brane absorption in a single pass is only 38% of the incoming power and, therefore, 
neglecting the multiple reflection is likely to be incorrect. Nevertheless, the reliability 
of the temperature estimation and especially the very low thermal budget used in the 
experimental condition are supported by FEM calculations. 

Despite the very low thermal budget of the exposure condition, we observe that a 
non-reversible modification of both the resonant mode wavelength and the PCN die- 
lectric environment occurs after laser illumination of about half an hour under high 
excitation power, that is, P = 2.7 mW for 2 = 780 nm and P = 0.7 mW for A = 
514 nm. In fact, when reducing the laser power and after waiting for the cooling of 
the sample, the PCN modes undergo sizeable blue-shift. The mode shift can be con- 
trolled either by varying the laser power or by iterating laser exposure steps. Figure 
(5.9 a) displays PL spectra acquired at low power (0.1 mW) after different exposure 
times with high power (0.7 mW) green Ar laser, showing a large blue-shift of both 
resonances without remarkable Q or intensity variations. To highlight the time de- 
pendence for both excitation lasers, the resonant wavelength variations of M1 are 
given as a function the square root of the time in Fig. (5.9 b). Although the 780 nm 
excitation induces a larger temperature rise, the comparison clearly indicates that the 
rate of blue-shift is much larger when he high power exposures are performed with 
A = 514 nm light, thus denoting a relevant wavelength dependence on the effect. The 
time evolution of the irreversible effect is a key aspect to assess its origin. Clearly, 
from the spectra of Fig. (5.9 a), it is evident that the blue-shift is not linear with the 
laser exposure time. But, as denoted in Fig. (5.9 b), it is linear with the square root of 
the time, thus indicating a fingerprint of a diffusive process [44]. Therefore, we attrib- 
ute the progressive mode shift with irradiation time to a non-thermal oxidation of the 
GaAs membrane. The oxidation acts by reducing the membrane effective thickness 
and/or by increasing the effective photonic crystal pores diameters. In both cases, the 
effect on the PCN mode is a shift toward smaller wavelengths. In order to demonstrate 
that the blue-shift depends on the presence of oxygen in the environment of the sample 
during laser irradiation, we perform the same confocal experiment in a controlled ni- 
trogen atmosphere. In this case, the blue-shift of the resonant modes is missing, indi- 
cating that oxygen is requisite for observing the PCN modes spectral shift. 
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It is well known that the local GaAs oxidation produces oxide protrusions due to 
the larger molar volume of the oxide with respect to the GaAs. A topography image 
of the D2 cavity, collected by a dielectric SNOM probe acting as shear-force sensor, 
after a laser exposure corresponding to an overall shift of 13 nm of the photonic modes 
is shown in Fig. (5.10 a). The presence of an oxide protrusion in a spatial region cor- 
responding to the laser spot is clearly observed and the central horizontal line profile 
is displayed in Fig. (5.10 b) to depict the height of the oxidized area. A value of about 
20 nm height is found with respect to the non oxidized region. This data could in 
principle be used to estimate the structural changes of the membrane. From the ratio 
of the molar volume of stoichiometric oxide, found in ref. [47], for the investigated 
D2 nanocavity we expect a reduction of the membrane thickness of the order of 10 nm 
on the top interface, and this could be used for calculating the mode shift. 

The final point of our oxidation analysis is to address the physical mechanisms 
for explaining the laser assisted non-thermal oxidation of the GaAs. Effects on the 
GaAs surface after light illuminations have already been reported when the energy of 
the incident radiation is larger than the GaAs band gap [48]. This indicates that a 
carrier generation at the crystal surface is essential to observe the enhanced oxidation 
effect. Thus, the photo-generated electron-holes pairs can catalytically enhance the 
breaking of the substrate bonds and the dissociation of the adsorbate molecules [49]. 
The experimental evidence that the reaction rate per incident photon is strongly de- 
pendent on photon energy has been already observed and related to the interaction of 
hot electrons with O2 via resonant tunneling [50]. Therefore, the catalytic effect of 
electron-hole pairs photo-generated at the surface of PCNs produces a large and fine 
controlled tuning of the resonances without drastically increasing the sample temper- 
ature. The laser assisted oxidation of the sample determines a local modification of 
the nanocavity dielectric environment, reducing the effective membrane thickness and 
increasing the effective photonic pores diameter. Furthermore, due to photo-enhanced 
reactivity of GaAs surface, green light laser is much more efficient in inducing the 
oxidation with respect to light wavelengths closer to the semiconductor band gap en- 
ergy. 

Finally and remarkably, post-growth local control of the photonic modes with a 
low thermal budget that avoids annealing and/or damage of the embedded light 
sources, such as quantum wells or quantum dots, can be exploited also for changing 
on demand the coupling strength between emitters and optical modes or between in- 
teracting optical modes in systems consisting on coupled PCNs. 


5.2.3 Antibonding ground state 


We also take advantage of the laser-assisted oxidation to directly measure the 
character of the weak coupled modes in K-aligned photonic molecule, thus experi- 
mentally demonstrating the existence of a photonic molecule antibonding ground 
state. This configuration, as theoretically predicted in Sec. (5.2.1) in the zero detuning 
case, provides an antisymmetric lower energy state. First of all, we need to compen- 
sate the fabrication induced disorder, which can be relevant in weak coupled modes. 
Therefore, we act on the detuning without significantly modify the coupling strength. 
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Since the oxidation produces a blue-shift of the photonic modes, we expose the laser 
beam on the PCN in which the ground state molecular mode (P1) is initially more 
confined. The resolved confocal PL is used to image with low spatial resolution both 
coupled modes, thus roughly monitoring the spatial distributions. As reported in the 
inset of Fig. (5.11 a), we expose the right-side of a K-aligned photonic molecules at 
steps of the order of half an hour with Ar at power of 1 mW. As a function of the 
oxidation induced detuning the mode splitting (of 1.2 nm) undergoes an anticrossing 
behavior, see Fig. (5.11 a). The initial splitting is of the order of 2.1 nm and reaches 
a minimum value of 1.2 nm after 150 minutes of laser exposure. When the condition 
of minimum 2 is satisfied, a far-field angular emission pattern of both coupled modes 
is performed by means of modified confocal experimental setup developed in 
ref. [21]. This procedure allow us to retrieve constructive (destructive) interference, 
fingerprint of the bonding (antibonding) character. In particular, in the experimental 
far-field patterns of Fig. (5.11 b and d) we observe for the lower energy mode (P1) a 
destructive interference along the wavevector axis perpendicular to the alignment di- 
rection, while for the higher energy mode (P2) a vanishing signal is found along the 
same wavevector axis. These results are in agreement with the numerical FDTD cal- 
culations of the angular pattern distributions for both modes, reported in Fig. (5.11 c 
and e). This is a further signature of the antibonding ground state in K-aligned pho- 
tonic crystal molecule. 


5.2.4 Coupling strength control 


The post-fabrication control of evanescent tunnelling in a D2 based photonic crys- 
tal molecules K-aligned is demonstrated through the combination of local and selec- 
tive tuning techniques, such as infiltration of water, micro-evaporation, and laser in- 
duced non thermal oxidation [26, 42]. The value of the photonic coupling is directly 
measured by the molecular mode splitting at the anticrossing point. In fact, we are 
able to either increase or decrease the detuning (A) and the photonic coupling (g), 
independently, by inducing a local asymmetric or symmetric perturbation with respect 
to the molecule environment, respectively. From Eq. (5.2), neglecting the losses, the 
mode splitting (2) between the molecular modes results: 


Q= 4g? + 42 (5.10) 


The coupling strength described by Eq. (5.7) depends on the overlap between the 
single atomic-like modes localized in each resonator. Being directly proportional to 
(5.9), it is therefore sensitive to the fabrication-induced disorder present in the central 
region of the molecule. This would limit the possibility of developing extended array 
of identically coupled photonic molecules, unless a post-fabrication control of g is 
implemented. In order to engineer g parameter, we make the most of the non-thermal 
laser assisted oxidation presented in Sec. (5.2.2), by illuminating the sample with Ar 
laser (514 nm, P = 2.5 mW) coupled to a SNOM dielectric probe, that guarantees a 
high degree of localization [51]. The laser high-power exposure is performed by keep- 
ing the tip in the middle of the photonic molecule in order to achieve a symmetric 
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oxidation that prevents mode detuning. After each laser exposition step the PL signal 
of photonic molecule is investigated by performing SNOM scan employing low exci- 
tation Ar laser power (0.1 mW). The topography detected by means of SNOM meas- 
urements allows us to precisely induce the local laser exposure in the center of the 
photonic molecule. In addition, we choose a photonic molecule with a very small 
splitting between P1 and P2, indicating that the disorder induced detuning is quite 
small. Furthermore, to assure a negligible detuning between the two PCNs we focus 
on the separation between the strongly coupled modes P3 and P4. The main effect of 
the laser exposure is to blueshift all the photonic modes, as shown for the case of P3 
mode reported in Fig. (5.12 a) as a function of exposure time. More interestingly, we 
observe that 2 undergoes a large and monotonic reduction from (14.7 + 0.2) nm to 
(10.2 + 0.2) nm, as shown in Fig. (5.12 b). Since P3 and P4 arise from the strong 
coupling regime, the observed splitting reduction is mainly attributed to an effective 
mode interaction variation rather than to an induced detuning. This assumption is con- 
firmed also by the overlap map between the two strongly coupled modes P3-P4 re- 
ported in Fig. (5.12 e) and calculated from the difference between the FDTD calcu- 
lated electric LDOS of P3 and P4 shown in Fig. (5.12 c-d), respectively. Therefore, 
the experimental N reduction demonstrates that the laser induced oxidation is a pow- 
erful tool for adjusting the photonic tunneling rate; in this particular case for decreas- 
ing g about 30% of its initial value. 

The laser induced oxidation of GaAs, due to the smaller index of refraction of the 
oxide, is equivalent to an effective reduction of the GaAs material that provides a 
blueshift of all the resonant modes. Therefore, as pointed out in Sec. (5.2.2), this pro- 
cess could imply a reduction of the effective slab thickness and/or an enlargement of 
the effective pore sizes. In order to understand the physical origin of the observed g 
reduction, we perform three-dimensional FDTD calculations either by reducing the 
slab thickness or by enlarging the diameter size of the eleven central pores of the K- 
aligned photonic molecule, as reported in Fig. (5.13 a-d). In both cases, all the pho- 
tonic modes experience a blue-shift, as shown for the notable case of P3 in Fig. (5.13 a 
and c). However, a difference emerges in investigating the mode splitting reported in 
Fig. (5.13 b and d). On one hand, the slab reduction is almost ineffective in changing 
2, on the other hand a reduction of the calculated 2 comparable with the experimental 
data is obtained when enlarging the central eleven pore diameter sizes. A decrease of 
(4.5 + 0.4) nm, which matches the experimental value, is found. On the contrary, the 
only enlargement of the pore diameter sizes cannot explain the P3 mode experimental 
large blue-shift of (28.0 + 0.1) nm. So, we need to take into account also the slab 
thickness reduction, whose main effect, however, is only to blue-shift the resonant 
modes. Therefore, the laser induced oxidation is effective in decreasing the photon 
coupling as a consequence of the filling fraction reduction in between the two adjacent 
PCNs. Moreover, we also expect that the refractive index increase in the middle of 
the photonic molecule would result in increasing the coupling strength g. 

In order to confirm this prediction, FDTD calculations are implemented assuming 
that the central pores of the molecule are completely filled by a medium with n = 
1.33 (the same refractive index of water). We perform different simulations by infil- 
trating 1, 5, and 11 central pores of the K-aligned molecule, respectively. The more 
the holes are infiltrated, the more the resonant modes are red-shifted, as highlighted 
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for P3 in Fig. (5.14 a). The calculated Q is reported as a function of the number of 
infiltrated pores in Fig. (5.14 b), exhibiting almost the same value in all the three 
cases. Therefore, g can be increased by local infiltration, but the main contribution is 
given by filling the single central pore, as highlighted by Fig. (5.12 e). 

We therefore perform a micro-infiltration experiment with a diluted solution of 
water and rhodamine, which helps in removing air bubbles inside the pipette (as can 
be seen by optical inspection). Heuristically, we observed that this solution strongly 
improves in filling efficiently the photonic pores. Micro-infiltration experiments are 
more difficult to control since the observed increase of N can also be produced by a 
nonsymmetric infiltration, with consequent increase of 4. So, we need to directly 
measure g in an anticrossing experiment, by acting on A in order to reach the zero- 
detuning condition. Only in this condition, from Eq. (5.10), R = 2g holds, and thus a 
direct measure of the coupling strength is obtainable. For this purpose, we implement 
a detailed experimental procedure schematized in Fig. (5.14 c) and employing the 
confocal PL oxidation-detection setup, since we do not require sub-wavelength reso- 
lution for the oxidation process. After each infiltration process, or laser exposure step, 
the PL signal of photonic molecule is investigated by performing a scan with low 
power excitation Ar laser (0.1 mW). In the first step (1) we infiltrate the pores in the 
external region adjacent to the right-side PCN. Then with the confocal setup we ex- 
pose the infiltrated region with the Ar laser at 2.5 mW, forcing the water evaporation 
and inducing local photo-oxidation. This allows to measure the anticrossing curve and 
thus to extract the initial g value. In the second step (2) we perform an infiltration in 
the central region of the photonic molecule in order to change g. At the final step (3), 
we induce a second non-symmetric oxidation on the left side of the photonic molecule 
in order to determine by a second anticrossing curve the final value of g. Experimen- 
tally, the first infiltration step is implemented on the right side of the photonic mole- 
cule because the higher-energy weak-coupled mode P2 is initially found to be local- 
ized on the right cavity, as retrieved by low-power confocal PL spatial distribution 
(not shown here). After the right-side infiltration, the weak-coupled mode localized 
on the right cavity becomes red-shifted by 4 nm, while the left cavity one remains 
almost unchanged, confirming that the micro-infiltration is local and selective. 

Reporting the mode splitting between P3 and P4 (2) as a function of the laser 
exposure time in Fig. (5.14 d) we find an increase of N after the right side micro- 
infiltration that has to be ascribed to the induced detuning. Then, we force the water 
evaporation and at the same time induce oxidation by illuminating the infiltrated right- 
side region (1). By doing this, we still act on a single PCN therefore on A. As a func- 
tion of the exposure time during step (1) we observe in the blue dots of Fig. (5.14 d) 
an initial reduction of 2. Then, once the detuning is perfectly compensated, 2 starts 
to increase, thus providing an anticrossing behavior. We can precisely extract the 
value of the coupling strength by measuring the minimum 202, for an exposure time of 
about 130 min, that results in g = (6.6 + 0.1) nm. 

Now, in such a configuration, we infiltrate the central region of the molecule per- 
forming step (2) in Fig. (5.14 c) in order to increase the coupling strength. After the 
central infiltration 2 increases from 14.1 nm to 17.4 nm, see the comparison between 
the last blue dot and the first black dot in Fig. (5.14 d). However, such enlargement 
cannot be simply ascribed to a coupling increase, since the infiltration could be not 
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perfectly symmetric. Therefore, g value after the infiltration (2) is determined by 
means of a second anticrossing curve obtained through an oxidation of the left-side 
region of the photonic molecule, see step (3) in Fig. (5.14 c). The results are reported 
by the black dots in Fig. (5.14 d), where the oxidation time is kept the same as meas- 
ured by the beginning of the three steps. The experimental minimum splitting after 
the central pores infiltration gives g = (8.5 + 0.1) nm, with an increase of 28% with 
respect to the value measured before the central dielectric modification. The values of 
the photonic coupling, both before and after the infiltration (2) are in good agreement 
with the FDTD calculations reported in Fig. (5.14 d). 

The presented combination of post-fabrication methods allows us to control with 
excellent accuracy the coupling strength in photonic molecules. The deterministic 
control of the photon evanescent tunneling rate can be extended to any photonic crys- 
tal molecule geometry, thus making possible the realization of array of photonic cav- 
ities with on demand tunneling rate between each pair of coupled resonators. Finally, 
this peculiar engineering of photonic crystal molecules would open the road to possi- 
ble progress in the exploitation of coherent interference between coupled optical res- 
onators both for quantum information processing and optical communication. 


5.3 Array of photonic molecules 


Photonic molecules permit the coupling of the identical modes of the bare cavities 
giving rise to spatially delocalized and spectrally separated modes. Nowadays, the 
nanotechnology research interest is deeply focused on devices, which are formed by 
an arrangement of nano-resonators in a long array, either one-dimensional or two- 
dimensional, since these systems represent a very promising platform for manipulat- 
ing and controlling the light propagation and the modal dispersion at the nanoscale 
[52, 53]. 

Photonic molecules formed by an array of nominally identical cavities, like cou- 
pled resonator optical waveguide (CROW), exhibit spatially delocalized optical 
modes and spectral minibands, which describe the light “hopping” between adjacent 
resonators [23, 54]. The CROW systems are therefore highly promising elements for 
building large scale photonic integrated circuits and exploit “slow-light” behavior. 
However, despite the large interest in the field, very little is known on the role of 
different neighboring interaction between PCNs array molecules. 

Here, we approach the photonic molecule of three coupled PCNs in two configu- 
ration geometries, leading to different behaviors in terms of mode splitting and spatial 
distribution, as the inter-cavity alignment changes. We demonstrate by near field op- 
tics that the photonic hopping between PCNs can be driven by the interplay between 
the nearest-neighbour and the next-nearest-neighbour interactions. When the nearest- 
neighbour coupling coefficient is the dominant term, the photonic molecular modes 
behave very similarly to quantum system, such as coupled quantum wells, concerning 
the mode envelope distributions [55]. On the contrary, by designing a photonic mol- 
ecule where the role of the next-nearest-neighbour coupling is not negligible, we ob- 
tain a non-intuitive picture for the mode splitting and spatial distributions. 
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For sake of simplicity, in the following analysis we focus only on the coupling 
between the lower energy mode of each single D2 nanocavity (M1). The theoretical 
investigation of the differences between M-alignment and K-alignment is performed 
by two-dimensional plane wave expansion calculations, by modelling an infinite array 
of D2 PCNs, whose schematics are reported in the insets of Fig. (5.15 a-b), respec- 
tively. We obtain the energy band diagram projected along the wavevector component 
parallel to the coupling axis, within the first Brillouin zone of the photonic crystal for 
both the array alignments, as reported in Fig. (5.15 a-b). The finite width of both dis- 
persion curves demonstrates the effective photonic hopping between the nanocavities. 
The expected interaction difference between the two kinds of alignments is confirmed 
by modal dispersion analysis, which shows clear minibands for both configurations 
but with different trends. For the M-alignment, a monotonic and large increase of the 
energy of the photonic guided mode as a function of k, (wavevector component par- 
allel to the coupling axis) is found, as reported in Fig. (5.15 a). This pattern corre- 
sponds to a photonic guided mode with positive group velocity. On the contrary, for 
the alignment along the K-axis a small and non-monotonic variation of the array pho- 
tonic guided mode energy is found , as reported in Fig. (5.15 b). This is a signature of 
anomalous light propagation that shows a group velocity sign change, from negative 
to positive, along the modal dispersion. In the tight binding picture for a photonic 
waveguide consisting of identical coupled resonators, it is possible to estimate the 
amplitude of each neighbour coupling term for both array alignments. In fact, the 
mode interaction in CROW is the optical analogue of the tight-binding limit in con- 
densed-matter physics, in which the overlap of atomic wave functions is large enough 
that corrections to the picture of isolated atoms are required, yet, at the same time, it 
is not large enough to render the atomic description completely irrelevant [56]. There- 
fore, the calculated minibands can be well reproduced by [24]: 


hw(k) = hwo — Yim=1 Jmcos(mkD) (5.11) 


where k represents the wavevector component considered, hwy is the energy of 
the isolated single cavity, D is the spatial separation between adjacent resonators and 
the parameters Jm are the nearest-neighbour (m = 1), the next-nearest-neighbour 
(m = 2) and successive neighbouring coupling terms. Therefore, J, corresponds to 
the nearest-neighbour interaction defined as g in Sec. (5.1) and Sec. (5.2). Each cou- 
pling term is expressed by the generalization of (5.7) [24]: 


= J Ac(r)Eo(r)-Eo(r+mD)dr 
Jm = hoo Tesi ier Emer as 

where D is represented by 3V3a9 for the M-alignment and by 3a% for the K- 
alignment. By fitting the calculated energy dispersion with Eq. (5.11), where only the 
first three Jm terms are considered, we find the coupling coefficient values reported 
in Table (5.1). In the M-alignment the modal dispersion of the miniband is governed 
by the coupling between adjacent PCNs (J, /a = —12.8 meV), while the amplitudes 
of the higher order contributions J, and J3 are negligible. The case of K-aligned mol- 
ecule is quite different. The values of J, and Jz are of the same order of magnitude 
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(J1/a = 1.16 meV, J,/a = 1.55 meV), while J} is negligible. Therefore, in the K- 
aligned molecule the approximation of nearest-neighbour interaction fails and in any 
finite array a complex mode splitting and non-intuitive delocalization of the resonant 
fields is expected. Moreover, looking at the /, signs, both the coupling terms in K- 
alignment show the opposite sign with respect to J, of the M-aligned array. This sign 
inversion is in agreement with the behavior of the two coupled PCNs molecule dis- 
cussed in Sec. (5.1.2). Therefore, the negative (positive) coupling strength gives rise 
to a bonding (antibonding) ground state in M-coupled (K-coupled) arrays of any 
length. 

In order to validate the properties extracted from the array modal dispersion anal- 
ysis, we experimentally investigate the photonic molecule array composed by three 
D2 nanocavities aligned either along the M or the K-axis. The SEM image of an M- 
aligned array with a = 311 nm is reported in the inset of Fig. (5.16 a). The typical 
near-field PL spectrum, spatially averaged over the array, is reported in Fig. (5.16 a) 
and it shows three peaks almost equally spaced by 10 nm, hereafter denominated Tn 
(n = 1,2,3 for decreasing wavelength), in agreement with the theoretical FDTD cal- 
culations shown in Fig. (5.16 c). The spectral offset between experiment and theory 
has to be ascribed to the fabrication-induced disorder (not considered in the simula- 
tions) and to the slab refractive index, which is evaluated as GaAs refractive index at 
the wavelength of interest neglecting any possible oxidation or dispersion of the ma- 
terial. The summary of Tn peak positions, evaluated by fitting the spectrum with a 
three Lorentzian peaks function, is reported in Fig. (5.16 b) for different samples, 
which possess nominal identical design and so differ only for the fabrication induced 
disorder. A spread of the order of few nm in the peak positions is clear, but the sepa- 
rations between T1-T2 and T2-T3 remain almost constant for every sample. This 
means that disorder has a limited impact on the mode coupling, since for the M-align- 
ment the interaction strength is larger than the detuning. 

The case of the K-aligned array molecule is more complex. A corresponding typ- 
ical SEM image is reported in the inset of Fig. (5.16 d). The summary of the Tn peak 
positions for nominally identical samples are shown in Fig. (5.16 e). The overall re- 
sulting splitting between T1-T3 is much smaller with respect to the M-alignment, de- 
noting a weaker mode coupling. We observe a large variation of the peak splitting as 
shown in Fig. (5.16 e) indicating that the disorder induces a mode detuning, which 
can be comparable with the mode coupling. The structure where the role of disorder 
is less pronounced is likely the array#1, whose PL spectrum is reported in 
Fig. (5.16 d). It shows an almost degeneration of T1 and T2 modes. In fact, we find 
that array#1 corresponds to a system of almost identical cavities, as inferred by the 
comparison between the experimental splittings of Fig. (5.16 d) and the FDTD calcu- 
lated ones, reported in Fig. (5.16 f), relative to A = 0 configuration. 

In order to understand the physics underlying the mode splitting of the two kinds 
of array we develop a minimal analytic model of three coupled resonators. The be- 
havior of the M-aligned array can be reproduced by considering only the nearest- 
neighbour negative coupling —| gm |. Assuming three identical resonators with energy 
hwo, thus neglecting A, the coupled system mode energies and distributions are found 
by solving the eigenvalue problem of the M-matrix given by: 
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M=(-|gyl foo —gyl (5.13) 
0 —l9ul hwo 


whose solutions are the eigenvalues E, and the eigenvectors cy: 


E, = hay — V2|gy| c = (1,V2,1) (5.14) 
E, = hw c, = (-1,0,1) (5.15) 
Ez = hw +vV2|gy] c3= (1,-v2,1) (5.16) 


where the three components of c,, describe the envelope mode amplitude on the 
basis composed by the three single adjacent resonators. The energy splitting between 
the first and the last resonant mode (2V2 |g M |) is V2 times larger than the mode split- 
ting of a system composed of only two coupled resonators. This prediction is in good 
agreement with the comparison between the average splitting of 20 nm between T1- 
T3 found in the investigated M-coupled array and the 13 nm observed in the two 
PCNs M-aligned photonic molecule, with an identical coupling region [6]. In addition, 
the FDTD calculated splitting between T1-T3 together with the analytic splitting ex- 
pression from (5.14 and 5.16), allows us to retrieve the coupling strength, |gy| = 
5.6 meV for a three-dimensional system. This result is almost a factor of two smaller 
than J,. Such discrepancy can be explained by total suppression in two-dimensional 
calculations of the out-of-plane losses that enhance the field confinement strength, 
thus increasing the overlap integral of Eq. (5.12). The near-field spatial distribution 
of the array modes is addressed by imaging the PL intensity detected by the SNOM 
setup reported in Fig. (3.2). The comparison between the PL intensity map of the M- 
aligned array#2 and the electric field intensity map calculated by FDTD for the same 
structure is reported in Fig. (5.17) for every Tn modes. It is clear that TI and T3 modes 
are mainly localized in the central PCN, thus corresponding to the c, and c3 states, see 
(5.14 and 5.16). While the T2 mode is almost completely localized in the external 
PCNs, as expected for the c, state, see (5.15). Finally, we observe that, since the ex- 
perimental average mode splitting is of the order of 10 nm (corresponding to 
7.3 meV), the role of disorder induced detuning is negligible in the investigated M- 
aligned arrays. 

Concerning the K-aligned geometry a quite different behavior is observed. In this 
configuration the next-nearest-neighbour interaction plays a crucial role, as stated 
from the calculations of the two-dimensional infinite K-coupled array being, for 
which Jz/J, = 1.34, as reported in Table (5.1). The photonic molecule denoted as 
K-aligned array#1, is the structure where the disorder induced detuning is mostly neg- 
ligible, as inferred by the comparison between the experimental and the FDTD calcu- 
lated spectrum of Fig. (5.16 d and f), respectively. PL SNOM measurements are ex- 
ploited to probe the mode intensity distributions, as reported in Fig. (5.18). In partic- 
ular, the higher quality factor mode (T3) exhibits an intensity distribution mainly de- 
localized over the two external PCNs. This distribution is in agreement with the cal- 
culated FDTD electric field intensity map, as highlight by the comparison between 
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Fig. (5.18 e and f). Concurrently, the T1 and the T2 experimental mode distributions, 
whose peak wavelength is evaluated by a three Lorentzian functions fit of the entire 
spectrum, are delocalized over both central and left side PCNs, as reported in 
Fig. (5.18 a and c), respectively. The same results are retrieved by the FDTD calcu- 
lated electric field intensity maps relative to the T1 and the T2 modes shown in 
Fig. (5.18 b and d). With the exception that T1 and T2 modes are exchanged with 
respect to the experimental case. This behavior can be driven by a little disorder in- 
duced detuning (even much smaller than the coupling strength), that slightly shifts T1 
with respect to T2 or vice-versa, nevertheless resulting in a very small splitting. 

In order to confirm this description of the K-aligned molecule, we employ the 
model of three coupled resonators including also the contribution of the next-nearest- 
neighbour coupling dx = |dy| besides the nearest-neighbour term gx = |g,|, which 
have to be positive. The resonant energies and distributions are found by solving the 
eigenvalue problem of the K-matrix given by: 


hwo = A, Ix dx 
K= CA hw-4c Ik (5.17) 
dx Ix hwo 


where we also consider the possible mode detuning A, and A, relative to one lat- 
eral resonator and to the central resonator, respectively. Pushing further the analysis 
we are able to estimate the nearest and next-nearest neighbour coupling strength by 
firstly solving the eigenvalues problem for the K-matrix (5.17), when the condition 
A, = 0 =A, holds. Subsequently, by calculating the mode splittings of the FDTD 
spectrum reported in Fig. (5.16 f), which corresponds to 0.35 nm and 2.94 nm split- 
ting between T1-T2 and T1-T3, respectively. We therefore find that these wavelength 
splittings are exactly reproduced by the analytic model in the condition gy = 
0.69 meV, dg = 0.88 meV. The obtained values, like the ones for the M-coupling, 
are still two times smaller with respect to the two-dimensional estimation, but exhibit 
a ratio equal to 1.28 that is quite similar to the ratio J,/J,. The corresponding eigen- 
values, evaluated in wavelength units as 4, = cE% +, and eigenvectors are: 


A, = 1306.26 nm cı = (—0.71, 0, 0.71) (5.18) 
A, = 1305.92 nm c, = (0.38, —0.84, 0.38) (5.19) 
A; = 1302.98 nm c, = (0.59, 0.54, 0.59) (5.20) 


While the peak positions can be reproduced quite accurately, the discrepancy be- 
tween cz and both the experimental and FDTD calculated spatial distributions re- 
ported in Fig. (5.18 e and f), likely suggests that the predictions of the analytic model 
are seminal but not fully correct, since the real system is more complex than three 
coupled oscillators. In any real device the dielectric disorder is an unavoidable draw- 
back. A systematic analysis of the role of mode detuning is therefore needed and it 
could represent an outlook for nanophotonic devices to engineer the photon hopping 
between adjacent resonators in high density optical circuits. For this purpose, and in 
order to deeply study the role of detuning in both array configurations, we exploit the 


133 


Imaging and engineering optical localized modes at the nanoscale 


noteworthy eigenvalues prediction coming from the analytic model of M-matrix and 
K-matrix in presence of both A, and A,. 

In Fig. (5.19 a-d) we report the 2, in wavelength units as a function of central 
detuning and as a function of lateral detuning, respectively. In the calculations we use 
the parameters: hw = 953.7 meV (Ap = 1300 nm), gy = 5.6 meV and dy = 
0 meV in order to model the M-coupled array, see Fig. (5.19 a-b), and gg = 
0.69 meV and dg = 0.88 meV in order to model the K-coupled array, see Fig. (5.19 
c-d). The results for the M-aligned array show, as reported in Fig. (5.19 a), that Az 
remains constant and equal to A, as a function of the central detuning. While only the 

lower energy and higher energy states shift, exhibiting a minimum A; — A3 split- 
ting for A, = 0, fingerprint of an anticrossing behaviour. In fact, Az is not dependent 
form central detuning since the corresponding eigenvector cz is given by (5.15). 
Therefore, the envelope of this state vanishes in the middle cavity and any local mod- 
ification in the central molecule region does not affect the state energy or distribution. 
On the other hand, as reported in Fig. (5.19 b) with the same coupling terms of 
Fig. (5.19 a), all the eigenvalues change with A). In this case the minimum A, — 43 
splitting is still at zero detuning, while both 4, — A, and 4, — 43 exhibit minimum 
separation for A, + 0. By analyzing the model (5.17) that reproduces the K-coupled 
array interaction we find that, differently from the two coupled cavities system, the 
minimum splitting condition between every pair of resonant states as a function of 
detuning occurs in A; + 0 and A, # 0. In Fig. (5.19 c) is highlight the accidental de- 
generacy between A, and 4, driven by a finite central detuning. Moreover, for A, 
lower than the one that corresponds to the accidental degeneracy A, is not affected by 
the perturbation and its envelope distribution vanishes in the central cavity, while for 
increasing A, these properties are gained by 44. 

This inversion between the two lower energy modes as a function of A, can thus 
explain the inverted intensity distributions observed between the experimental and the 
calculated P1 and P2 modes, reported in Fig. (5.18 a-d). The trend evaluated in 
Fig. (5.19 d) as a function of lateral detuning is similar to the one reported in Fig. (5.19 
c), with the exception that no accidental degeneracy is observed and an anticrossing 
behavior between A, and 4, occurs. All these results highlight the quite complex de- 
tuning dependence of the relative mode splittings. 


In conclusion, we show that different PCNs array configuration geometries lead 
to different results in terms of mode splitting and spatial distribution. When both the 
nearest and next-nearest-neighbour coupling terms are relevant, a non-intuitive pic- 
ture describes the resonant modes, in addition the photon hopping between adjacent 
nanoresonators is strongly affected. These findings could open the way to exploit ex- 
otic configurations of coupled PCNs to tailor the mode spatial distribution or the group 
velocity in CROW of any length by engineering the dielectric properties of adjacent 
resonators in high-density optical circuits. 
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Figure 5.1: SNOM measurements on M-aligned photonic molecule with a = 311 nm, 
performed by SNOM PL with dielectric tip. Figures adapted from [6]. a) Average 
spectrum of the molecule highlighted in the inset with the direction of the M-axis. 
b) Spectral shift maps of P1, P2, P3 and P4, respectively. The scale bar is 300 nm. 
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Figure 5.2: SNOM PL measurements on K-aligned photonic molecule with a = 
311 nm, performed with a dielectric tip. Figures adapted from [6]. a) Average spec- 
trum of the molecule highlighted in the inset with the direction of the K-axis. b) Spec- 
tral shift maps of P1, P2, P3 and P4, respectively. The scale bar is 300 nm. 
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Figure 5.3: FDTD calculations of the K-coupled photonic molecule without detuning. 
a) Normalized spectrum averaged over the molecule reported in the inset. The P1 
(lower energy) and P2 (higher energy) peaks correspond to the antibonding (A) and 
bonding (B) state, respectively. b)-c) Spatial distribution of the real part of the electric 
field and angular distribution in the wavevector space, relative to the antibonding 
mode (P1), respectively. d)-e) Spatial distribution of the real part of the electric field 
and angular distribution in the wavevector space, relative to the bonding mode (P2), 
respectively. The white circles in d) and e) indicate the aperture angles of 30° and 60°. 
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Figure 5.4: a) Map of the quantity (5.9) evaluated by FDTD. b) Scheme of the K- 
molecule where the size of the five central pores (red) is reduced for tuning the cou- 
pling strength. c) Calculated peak position of the B and A modes as a function of the 
diameter size of the five central pores highlighted in b). The blue circle corresponds 
to degenerate condition. d) Calculated spectral broadening of the ground state (P1) 
and the excited state (P2) as a function of the diameter size of the five central pores. 
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Figure 5.5: a)-b) Resonant energy of both bonding and antibonding modes in M- 
aligned and K-aligned photonic molecules, respectively, as a function of center-to- 
center separation (in units of the hole-to-hole distance). The calculations are per- 
formed by the guided-mode-expansion method [31]. In the insets are reported the de- 
sign of the photonic molecule exhibiting a center-to-center separation of 3V3a and 
3a for the M-aligned and K-aligned photonic molecules, respectively. 
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Figure 5.6: Left panel: schematics of the K-coupled molecule with five central pores 
diameter size reduced to 173.2 nm. First row, nominal sample with n = 3.484; second 
row, reduction to n = 3.416 in the central circular region (pale yellow circle); third 
row, reduction to n = 3.348 in the central circular region (yellow circle). Right panel: 
FDTD spectra relative to the structure in the left panel. The bonding or antibonding 
character is highlighted. 
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Figure 5.7: a) PL intensity map of the mode M1 of the PCN with a = 311 nm reported 
in the inset, obtained at low excitation power (P = 0.1 mW) with the 514 nm laser. 
b) M1 spectra acquired for increasing power excitation of the diode laser (780 nm) 
from 0.25 mW to 2.7 mW. c) Map of the M1 mode spectral shift measured with high 
excitation density (P = 2.7 mW) provided by the 780 nm laser. d) Spatial profile of 
the thermal induced spectral shift evaluated from the average line cut in c). 
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Figure 5.8: a) Spatial distribution of the temperature gradient calculated by FEM on 
the top surface of a D2 nanocavity with slab thickness 320 nm, employing laser irra- 
diation at 516 nm with a power of 0.7 mW. b) Radial spatial profile extracted from a). 
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Figure 5.9: a) PL spectra of the investigated D2 nanocavity acquired at low Ar laser 
power (0.1 mW) after three subsequent exposures to high Ar laser. power (0.7 mW). 
Starting from the bottom, the second and the third spectra are collected after 120 
minutes and 200 minutes of high power laser exposition with respect to the first spec- 
trum, respectively. b) M1 mode wavelength blue-shift plotted as a function of the 
square root of high power exposition time both for Ar and diode lasers. 
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Figure 5.10: a) Topography of the PCN after a laser exposure corresponding to an 
overall shift of 13 nm for the M1 mode. The map is collected by a dielectric probe 
that acts as shear-force sensor. A protrusion is observed inside the dotted circle due to 
the laser local oxidation. b) Line profile of a) (black circles) for evaluating the height 
of the oxidized region, being approximately 20 nm. The red dashed line is a guide to 
the eyes obtained through a FFT low pass filtering. 
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Figure 5.11: a) Splitting (2) between the weakly coupled modes P1 and P2 in K- 
aligned photonic molecules, as a function of the Ar laser exposure time. The local 
oxidation is performed by means of confocal setup illumitating with 1.0 mW on the 
right side of the molecule, as shown in the inset, since the ground state mode (P1) is 
initially localized on the right-side PCN. b)-c) Experimental far-field angular distri- 
bution and FDTD calculation, respectively, of the ground state molecular mode (P1) 
in the condition 4 = 0. d)-e) Experimental far-field angular distribution and FDTD 
calculation, respectively, of the excited state molecular mode (P2) in the condition 
A = 0. The white circles in b)-e) indicate the aperture angles of 30° and 60°. 
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Figure 5.12: a)-b) Experimental P3 mode wavelength and splitting 2 between P3 and 
P4, respectively, as a function of the Ar laser exposure time on K-aligned photonic 
molecule. The schematics of the exposition is reported in the inset of a). c)-d) Electric 
LDOS calculated by FDTD at the peak wavelength of P3 and P4, respectively. e) Dif- 
ference between the maps in c) and d) to highlight the P3-P4 superposition regions. 
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Figure 5.13: FDTD calculations performed on K-aligned photonic. a)-b) P3 mode 
wavelength and £2 as a function of the GaAs membrane thickness, respectively. c)-d) 
P3 mode wavelength and 2, respectively, as a function of the diameter size of the 
eleven central pores of the K-aligned photonic molecule reported in the inset of c). 
The error bars are provided by the finite spectral resolution of the FDTD calculations. 
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Figure 5.14: a)-b) P3 mode wavelength and 22, respectively, calculated by FDTD as a 
function of the number of infiltrated pores with n = 1.33, as schematically shown in 
the inset of a). c) Schematic view of the three experimental steps performed. d) 2 is 
reported as a function of the exposure time of subsequent lateral micro-evaporation 
and micro-oxidations, as described in c). The red dot at negative time corresponds to 
the value before the right-side infiltration. The blue circles represent 2 after the right- 
side infiltration and during the right-side laser illumination (1). The black dots repre- 
sent N after the central infiltration (2) and during the left-side laser exposure (3). 
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Figure 5.15: Two-dimensional calculations of the dispersion relation (black dots) of 
infinite array of D2 nanocavities, performed with plane wave expansion, reported in 
the energy range of the single cavity fundamental mode, whose energy is highlighted 
by horizontal dashed line. a) Calculation relative to the M-coupled infinite array (sche- 
matics in the inset) in the irreducible part of the first Brillouin zone, with D = 3V3a. 
The y-axis is parallel to the TM direction of the reciprocal lattice. The red curve pro- 
vides the monotonous dispersion fit performed with Eq. (5.11). b) Calculation relative 
to the K-coupled infinite array, (schematics in the inset) in the irreducible part of the 
first Brillouin zone, with D = 3a. The x-axis is parallel to the TK direction of the 
reciprocal lattice. The red curve provides the non-monotonous dispersion fit per- 
formed with equation Eq. (5.11). The scale bar in both insets is 600 nm. 
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Figure 5.16: a) Experimental PL near-field spectrum averaged over the three M-cou- 
pled array nanocavities (see SEM image in the inset). Three sharp peaks (T1,T2,T3) 
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are observed. b) Resonant modes wavelength in nominally identical M-coupled ar- 
rays, evaluated by fitting each PL mean spectrum with a three Lorentzian peaks func- 
tion. The second case in b) corresponds to the spectrum in a). c) FDTD spectrum 
averaged over the M-coupled array, see the inset schematics. d)-f) Same analysis of 
a)-c), but concerning the K-aligned array. e) Resonant modes wavelength for nomi- 
nally identical K-coupled arrays. The first case corresponds to the spectrum shown 
in b). The scale bar is 600 nm. 
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Figure 5.17: M-coupled array: PL near-field intensity map and FDTD calculated elec- 
tric field intensity relative to T1 mode a)-b), T2 mode c)-d) and T3 mode e)-f). The 
experimental maps are relative to the array investigated in Fig. (5.16 a), while the 
FDTD maps correspond to the peak wavelength of modes reported in Fig. (5.16 c). 
The scale bar in all the maps is 600 nm. 
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Figure 5.18: K-coupled array: PL near-field intensity map and FDTD calculated elec- 
tric field intensity relative to T1 mode a)-b), T2 mode c)-d) and T3 mode e)-f). The 
experimental maps are relative to the array investigated in Fig. (5.16 d), while the 
FDTD maps correspond to the peak wavelength of modes reported in Fig. (5.16 f). 
The disorder induced detuning could be the responsible for the exchange of the ex- 
perimental intensity distribution between T1 and T2, with respect to the FDTD calcu- 
lation. The scale bar in all the maps is 600 nm. 
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Figure 5.19: a)-b) Eigenvalues of M-matrix (5.13) as 
cavity detuning, normalized to |gy|, respectively. The used parameters are: Ey = 
953.7 meV (Ao = 1300 nm), gy = —5.6 meV and dy = 0 meV, reported in wave- 
length units. c)-d) Eigenvalues of K-matrix (5.17) as a function of central and lateral 
cavity detuning, normalized to |gx|, respectively, calculated with Ey = 953.7 meV, 
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Gx = 0.69 meV and dy = 0.88 meV, reported in wavelength units. 


Jı [a/a] J [a/A] J3 [a/a] 


M-alignment 


K-alignment 


+3.9- 1074 | -0.7- 1074 


Table 5.1: Values of the coupling terms Jm, in a/A units, calculated for the two-di- 
mensional infinite array aligned along M and K-direction, respectively, and obtained 
by fitting with Eq. (5.11) the dispersion relations calculated in Fig. (5.15 a-b). 


145 


a function of central and lateral 


Imaging and engineering optical localized modes at the nanoscale 


References 


[1] J. Joannopoulos, S. Johnson, J.Winn, and R. Meade. Photonic Crystals: Molding 
the Flow of Light. Princeton University Press, 2nd edn., (2008). 

[2] S. Noda, M. Fujita, and T. Asano. “Spontaneous-emission control by photonic 
crystals and nanocavities” In: Nat. Photon. 1 (2007), pp. 449-458. 

[3] M. Bayer, T. Gutbrod, J. Reithmaier, A. Forchel, T. Reinecke, P. Knipp, A. 
Dremin, and V. Kulakovskii. “Optical Modes in Photonic Molecules” In: Phys. Rev. 
Lett. 81, 12 (1998), pp. 2582-2585. 

[4] S. Ishii, K. Nozaki, and T. Baba. “Photonic Molecules in Photonic Crystals” In: 
Jpn. J. Appl. Phys. 45, 8R (2006), p. 6108. 

[5] F. Brossard, B. Reid, C. Chan, X. Xu, J. Griffiths, D. Williams, R. Murray, and R. 
Taylor. “Confocal micro-photoluminescence mapping of coupled and detuned states 
in photonic molecules” In: Opt. Express 21, 14 (2013), pp. 16934-16945. 

[6] S. Vignolini, F. Intonti, M. Zani, F. Riboli, D. Wiersma, L. Li, L. Balet, M. 
Francardi, A. Gerardino, A. Fiore, and M. Gurioli. “Near-field imaging of coupled 
photonic-crystal microcavities” In: Appl. Phys. Lett. 94, 15 (2009), p. 1103. 

[7] C. Schmidt, A. Chipouline, T. Kasebier, E.-B. Kley, A. Tiinnermann, T. Pertsch, 
V. Shuvayev, and L. Deych. “Observation of optical coupling in microdisk resona- 
tors” In: Phys. Rev. A 80, 4 (2009), p. 3841. 

[8] A. Nakagawa, S. Ishii, and T. Baba. “Photonic molecule laser composed of Galn- 
AsP microdisks” In: Appl. Phys. Lett. 86, 4 (2005), p. 1112. 

[9] Y. Rakovich, J. Donegan, M. Gerlach, A. Bradley, T. Connolly, J. Boland, N. 
Gaponik, and A. Rogach. “Fine structure of coupled optical modes in photonic mole- 
cules”. In: Phys. Rev. A 70, 5 (2004), p. 1801. 

[10] S. Vignolini, F. Riboli, F. Intonti, D. Wiersma, L. Balet, L. Li, M. Francardi, A. 
Gerardino, A. Fiore, and M. Gurioli. “Mode hybridization in photonic crystal mole- 
cules” In: Appl. Phys. Lett. 97, 6 (2010), p. 3101. 

[11] D. Smith, H. Chang, K. Fuller, A. Rosenberger, and R. Boyd. “Coupled-resona- 
torinduced transparency” In: Phys. Rev. A 69, 6 (2004), p. 3804. 

[12] Q. Xu, S. Sandhu, M. Povinelli, J. Shakya, S. Fan, and M. Lipson. “Experimental 
Realization of an On-Chip All-Optical Analogue to Electromagnetically Induced 
Transparency” In: Phys. Rev. Lett. 96, 12 (2006), p. 3901. 

[13] S. Yanik M.F.and Fan. “Stopping Light All Optically” In: Phys. Rev. Lett. 92, 8 
(2004), p. 3901. 

[14] A. Dousse, J. Suffezynski, A. Beveratos, O. Krebs, A. Lemaitre, I. Sagnes, J. 
Bloch, P. Voisin, and P. Senellart. “Ultrabright source of entangled photon pairs” In: 
Nature 466 (2010), pp. 217-220. 

[15] A. Imamoglu, D. Awschalom, G. Burkard, D. DiVincenzo, D. Loss, M. Sherwin, 
and A. Small. “Quantum Information Processing Using Quantum Dot Spins and Cav- 
ity QED” In: Phys. Rev. Lett. 83, 20 (1999), pp. 4204-4207. 

[16] J. Vasco, P. Guimaraes, and D. Gerace. “Long-distance radiative coupling be- 
tween quantum dots in photonic crystal dimers” In: Phys. Rev. B 90, 15 (2014), p. 
155436. 


146 


Niccolò Caselli 


[17] D. Angelakis, M. Santos, and S. Bose. “Photon-blockade-induced Mott transi- 
tions and XY spin models in coupled cavity arrays” In: Phys. Rev. A 76, 3 (2007), p. 
1805. 

[18] D. Gerace, H. Tiireci, A. Imamoglu, V. Giovannetti, and R. Fazio. “The quan- 
tumoptical Josephson interferometer” In: Nature Phys. 5 (2009), pp. 281-284. 

[19] M. Hartmann, F. Brandào, and M. Plenio. “Quantum many-body phenomena in 
coupled cavity arrays” In: Laser Photon. Rev 2, 6 (2008), pp. 527-556. 

[20] A. Tomadin, V. Giovannetti, R. Fazio, D. Gerace, I. Carusotto, H. Tiireci, and A. 
Imamoglu. “Signatures of the superfluid-insulator phase transition in laser-driven dis- 
sipative nonlinear cavity arrays” In: Phys. Rev. A 81, 6 (2010), p. 061801. 

[21] F. Intonti, F. Riboli, N. Caselli, M. Abbarchi, S. Vignolini, D.Wiersma, A. 
Vinattieri, D. Gerace, L. Balet, L. Li, M. Francardi, A. Gerardino, A. Fiore, and M. 
Gurioli. “Young’s Type Interference for Probing the Mode Symmetry in Photonic 
Structures” In: Phys. Rev. Lett. 106, 14 (2011), p. 143901. 

[22] M. Brunstein, T. Karle, I. Sagnes, F. Raineri, J. Bloch, Y. Halioua, G. Beaudoin, 
L. Le Gratiet, J. Levenson, and A. Yacomotti. “Radiation patterns from coupled pho- 
tonic crystal nanocavities” In: Appl. Phys. Lett. 99, 11 (2011), p. 1101. 

[23] A. Yariv, Y. Xu, R. Lee, and A. Scherer. “Coupled-resonator optical waveguide: 
a proposal and analysis” In: Opt. Lett. 24, 11 (1999), pp. 711-713. 

[24] F.-S. Chien, J. Tu, W.-F. Hsieh, and S.-C. Cheng. “Tight-binding theory for cou- 
pled photonic crystal waveguides” In: Phys. Rev. B 75, 12 (2007), p. 125113. 

[25] A. Chalcraft, S. Lam, B. Jones, D. Szymanski, R. Oulton, A. Thijssen, M. 
Skolnick, D. Whittaker, T. Krauss, and A. Fox. “Mode structure of coupled L3 pho- 
tonic crystal cavities” In: Opt. Express 19, 6 (2011), pp. 5670-5675. 

[26] F. Intonti, N. Caselli, S. Vignolini, F. Riboli, S. Kumar, A. Rastelli, O. Schmidt, 
M. Francardi, A. Gerardino, L. Balet, L. Li, A. Fiore, and M. Gurioli. “Mode tuning 
of photonic crystal nanocavities by photoinduced non-thermal oxidation” In: Appl. 
Phys. Lett. 100, 3 (2012), p. 3116. 

[27] C.-S. Kee, H. Lim, and J. Lee. “Coupling characteristics of localized photons in 
two-dimensional photonic crystals” In: Phys. Rev. B 67, 7 (2003), p. 073103. 

[28] M. Doty, J. Climente, M. Korkusinski, M. Scheibner, A. Bracker, P. Hawrylak, 
and D. Gammon. “Antibonding Ground States in InAs Quantum-Dot Molecules” In: 
Phys. Rev. Lett. 102, 4 (2009), p. 7401. 

[29] B. Gallinet, J. Kupec, B. Witzigmann, and M.-A. Dupertuis. “Analysis of pho- 
tonic crystal defect modes by maximal symmetrization and reduction” In: J. Opt. Soc. 
Am. B 27. 7 (2010), pp. 1364-1380. 

[30] K. Atlasov, K. Karlsson, A. Rudra, D. B., and E. Kapon. “Wavelength and loss 
splitting in directly coupled photonic-crystal defect microcavities” In: Opt. Express 
20, 16 (2008), pp. 255-264. 

[31] L. Andreani and D. Gerace. “Photonic-crystal slabs with a triangular lattice of 
triangular holes investigated using a guided-mode expansion method” In: Phys. Rev. 
B 73, 23 (2006), p. 5114. 

[32] K. Hennessy, C. Hégerle, E. Hu, A. Badolato, and A. Imamoglu. “Tuning pho- 
tonic nanocavities by atomic force microscope nano-oxidation” In: Appl.Phys. Lett. 
89, 4 (2006), p. 1118. 


147 


Imaging and engineering optical localized modes at the nanoscale 


[33] C. Chen, J. Zheng, T. Gu, J. McMillan, M. Yu, G.-Q. Lo, D.-L. Kwong, and C. 
Wong. “Selective tuning of high-Q silicon photonic crystal nanocavities via laser- 
assisted local oxidation” In: Opt. Express 19, 13 (2011), pp. 12480-12489. 

[34] I. Fushman, E.Waks, D. Englund, N. Stoltz, P. Petroff, and J. Vuckovic. “Ultra- 
fast nonlinear optical tuning of photonic crystal cavities” In: Appl.Phys. Lett. 90, 9 
(2007), p. 1118. 

[35] M. Loncar, T. Yoshie, A. Scherer, P. Gogna, and Y. Qiu. “Low-threshold pho- 
tonic crystal laser” In: Appl.Phys. Lett. 81, 15 (2002), pp. 2680-2682. 

[36] T. Tanabe, M. Notomi, S. Mitsugi, A. Shinya, and E. Kuramochi. “All-optical 
switches on a silicon chip realized using photonic crystal nanocavities” In: Appl. 
Phys. Lett. 87, 15 (2005), p. 1112. 

[37] K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. Atatiire, S. Gulde, S. Fält, 
E. Hu, and A. Imamoglu. “Quantum nature of a strongly coupled single quantum dot- 
cavity system” In: Nature 445 (2007), pp. 896-899. 

[38] F. Intonti, S. Vignolini, F. Riboli, A. Vinattieri, D.Wiersma, M. Colocci, L. Balet, 
C. Monat, C. Zinoni, L. Li, R. Houdré, M. Francardi, A. Gerardino, A. Fiore, and M. 
Gurioli. “Spectral tuning and near-field imaging of photonic crystal microcavities” In: 
Phys. Rev. B 78, 4 (2008), p. 041401. 

[39] A. Faraon and J. Vuckovic. “Local temperature control of photonic crystal de- 
vices via micron-scale electrical heaters” In: Appl.Phys. Lett. 95, 4 (2009), p. 3102. 
[40] K. Hennessy, A. Badolato, A. Tamboli, P. Petroff, E. Hu, M. Atatiire, J. Dreiser, 
and A. Imamoglu. “Tuning photonic crystal nanocavity modes by wet chemical digital 
etching” In: Appl. Phys. Lett. 87, 2 (2005), p. 1108. 

[41] X. Yang, C. Chen, C. Husko, and C. Wong. “Digital resonance tuning of high- 
Q/Vn silicon photonic crystal nanocavities by atomic layer deposition” In: Appl. 
Phys. Lett. 91, 16 (2007), p. 1114. 

[42] F. Intonti, S. Vignolini, F. Riboli, M. Zani, D. Wiersma, L. Balet, L. Li, M. 
Francardi, A. Gerardino, A. Fiore, and M. Gurioli. “Tuning of photonic crystal cavi- 
ties by controlled removal of locally infiltrated water” In: Appl.Phys. Lett. 95, 17 
(2009), p. 3112. 

[43] L. Midolo, P. van Veldhoven, M. Déndar, R. Nétzel, and A. Fiore. “Electrome- 
chanical wavelength tuning of double-membrane photonic crystal cavities” In: 
Appl.Phys. Lett. 98, 21 (2011), p. 1120. 

[44] H. Lee, S. Kiravittaya, S. Kumar, J. Plumhof, L. Balet, L. Li, M. Francardi, A. 
Gerardino, A. Fiore, A. Rastelli, and O. Schmidt. “Local tuning of photonic crystal 
nanocavity modes by laser-assisted oxidation” In: Appl.Phys. Lett. 95, 19 (2009), p. 
1109. 

[45] S. Vignolini, F. Intonti, L. Balet, M. Zani, F. Riboli, A. Vinattieri, D. Wiersma, 
M. Colocci, L. Li, M. Francardi, A. Gerardino, A. Fiore, and M. Gurioli. “Nonlinear 
optical tuning of photonic crystal microcavities by near-field probe” In: Appl.Phys. 
Lett. 93, 2 (2008), p. 3124. 

[46] A. de Rossi, M. Lauritano, S. Combrié, Q. Tran, and C. Husko. “Interplay of 
plasma-induced and fast thermal nonlinearities in a GaAs-based photonic crystal 
nanocavity” In: Phys. Rev. A 79, 4 (2009), p. 3818. 


148 


Niccolò Caselli 


[47] C.-H. Tsai, S.-R. Jian, and H.-C. Wen. “Tip-induced local anodic oxidation on 
p-GaAs surface with non-contact atomic force microscopy” In: Applied Surface Sci- 
ence 254, 5 (2007), pp. 1357-1362. 

[48] T. Suzuki and M. Ogawa. “Degradation of photoluminescence intensity caused 
by excitation-enhanced oxidation of GaAs surfaces” In: Appl.Phys. Lett. 31, 7 (1977), 
pp. 473-475. 

[49] V. Bermudez. “Photoenhanced oxidation of gallium arsenide” In: J. Appl. Phys. 
54, 11 (1983), pp. 6795-6798. 

[50] Y. Chen, J. Seo, F. Stepniak, and J. Weaver. “Visible-light-induced oxidation for 
O2 on GaAs(110): The role of hot electrons” In: J. Chem. Phys 95, 11 (1991), pp. 
8442-8448. 

[51] F. Intonti, V. Matarazzo, A. Nasir, O. Makarovsky, R. Campion, A. Patanè, S. 
Kumar, A. Rastelli, O. Schmidt, and M. Gurioli. “Nano-sized light emitting diodes by 
near field laser exposure” In: Appl.Phys. Lett. 98, 18 (2011), p. 3102. 

[52] N. Matsuda, H. Takesue, K. Shimizu, Y. Tokura, E. Kuramochi, and M. Notomi. 
“Slow light enhanced correlated photon pair generation in photonic-crystal coupled- 
resonator optical waveguides” In: Opt. Express 21, 7 (2013), pp. 8596-8604. 

[53] N. Matsuda, E. Kuramochi, H. Takesue, and M. Notomi. “Dispersion and light 
transport characteristics of large-scale photonic-crystal coupled nanocavity arrays” In: 
Opt. Lett. 39, 8 (2014), pp. 2290-2293. 

[54] J. Ma, L. Martinez, S. Fan, and M. Povinelli. “Tight-binding calculation of radi- 
ation loss in photonic crystal CROW” In: Opt. Express 21, 2 (2013), pp. 2463-2473. 
[55] D. Smith and C. Mailhiot. “Theory of semiconductor superlattice electronic 
structure” In: Rev. Mod. Phys. 62, 1 (1990), pp. 173-234. 

[56] N. Ashcroft and N. Mermin. Solid State Physics. Harcourt Inc., Orlando, FL, 
USA, (1976). 


149 


Chapter 6 


Engineering light confinement in disordered media 


Disordered photonic materials can diffuse and localize light through random mul- 
tiple scattering, offering opportunities to study mesoscopic phenomena, control light- 
matter interactions, and provide new strategies for photonic applications. Light 
transport in such media is governed by photonic modes characterized by resonances 
with finite spectral width and spatial extent. Considerable steps have been made to- 
wards control over the transport using wavefront shaping techniques. The selective 
engineering of individual modes, however, has been addressed only theoretically. 
Here, we experimentally demonstrate the possibility to engineer the confinement and 
the mutual interaction of modes in a two-dimensional disordered photonic structure. 
The strong light confinement is achieved at the fabrication stage by an optimization 
of the structure and an accurate and local tuning of the mode resonance frequencies is 
achieved via post-fabrication processes. To show the versatility of our technique, we 
selectively control the detuning between overlapping localized modes and observe 
both frequency crossing and anti-crossing behaviors, thereby paving the way for the 
creation of open transmission channels in strongly scattering media. 


6.1 Disorder driven light confinement 


The ability to mold the flow of light at the wavelength scale has been largely 
investigated in photonic-crystal-based devices, in which the light propagation is 
driven by interference between multiply Bragg scattered waves giving rise to Bloch 
modes, which can be engineered by varying the structural parameters of the material 
[1]. In disordered media, both the direction and phase of the propagating waves are 
randomized in a complex manner, making any attempt to control light propagation 
particularly challenging. 

Disordered media are at present investigated in several contexts, ranging from the 
study of collective multiple scattering phenomena to cavity random lasing and to the 
possibility to provide efficient solutions in renewable energy, imaging and spectros- 
copy-based applications [2-7]. Transport in such systems can be described in terms 
of photonic modes, or quasimodes, which exhibit characteristic spatial profiles and 
spectra [8]. In diffusive systems, these modes are spatially and spectrally overlapping 
whereas in the regime of Anderson localization, they become spatially and spectrally 
isolated [9]. Unlike Bloch modes in periodic systems, the precise formation of pho- 
tonic modes in a single realization of the disorder is unpredictable. Control over light 
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transport can be obtained by shaping the incident wave to excite only a specific part 
of the modes available in a given system [7,10]. To fully exploit the potential of dis- 
ordered systems, however, mode control is needed. It was shown theoretically that 
isolated modes could be selectively tuned and possibly coupled to each other by a 
local fine modification of the dielectric structure [11, 12]. 

We experimentally prove the ability to control the spectral properties of an indi- 
vidual photonic mode in a two-dimensional disordered photonic structure. The disor- 
dered pattern of air holes is fabricated on the same suspended GaAs membrane with 
InAs QDs embedded, already presented in Sec. (3.2) for the realization of PCNs. The 
difference relies on the disorder distribution of circular holes, as reported in the SEM 
images of Fig. (6.1 a-b). We design different disordered structures with holes filling 
fraction (f) ranging from 0.13 to 0.35 and hole diameters size (hg) ranging from 
195 nm to 270 nm. Fabrications induced deviations from the nominal structural pa- 
rameters values are measured with a scanning electron microscope and are found to 
be within 5% for hg and 10% for f. To avoid merging between adjacent holes during 
the growth process (proximity effects), we imposed a minimum distance equal to 
1.3 ha between the centers of the nearest neighboring holes. 

The PL of the sample is excited and detected by means of the SNOM setup in 
illumination-collection geometry presented in Sec. (3.1), equipped with dielectric 
probes. As excitation source we use a diode laser (780 nm) delivering to the near- 
field probe 60 uW power. The combination of spectrometer and photodiodes array 
allows us to investigate during a single SNOM scan the spectral range between 
1150 nm and 1380 nm wih a resolution of 0.45 nm. Figure (6.2) shows typical 
SNOM PL maps evaluated at a fixed emission wavelength for samples with f = 0.35 
and f = 0.13, respectively, whose intensity values are normalized by the average in- 
tensity of the sample. Several localized modes are identified by the bright spots in the 
intensity distributions and by the peaks in the wavelength spectrum. For instance, two 
spectra corresponding to the maximum intensity observed in Fig. (6.2 a-b) are re- 
ported in Fig. (6.2 c-d). The main difference between the two spectra is the peak am- 
plitude value. In fact, the PL spectrum of the sample with f = 0.35 supports a peak 
more than ten times more intense than the average intensity value of the entire sample; 
whereas the PL spectrum of the sample with f = 0.13 exhibits a peak with an inten- 
sity 3.6 times the average value. Although such intensity distributions and spectra 
provide useful local information about two single peaks, a statistical analysis of all 
the investigated wavelength and positions is required to understand the mesoscopic 
properties of every sample, and in particular the influence of the filling fraction and 
the hole diameter size on the localized modes. 

Figures (6.3 a-b) shows the probability density function (P) of PL normalized in- 
tensity (Ip, /(IpL)) in linear and semi-logarithmic scale, respectively. Each distribu- 
tion is evaluated taking into account all the collected data during the SNOM scan. By 
increasing f the P distributions become more asymmetric, see Fig. (6.3 a), and heavy- 
tailed, see Fig. (6.3 b), with a notably higher probability to observe high-intensity 
peaks. For instance, peaks with a normalized intensity about 4 are roughly 100 times 
more probable in samples with f = 0.35 than in those with f = 0.20. The onset of 
the occurrence of rare, yet very bright peaks can be quantified by the normalized var- 
iance 02/1? of P, where o and u = (Ip) are the variance and the mean intensity of 
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the distribution, respectively. Figure (6.3 c) shows that 02/1? increases by at least a 
factor of five on increasing the filling fraction from f = 0.13 to f = 0.35. The influ- 
ence of hg on the normalized variance is less evident, as shown in Fig. (6.3 d). 

The PL fluctuations of the integrated quantum sources are primarily driven by 
LDOS fluctuations, which in turn are related to the light transport properties of the 
system as well as to the near-field interaction between sources and scatterers [13-16]. 
In our experimental configuration the interplay between the involved processes, such 
as the excitation efficiency of QDs, the collection efficiency of the SNOM tip and the 
light transport in the disordered medium, does not allow us to quantitatively relate the 
slope of the P distributions tails to the dimensionless conductance, which is the uni- 
versal scaling parameter that quantifies the wave transport in a random medium [17]. 
Nevertheless, the appearance of longer tails for samples with increasing f is a signa- 
ture of the enhanced confinement of the photonic quasi-modes. On the other hand, the 
weak dependence of the fluctuations on f = 0.35 shows that the investigated disor- 
dered systems are structurally independent of possible dispersion of the scatterer size, 
as highlighted in Fig. (6.3 d). To retrieve the average spatial extension of the photonic 
modes, we exploit the quality factor decrease in proximity of the sample borders due 
to in-plane leakages. In fact, the more confined modes exhibit a Q reduction at shorter 
distances from it. The average mode decay length is then evaluated by investigating 
the mode spectral width dependence on the distance to the border [18]. Figures (6.4 
a-b) show the FWHM of the averaged wavelength autocorrelation function of the 
spectra, collected in spatial strips parallel to the sample border as a function of dis- 
tance to the sample border. The average wavelength autocorrelation FWHM is calcu- 
lated within a vertical strip (200 nm width 9 um height) parallel to the sample border 
and contains 45 near-field spectra. The averageFWHMdecreases very rapidly for 
samples with f = 0.35 (hg = 200 nm), reaching a constant value at a distance of 3.4 
um away from the border, whereas for sam ples with f = 0.13 (hg = 200 nm) the 
decay is much slower, reaching a constant value at a distance of approximately 10 — 
12 um. The slight increase of the FWHM in the bulk of the sample with f = 0.35 is 
probably due to a small increase in the hole diameter sizes caused by a proximity 
effect of the electron beam during the writing stage. Figure (6.4 c) shows the decay 
length, calculated from exponential fits of the FWHM decrease, as a function of the 
density of scatterers. The f = 0.35 samples support modes with a decay length of 
approximately È = 0.8 um with a ratio L/2È = 15, where L = 25 um is the sample 
size; whereas for f = 0.13 the decay length is approximately È = 3.5 um with a ratio 
L/2E = 3.5. 

The average FWHM of frequency autocorrelation in the middle of the samples is 
related to the spectral width of localized modes, which provide the average intrinsic 
quality factor Q = 440 + 80 for f = 0.35. This value is bounded by the out-of-plane 
losses, which are roughly four times larger than the in-plane losses (inherent to the 
finite size of the sample), as we verified by a comparison between two- and three- 
dimensional FDTD calculations (not shown here). 
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6.2 Engineering disordered photonic modes 


Having identified a proper set of design parameters to create disordered structures 
that support strongly confined photonic modes, we implement a post-processing step 
for achieving a fine-tuning of the resonant modes wavelength. The technique is based 
on the local modification of the dielectric environment in the proximity of the maxi- 
mum intensity of the modes. In a sample found to be in the regime of strongest con- 
finement, that is with f = 0.35 and hg = 200 nm, we demonstrate the ability to per- 
form a gentle and reversible spectral tuning on a single localized mode. We also engi- 
neer the spectral splittings of selected pairs of modes, which show a crossing or an 
anti-crossing feature, depending on the mutual spatial overlap and therefore on the 
mode interaction strength, as discussed for coupled PCNs in Sec. (5.1.1). 

A reversible and smooth engineering is achieved by exploiting the perturbation 
induced by the dielectric SNOM tip, which modifies the local dielectric environment 
where multiple light scattering takes place, inducing a red-shift of the modes propor- 
tional to the intensity of the confined electric field, as deeply discussed in Sec. (2.5). 
Figure (6.5 a) shows the intensity image of an isolated mode, whereas Fig. (6.5 b) 
shows the normalized intensity spectra taken close to the PL maximum intensity of 
the mode (point A) and on the PL low intensity tail (point B). The spectral shift map 
is reported as three-dimensional function in the inset of Fig. (6.5 b) and as the ordinary 
two-dimensional distribution in Fig. (6.5 c). Therefore, the localized electric field in- 
tensity distribution clearly describes the speckle patterns typical of disorder induced 
localization [18-20]. Moreover, the amplitude of the observed maximum spectral shift 
is 0.43 nm, which represents a small fraction of the mode linewidth (about 3 nm), but 
demonstrates that the dielectric perturbation in the near-field is able to slightly modify 
the resonant frequency of a disordered resonant mode. The spectral shift maximum 
value also allows us to give an estimation of the modal volume V of the localized 
mode. The polarizability of the tip employed, and therefore a direct calibration-meas- 
urement of the spectral shift method used to retrieve V, is determined by measuring 
the spectral shift in a photonic crystal nanocavity with known V. This allowed us to 
give an accurate estimation of the modal volume, V = 0.2 um”, that is a value con- 
sistent both with the spatial extension retrieved by SNOM imaging of Fig. (6.5 c) and 
by statistical estimation based on the decay length reported in Fig. (6.4 c). 

To control the modes over a broader spectral range, we exploit the post-fabrica- 
tion technique based on laser-assisted micro-oxidation through SNOM probe pre- 
sented in Sec. (5.2.2). This technique permanently and locally modifies the dielectric 
environment by reducing the effective GaAs membrane thickness and increasing the 
effective pore diameter, thus shifting the localized photonic modes towards lower 
wavelengths. In order to demonstrate the ability to engineer the spectral and spatial 
properties of disorder based modes we carefully choose a set of three modes that are 
spatially close but detuned spectrally. Figures (6.6 a-c) show the PL intensity distri- 
bution of the three modes (labeled A, B and C) before starting the laser assisted oxi- 
dation process. Each mode is characterized by a main intensity peak surrounded by 
speckles with lower intensity. The spatial regions highlighted by the white (red) dotted 
lines in Fig. (6.6 a-c), identify the spatial overlap between A and B modes (B and C 
modes). The overlaps indicate a stronger expected interaction between A and B rather 
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than between B and C. The laser exposition is made in correspondence to the maxi- 
mum intensity of mode B (white circle). Figures (6.6 e and g) show the evolution for 
increasing laser exposure time of the normalized PL spectra averaged over the B-C 
and A-B overlap regions, respectively. The lower PL intensity in BC with respect to 
AB (demonstrated by the larger signal to noise ratio in the normalized spectra) reflects 
the smaller spatial overlap between the modes. 

Figure (6.6 f) shows the evolution of the resonant wavelengths of modes A, B and 
C as a function of the laser exposure time. Let us first consider the variation of mode 
B with respect to mode C. The wavelength of mode B blue-shifts monotonically 
whereas that of mode C stays at the unperturbed value within the experimental uncer- 
tainty. After about 310 min of exposure, the spectral overlap between the two modes 
is maximized. By further oxidizing the structure, the two resonances spectrally cross 
until they go out of resonance at At = 465 min. Therefore, the coupling strength mode 
B and C results much smaller than the broadening of the individual modes. The two 
modes are likely to be in the weak coupling regime. One possible way to get a greater 
interaction strength is to choose two modes with a larger spatial overlap, as for the 
case of mode A and B. The time evolution of the A and B wavelengths reported in 
Fig. (6.6 f) shows that the local oxidation undoubtedly induces a mode anti-crossing. 
The position of the anti-crossing occurrence is highlighted by reporting the wave- 
length splitting Q,g as a function of the exposure time, as shown in Fig. (6.6 d). Ini- 
tially Q4g decreases, reaching the minimum value around At = 170 min. Further- 
more, the amplitudes of the two interacting modes reach the same value at the time of 
minimum splitting, as highlighted by the the evolution of the PL spectra in Fig. (6.6 g). 
The black solid line in Fig. (6.6 d) is the result of the fit of the experimental splitting 
with the expected strong coupling behavior [21]. Following Eq. (5.10), the minimum 
splitting evaluated by the fit of the anticrossing experimental data directly gives the 
mode interaction strength. In fact, we get that Q4g min = 2g. We find that O42 min = 
(3.55 + 0.15) nm, that is larger than the mode spectral FWHM, equal to (3.00 + 
0.10) nm, and thus the strong coupling condition is satisfied. 

The observation of both resonant wavelength crossing and anticrossing behaviors 
represents the first step towards the creation of chains of hybridized localized ran- 
dombased modes, better known as necklace states, observed in 1D-systems and pre- 
dicted also for higher dimensions [22, 23]. 

In conclusion, we proved how it is possible to obtain and optimize strongly con- 
fined optical modes in two-dimensional disordered structures. Moreover, the near- 
field imaging of the localized photonic modes allowed us to statistically evaluate the 
degree of light confinement. The investigated two-dimensional disordered structures 
are particularly suited for improving the absorption efficiency of thin-film solar cells 
and can provide a new approach for high-extraction-efficiency in light-emitting di- 
odes. In fact, they sustain an efficient light trapping due to coherent optical effects 
along with broadband and wide-angle properties due to disorder. More strikingly, in 
this work we demonstrated the possibility to largely engineer the spectrum and the 
interaction of random modes by a fine-tuning of the disordered system in post-fabri- 
cation processes. Our achievements constitute an important step towards the control 
of disorder-based photonic resonators and traps. Therefore, they could be of great ad- 
vantage for developing random laser devices. Furthermore, the presented findings 
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may open new routes for addressing mesoscopic transport phenomena in random me- 
dia. As an example, the control of the optical confinement and the coupling between 
modes could be used to exploit two-dimensional necklace states in designing open 
transmission channels in disordered structures. 
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Figure 6.1: a)-b) SEM images of disordered samples with f = 0.35 and f = 0.13 and 
hole diameters hg = 215 nm and hg = 185 nm, respectively. The inset of a) shows 
the suspended photonic membrane clamped at its edges. 
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Figure 6.2: a)-b) SNOM PL intensity maps for samples with f = 0.35 and f = 0.13, 
respectively, evaluated at the wavelength of the highest intensity peak of each scan. 
c)-d) Near-field spectra of f = 0.35 and f = 0.13 detected in the position corre- 
sponding to the maximum value of the images shown in a) and b), respectively. The 
reported intensities as a function of position and as a function of wavelength are nor- 
malized to the average PL of each sample. 
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Figure 6.3: a)-b) Probability density function (P) of the PL measured intensities for 
f = 0.00, 0.13, 0.20, 0.26, 0.35 and hg = 200 nm, presented on linear and semi-log- 
arithmic scales, respectively. The PL intensities Ip, are normalized to the PL average 
value of the respective sample (u = (Ip,)). For the linear scale a) we also show the 
distribution for the slab without any holes (f = 0) divided by a factor of 10 that ex- 
hibits the minor fluctuations only due to the InGaAs quantum dots. c)-d) Normalized 
variance (02/1?) of the probability density function distributions as a function of f 
(with hg = 200 nm) and as a function of hd (with f = 0.35). The error bars of a? /u? 
have been evaluated by calculating this quantity for different wavelengths ranges, 
whereas for the filling fractions the uncertainty comes from the dispersion of the hole 
diameter sizes evaluated using SEM images. 
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Figure 6.4: a)-b) Average spectral FWHM ofthe peaks as a function of the penetration 
depth inside the sample (black line) for filling fractions f = 0.35 and f = 0.13, re- 
spectively, superimposed on a SEM image of the samples. The border of the samples 
is located at position x = 0. The red dotted line is an exponential fit to the data. c) 
Average decay lengths of the photonic modes versus filling fraction, evaluated by the 
exponential fits. The errors are evaluated by performing the above analysis along the 
four edges of each sample. 
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Figure 6.5: a) PL amplitude map of a localized mode in f = 0.35 and hg = 200 nm 
disordered sample, whose SEM image is superimposed. The amplitude is obtained by 
a Lorentzian fit of the collected spectra. b) Normalized spectra of the mode shown 
in a) evaluated close to the maximum amplitude (A) and at a point showing a low PL 
signal (B). The mode shift is caused by the presence of the dielectric SNOM tip in the 
near-field of the mode and its three-dimensional distribution is reported in the inset. 
c) Two-dimensional spectral shift map of the mode reported in a)-b) with superim- 
posed the topography map acquired during the SNOM scan. 
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Figure 6.6: a)-c) PL intensities of modes A, B and C before starting the laser-assisted 
local oxidation, respectively. The white circle indicates the position of the SNOM tip 
performing the laser-induced oxidation. The the white (red) dotted-lines regions iden- 
tify the spatial overlap between the modes A and B (B and C). The resonant wave- 
length of mode C is considered as the reference value defined Ap. d) Separation be- 
tween the resonant wavelengths of A and B (AB) as a function of the exposure time 
At showing an anticrossing behavior. The black solid line shows the fit of the exper- 
imental splitting with the expected trend for strong coupling. e) and g) Evolution of 
the normalized PL spectra mediated within the overlap regions B-C and A-B, respec- 
tively, for increasing laser exposure time. The average intensity of the modes within 
these regions is typically 10 — 15 times smaller than the maximum value of each 
mode. f) Resonant wavelength of the three modes evaluated in the overlap regions as 
a function of At. The error bars are provided by the two-peak Lorentzian fits of each 
spectrum. 
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Conclusions 


In conclusion, in this thesis we expanded the horizons of nano-photonics by de- 
veloping novel near-field probes, along with groundbreaking imaging techniques and 
post-fabrication methods. The presented experimental results are carefully supported 
by theory and numerical calculations. 


The introduction of the campanile tip nanoantenna on top of an optical fibre rep- 
resents a new paradigm for scanning near-field probes and, more in general, for sub- 
wavelength investigation. We demonstrated the strong field confinement and en- 
hancement at the campanile tip nanogap apex, which gives a very efficient throughput 
and an outstanding sub-wavelength spatial resolution, by investigating light emitting 
InP nanowires. The novel properties of the campanile tip are exalted by the perturba- 
tion imaging technique applied to photonic crystal nanocavities. In fact, a high-fidelity 
simultaneous detection of the electric and the magnetic field components of the local- 
ized light states was clearly observed. By means of the campanile tip perturbation 
imaging we obtained a comparable strength between the observed red and the blue 
spectral shifts, which are the fingerprints of the electric and magnetic interaction, re- 
spectively. In addiction we gained a dramatic improvement of the near-field imaging 
spatial resolution down to 2/30. The campanile near-field probe architecture could 
therefore allow the nanoscale investigation of light behaviour in optical metamaterials 
with an unprecedented resolution, thus triggering the quest for novel sensing applica- 
tions and field enhancement effects. 


It this work we also defined a new near-field imaging technique. The developed 
Fano phase-sensitive near-field investigation could revolutionize the sub-wavelength 
imaging research field for two reasons. Firstly, the optical imaging method is appli- 
cable to every kind of nano-resonator which exhibits sharp resonances, even realized 
on optically non active materials, such as silicon based as well as glass, polymer or 
metallic materials. Secondly, exploiting a detailed analysis of the modification of the 
Fano lineshapes, the technique allows the phase retrieval of the localized modes, since 
it is based on an integrated optical interferometer. Therefore, our developments opens 
up new strategies to investigate with a deep sub-wavelength spatial resolution the 
electric field intensity and the phase distribution of modes localized in a wealth class 
of nanophotonic and nanoplasmonic resonators as it does not need any special design 
or challenging fabrication steps and it can also image different polarization compo- 
nents. 
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Concerning the topic of photonic molecules, we characterized the coupling 
strength and the symmetry character of the supermodes in coupled nano-resonators 
with different alignment. The notably systems of two and of three closely aligned 
photonic crystal nanocavities were investigated. We found that the effective mode 
interaction is not limited to the nearest-neighbour coupling. 

We also developed a laser-assisted oxidation technique that enables a local post- 
growth modification of the underlying photonic crystal dielectric environment. Such 
post-fabrication technique was exploited to adjust the resonant wavelength of a given 
nano-resonator to a target value, and, more significantly, in combination with a local 
nano-infiltration procedure, was used to increase or decrease on demand the coupling 
strength in photonic molecules. It also allowed to probe the antibonding character of 
the K-aligned molecule ground state by means of an accurate minimization of the 
disorder induced detuning. The coupling between adjacent optical nano-resonators is 
a topic of vast attraction both for fundamental physics that exploits cavity quantum 
electrodynamics and for high density photonic integrated circuits applications. Our 
findings could open the way to exploit exotic configurations of coupled nanocavities 
to tailor the mode spatial distribution or the group velocity in CROW of any length 
by engineering the dielectric properties of adjacent resonators. 


In the last part of the thesis we investigated disordered photonic systems by scan- 
ning near-field microscopy. The presented analysis led to an optimization of the struc- 
ture parameters in order to achieve the highest degree of light localization. Moreover, 
we were able to fine tune the observed photonic modes and to probe the effective 
mode interaction between them by employing the laser-driven nano-oxidation tech- 
nique. Our achievements constitute an important step towards the control of disorder- 
based photonic resonators and traps. Therefore, they could be of large interest for 
random laser applications. Furthermore, they could allow for addressing mesoscopic 
transport phenomena in random media. In fact, the control of the coupling between 
photonic modes could be used to create two-dimensional necklace states and design 
open transmission channels in disordered structures. 

Finally, all the presented results put the basis for an unprecedented imaging in- 
vestigation of light behaviour in optical nano-resonators. Therefore, they would pos- 
sibly pave the way for novel interactions and devices which exploit the strong cou- 
pling regime between single light emitter and localized optical modes or even between 
complex interacting light states. 
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Appendix: 


Coupling strength calculation in photonic molecule 


In order to give an analytical expression of the coupling strength (g) between two 
nanocavities that constitute a photonic molecule, in Sec. (5.1.1) we used the result of 
tight-binding theory: 


JAe(MEo@)}-Eo(r+D)dr 
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(A.1) 


In this Appendix we will discuss the exact expression of the interaction between 
two optical resonators that results from Maxwell’s equations evaluated in a photonic 
crystals molecule. The following analysis will also strengthen the analogy between 
photonic molecules and coupled harmonic oscillators employed in Sec. (5.1.1). 

The usual employed hypothesis are the absence of e.m. sources, homogeneous 
and neutral linear dielectric media (with a magnetic permeability 4, = 1) so that the 
dielectric constant (e = n?) is a scalar depending only to the spatial coordinate. After 
factorizing the temporal evolution, it follows the Master equation (in c.g.s units) can 
be equivalently expressed both for the magnetic and for the electric field: 


fv ^ 5 VA} H(r) = (e) H(r) (A.2) 


(WAV AEC) = (2) «MEO (A3) 


The frequencies related to any of photonic problem for which e(r) is defined, can 
be found by solving the eigenvalues of the hermitian operator fv A se V A} But this 
is in general a difficult task. Therefore, for the photonic molecule case, we take a 
diferent approach. We choose the Master Eq. (A.3) to exploit the linearity of the de- 
rivatives. Then, the magnetic field can be evaluated as: 


i 
H(r)= -uV ^E) (A.4) 
while the electric filed has still to satisfy the scalar Maxwell equation: 


V-e(r)E(r)=0 (A.5) 
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Then, considering for instance the K-aligned photonic molecule reported in 
Fig. (A.1 c), we apply Eq. (A.3) to the single nanocavities, whose dielectric environ- 
ment is shown in Fig. (A.1 a-b), respectively, and to the photonic molecule as well: 


{VAVA}E;(r) = (ey e(NE(r) i=12; (4.6) 
VAVAE) = (2) a MEM AD 


where the frequency, the electric field and the dielectric constant are w;, E; and 
£i, respectively for the single nanocavities and w;, E; and €, for the molecule. There- 
fore Eq. (A.7) holds for the molecule case. 

Moreover, it is useful to express the overall dielectric constant as: 


(1) = & (r) + As, (r) = e2(7) + Aen (r) (A.8) 


where €, and €, are defined in Fig. (A.1 a-b). While Ag; corresponds to the four 
holes that form the D2 single nanocavity with refractive index equal to (Estap — Eqir)- 

Then, in the tight-binding hypothesis the molecular modes can be expressed as a 
linear combination of the single nanocavity modes, considering the coupling at a weak 
perturbation. This assumption gives: 


E;(r) = (,E\(r) + c2E2(r) (A.9) 


Therefore, putting together Eq. (A.6) and Eq. (A.7) thanks to the linearity of the 
derivatives and exploiting the relations (A.8) and (A.9), in the zero detuning approx- 
imation (w; = wz = Wo) we obtain the equation: 


Wo 2 Wt 2 
(=) (CEA + €C2 E>) = (=) [Cs + As) E1 + (ez + Ae2)c2E2|(A.10) 


c 


As often done in these kind of calculation we project (A.10) on the basis of the 
single cavity modes E; and E}, thus finding the matrix expression of the problem: 


IO: E WF )a4 1 = w?Aa,,| IO: = w? )a,2 = 4) Fa = 0(A.11) 
[(w3 = w? )azı = w? Aa, | [(w3 = w? )a2,2 = w? Aa, | c2 l 


where the coefficient a and Aa for i, j = 1, 2 are given by: 


aij = JEi@)-&(ME;(r) dr (A.12) 

Aa; j = | Ei(r) - Ag (r)E;(r) dr (A.13) 

In the framework of the perturbation method we can safely neglect some terms in 
(A.11). Moreover, by exploiting the fact that the two considered cavities are identical 


we therefore assume that: 


|we — w2| < wi (A.14) 
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Q42 = Q21; A414 = Q22; Ady, = Aaz,; 4a, = Aa, (A.15) 
Ay. © Ady2 K ayy (A.16) 


The condition (A.14) allows us to make the approximation w7Aa,, = w6Aa,, in 
the diagonal terms and w7Aa,2 = wéAa,. in the off-diagonal terms. The rela- 
tions (A.15) are due to the two considered identical cavities. The validity of (A.16) 
relies on the fact that the integral a4, involves the superposition of the same field, 
while the other terms in (A.16) are proportional to the superposition of different fields, 
which we assume to give a lower contribution. Therefore, (A.11) can be written as: 


@ =» Wi)a,, = w3Aa, | [-w3Aa, 2] (£) = 0(A.17) 
[-w3Aa,-] IO — wi)a,1 — w3Aa, | c2 i 


The solutions for w;, can be determined by solving: 


ie = WF )ay1 = wp Aa, 1] [-w3Aa,-] =0 (A.18) 
[-w3Aa, 2] [Co z WF )a11 z w3Aa, |] l 
which gives: 
w? = w? (1 = Min tha) (A.19) 


This finding means that the molecular modes are slighted spectrally shifted from 
the single cavity resonances. In particular, the term Aa, gives a shift due to the pres- 
ence of another cavity in the dielectric environment (Ag;) where the single cavity 
mode (E;) is present. Therefore, it is as a self-sustained perturbation. On the other 
hand, the term Aa; accounts for an effective interaction between the two cavities. In 
fact, as stated in the expression (A.13) it depends on the overlap integral between E, 
and E>. 

In order to add more clarity, we can also define the parameters which allow us to 
write a matrix equation for the photonic molecule problem that is equivalent to that 
of two classical coupled harmonic oscillators: 


P = wi (1 - Acs) (A.20) 
1,1 
g = "i (A.21) 
in fact, (A.17) becomes: 
p-o? ge L 
( n yg (e) =0 (A.22) 
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Coupled harmonic oscillators 


The motion equations of the system of two coupled classical harmonic oscillators, 

as the one reported in Fig. (A.2), are: 
ï + wax, + 97(x, — x2) = 0 (A.23) 
Xo + wkx, + g?(X2 — x) = 0 (A.24) 


where x, and x; are the position of the two masses, w3 = k/m and g? = kj/m. 
Then, by looking for wave solutions of the form x; = cje ‘©* we obtain in the basis 
of the two single oscillators: 


2 pg? 2 
ww +g =g Gay 
(Ge 


Finally, the Eq. (A.25) has the same structure and solutions of Eg. (A.22) if we 
define B = wê + g?. 
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Figure A.1: a) Dielectric environment £4 of the single left-side PCN. b) Dielectric 
environment £, of the single right-side PCN. b) Dielectric environment €; of the re- 


sulting photonic molecule. 


Figure A.2: Schematics of two classical coupled one-dimensional harmonic os- 
cillators made by two equal masses and springs, which interact through a third 


spring k, 2. 
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